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THE  DETERMINATION  OF  RATE-LIMITING  STEPS  DURING  SOOT  FORMATION 

Annual  Report 

I.  Summary 

A  single-pulse  shock  tube  has  been  used  to  examine  toluene  pyrolysis  and  the  rich  oxidation 
of  benzene  {<f>  >  7.5)  diluted  in  argon  over  the  temperature  range  of  1200  to  2000K  and  at  total 
pressures  of  ten  to  thirteen  atmospheres.  Dwell  times  were  about  500  microseconds.  Collected 
gas  samples  were  analyzed  using  gas  chromatography  for  hydrogen,  carbon  oxides  and  Ci-  to  Ch- 
hydrocarbons.  Low  molecular  weight  products  during  benzene  oxidation  include  cyclopentadiene 
and  vinylacetylene  and  support  literature  proposals  for  oxidation  of  benzene.  High  molecular  weight 
products  are  dominated  by  species  containing  mixtures  of  five-  and  six-membered  rings.  During 
toluene  pyrolysis  very  rapid  production  of  polyaromatic  hydrocarbons  occurs  between  1400  and 
1450K,  consistent  with  temperatures  at  which  particle  inception  occurs  in  diffusion  flames. 

Simplified  modeling  techniques  of  soot  formation  have  been  explored.  The  most  successful 
model  is  one  based  on  the  MAEROS  code.  Using  experimentally  determined  benzene  production 
rates  from  a  premixed  ethylene,  laminar  flame  (Harris,  Weiner  and  Blint)  and  known  rates  for 
carbon  (as  acetylene)  addition  to  particles,  production  rates  for  soot  are  predicted  accurately. 

II.  Introduction 

There  is  increasing  experimental  evidence  that  important  rate-limiting  steps  to  soot  formation 
are  the  production  and  growth  of  aromatic  rings.  Detailed  modeling  and  comparison  to  experimental 
data  has  led  to  a  good  understanding  of  mechanisms  and  rates  for  the  production  of  benzene  (and 
phenyl  radical).  Existing  mechanisms  for  the  formation  of  naphthalene  (plus  related  species)  and 
higher  order  aromatics  are  principally  thermochemical  estimates  and  no  experimental  verification 
of  these  steps  (mechanisms  or  rate  constants)  are  available.  Another  large  uncertainty  in  soot 
formation  mechanisms  is  the  effect  of  oxidation.  Under  many  conditions,  growth  of  polyaromatic 
hydrocarbons  (PAH’s)  may  be  quite  fast  but  is  counterbalanced  by  oxidation.  Therefore,  the  net 
production  rate  of  PAH’s  and  therefore  of  soot  is  only  a  small  fraction  of  the  “gross”  production 
rate  of  PAH’s.  Oxidation  processes  (mechanisms  and  rates)  of  aromatic  hydrocarbons  at  elevated 
temperatures  (>  1200K)  are  unfortunately  very  poorly  known. 

To  investigate  these  important  and  competitive  phenomena,  an  experimental  and  modeling 
program  is  being  performed.  The  experimental  phase  of  the  program  involves  the  use  of  a  single- 
pulse  shock  tube  to  thermally  stress  mixtures  of  gases,  which  are  then  quenched,  and  the  reaction 
products  are  analyzed.  This  program  is  focusing  on  detection  of  two-  and  three-ringed  products  in 
order  to  test  mechanisms  proposed  for  the  formation  of  these  species.  In  support  of  this  work  and 
under  corporate  sponsorship,  the  gas  sampling  system  has  been  modified  to  collect  and  detect  high 
molecular  weight  species.  To  examine  the  competitive  oxidation  processes,  rich  oxidation  of  benzene 
has  been  investigated.  Furthermore,  in  order  to  help  confirm  detailed  chemical  kinetic  modeling, 
species  concentrations  wiM  be  meas'"’‘*d  in  real  time  by  '•oupling  a  time-of-flight  mass  spectrometer 
to  the  end  wall  of  the  shock  tube. 

Two  types  of  modeling  a^e  being  performed  as  part  of  this  contract.  Detailed  chemical  kinetic 
modeling  is  being  used  for  comparison  to  experimental  data  and  to  test  reaction  mechanisms  and 
rate  constants  for  formation  of  multi-ringed  species  and  for  oxidation  processes.  Secondly,  using 
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concepts  developed  from  this  work  and  from  the  literature,  simplified  models  are  being  examined 
which  may  be  useful  in  describing  the  formation  and  growth  of  soot.  Since  this  latter  work  is 
relatively  new  and  the  results  are  extremely  encouraging,  details  of  the  calculations  are  provided  in 
this  report. 

III.  Results 

A.  Modifications  to  Experimental  Facilities 

A  primary  goal  of  this  program  is  to  collect  experimental  data  on  the  production  of  high  molec¬ 
ular  weight  species  during  pyrolysis  and  “oxidative  pyrolysis”  of  hydrocarbons.  In  order  to  collect 
this  data,  the  sampling  system  was  modified.  The  system  was  modified  under  corporate  sponsor¬ 
ship  and  several  designs  and  calibration  procedures  were  tested.  Schematics  of  the  original  and 
final  sampling  system  are  compared  in  Figures  la  and  lb.  In  the  original  system,  lines  were  heated 
to  80°C,  whereas  for  the  modified  system,  all  lines  were  heated  to  150°C,  and  better  monitoring 
procedures  were  established  to  avoid  “cold”  spots.  In  addition,  high  temperature  solenoids  were 
installed  in  the  modified  system.  An  important  change  was  moving  the  sampling  collection  volume 
to  after,  rather  than  before,  the  gas  sampling  loops.  This  was  required  in  order  to  minimize  the 
number  of  valves  between  the  shock  tube  and  the  sampling  loops,  since  each  valve  resulted  in  some 
loss  of  the  high  molecular  weight  products.  This  design  is  not  ideal  since  it  does  not  allow  for 
mixing  any  non-uniformities  in  the  collected  sample.  Nevertheless,  comparison  of  data  obtained 
using  the  modified  system  compares  well  with  data  obtained  using  the  original  sampling  arrange¬ 
ment.  In  addition,  an  enhancement  in  the  calibration  technique  wm  developed.  All  species  with 
molecular  weights  of  about  100  g/mole  and  below  could  be  adequately  calibrated  using  premixed 
gaseous  samples  in  bottles  which  could  be  interchanged  using  quick  disconnect  fittings.  Accuracy 
of  analysis  for  these  species  is  approximately  ±  3  to  6  percent.  This  uncertainty  is  size  dependent 
with  the  larger  uncertainty  assigned  to  the  heavier  species  (e.g.,  benzene,  toluene,  styrene,  etc.). 
For  yet  larger  species,  solid  samples  were  dissolved  in  toluene  and  then  injected  into  an  atmosphere 
of  argon  and  allowed  to  mix  at  high  temperature  in  a  mixing  volume  for  about  one-half  hour  (See 
Fig.  lb).  Subsequently,  the  gas  was  analyzed  and  calibration  factors  were  determined  relative  to 
the  toluene  response.  Uncertainties  were  approximately  ±  10  to  30  percent,  again  depending  on  the 
size  of  the  molecule.  The  largest  species  for  which  a  calibration  was  obtained  was  phenanthrene. 
Attempts  to  observe  pyrene  failed.  The  very  large  uncertainties  for  the  heaviest  species  suggests 
that  profiles  of  these  species  are  semi-quantitative. 

A  second  goal  of  this  program  is  to  observe  species  concentrations  in  real-time  rather  than 
just  analyzing  end  geises  from  the  single-pulse  shock  tube.  In  order  to  perform  this  research,  a 
time-of-flight  mass  spectrometer  (CVC)  will  be  coupled  to  the  end  wall  of  the  shock  tube,  and  data 
collected  using  a  high  speed  data  acquisition  system  (LeCroy).  The  development  and  installation 
of  this  equipment  is  behind  schedule  due  to  ordering  problems,  and  further  delayed  since  the  mass 
spectrometer  has  not  yet  been  delivered,  although  it  is  now  well  past  the  manufacturer’s  delivery 
date.  This  work  is  being  performed  under  corporate  sponsorship. 

B.  Experimental  Results 

Since  the  end  of  the  previous  AFOSR  contract  in  November  1987,  twenty  separate  mixtures 
have  been  shock  heated  and  products  analyzed  by  gas  chromatography.  The  mixtures  are  listed  in 
Table  I.  Six  of  the  aiixtures  were  examined  as  part  of  the  first  year  of  this  AFOSR  program;  the 
remainder  were  tested  under  corporate  sponsorship.  For  these  tests,  dwell  times  were  typically  500 
microseconds  and  total  pressures  (balance  argon)  were  9  to  13  atmospheres. 

In  the  following  paragraphs,  results  on  toluene  pyrolysis  and  benzene  oxidation  are  discussed. 
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Figure  la.  Original  Sampling  System 
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Table  1 


Series  of  Experiments  Completed  During  the  First  Year  of  AFOSR  Program 


Reactants 

Initial  Concentrations 

Toluene 

1% 

Benzene 

1.1% 

Benzene/oxygen 

1.1/0.22% 

Benzene/oxygen 

0.11/0.022% 

Benzene/oxygen 

1. 1/1.1% 

Dicyclopentadiene 

0.125% 

*  Acety  lene /ethene 

1. 1/2.2% 

*  A  cety  lene /ethene 

0.125/0.25% 

*Toluene  /  methanol 

1.0/0.15% 

*Toluene /hydrogen 

1.0/3. 1% 

*Methylcyclohexane 

1% 

*Methanol 

4% 

*Natural  Gas/oxygen 

4/1% 

*Methane/oxygen 

4/1% 

*Cyclopentadiene 

1.4% 

*Cyclopentadiene 

0.1% 

*  Performed  under  corporate  sponsorship 


The  data  on  dicyclopentadiene  has  not  been  analyzed  due  to  the  complexity  of  the  gas  chro¬ 
matogram.  Corporate  sponsored  work  subsequently  demonstrated  that  dicyclopentadiene  could 
be  successfully  decomposed  (by  refluxing)  into  cyclopentadiene,  which  is  then  mixed  with  argon 
and  shock  heated  in  the  single-pulse  shock  tube.  Work  on  cyclopentadiene  will  continue  in  the 
second  year  of  the  AFOSR  program. 

1.  Toluene  Pyrolysis 

One  of  the  first  mixtures  analyzed  for  this  program  was  1%  toluene,  since  several  previous  runs 
with  this  mixture  have  been  performed  and  provide  a  reference  set  of  data.  The  newly  obtained 
data  is  in  excellent  agreement  with  the  results  presented  previously  (Ref.  1)  for  the  lower  molecular 
weight  compounds  (<  128  g/mole).  In  addition,  because  of  the  enhanced  sampling  system,  profiles 
on  a  variety  of  PAH’s  were  also  obtained  in  this  recent  series  of  experiments.  Some  of  the  aromatic 
compounds  produced  from  pyrolysis  of  1%  toluene  are  shown  in  Figure  2.  Profiles  of  toluene, 
benzene,  indene,  and  naphthalene  shown  in  Figure  2  agree  well  with  the  earlier  data.  The  remaining 
profiles  are  new  information.  The  low  temperature  formation  of  bibenzyl  is  from  recombination 
of  benzyl  radicals;  although,  due  to  uncertainty  in  the  thermochemistry  of  benzyl,  it  is  unknown 
whether  the  bibenzyl  could  be  produced  at  high  temperatures  or  in  the  quenching  waves.  Hopefully, 
detailed  modeling  will  assist  in  deciphering  this  data.  The  other  products  shown  in  Figure  2 
represent  a  relatively  small  portion  of  the  total  number  of  high  molecular  weight  species,  although 
the  mass  contained  in  these  selected  species  represents  about  75%  or  more  of  the  total  mass  of  the 
species  observed.  Possible  reactions  describing  the  production  of  selected  species  are  proposed  in 
Figure  3.  Analysis  of  this  data  and  detailed  modeling  is  continuing.  The  result  that  PAH’s  are 
produced  very  rapidly  at  1400  to  1500K  is  consistent  with  concepts  developed  from  diffusion  flames 
which  indicate  inception  occurs  at  1300  to  1500K. 

2.  Benzene  Oxidation 

The  oxidation  of  benzene  was  also  examined  during  this  past  year.  Prior  to  the  oxidation  runs, 
a  pyrolysis  experiment  was  performed.  However,  based  on  the  product  distribution  and  comparison 
with  previously  obtained  data,  it  is  believed  that  an  impurity  of  a  small  amount  of  oxygen  (~  1000 
ppm)  affected  the  recent  “pyrolysis”  experiments.  Benzene  decay  and  light  product  formation 
(methane,  acetylene,  diacetylene,  cyclopentadiene,  and  vinylacetylene)  are  shown  in  Figures  4  -  6 
for  the  two  oxidation  runs  at  high  benzene  concentrations  and  for  the  previous  benzene  pyrolysis 
experiments.  In  these  figures,  the  measured  concentrations  of  reactants  rather  than  the  percent 
calculated  from  partial  pressure  mixing  (Table  1)  are  listed.  The  effect  of  oxygen  clearly  is  to  cause 
decay  of  benzene  at  lower  temperatures  as  well  as  to  enhance  the  low  temperature  production  of 
the  dominant  product,  acetylene.  In  addition,  oxygen  causes  the  low  temperature  formation  of 
cyclopentadiene  and  vinylacetylene,  which  strongly  supports  the  mechanistic  arguments  proposed 
by  Venkat  et  al  (Ref.  2)  for  benzene  oxidation.  Oxygen  only  slightly  affects  the  production  of 
diacetylene.  This  species  is  a  byproduct  of  the  thermal  decomposition  of  phenyl  radical  which  is 
slow  due  to  its  high  activation  energy  (~  70  kcal/mole).  The  slight  enhancement  in  the  production  of 
diacetylene  with  increasing  oxygen  presumably  is  due  to  an  increase  in  phenyl  radical  concentration 
(via  radical  attack  on  benzene). 

Dominant  polyaromatic  hydrocarbons  produced  during  the  rich  oxidation  of  benzene  (1.25% 
benzene,  0.3%  oxygen)  are  shown  in  Figure  7.  A  large  number  of  other  high  molecular  weight 
products  were  also  produced,  although  many  of  these  are  unidentified.  It  is  estimated,  however, 
that  the  total  mass  contained  in  unidentified  species  is  less  than  30%  of  the  mass  for  all  the  high 
molecular  weight  compounds.  The  dominant  PAH  products  are  similar  to  those  observed  during 
pyrolysis  of  toluene.  In  both  cases,  a  predominance  of  both  five-  and  six-membered  rings  is  ob¬ 
served.  The  presence  of  species  containing  five-membered  rings  is  not  surprising,  based  on  previous 
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Formation  of  Selected  Aromatic  Products 
SPST  Pyrolysis  of  1%  Toluene 


Figure  2 
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Benzene  Pyrolysis  11500  PPM 
Liqht  Products 


Figure  6 
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studies  in  flames;  however,  existing  models  for  PAH  growth  and  soot  formation  neglect  the  potential 
importance  of  these  species.  Profiles  and  identities  of  PAH  species  during  the  oxygen-contaminated 
experiment  on  benzene  pyrolysis  were  similar  to  those  in  Figure  7,  although  they  were  shifted  to 
higher  temperatures.  The  profiles  of  PAH’s  for  the  less  rich  oxidation  (1%  oxygen)  were  severely 
depressed,  although  a  carbon  balance  (which  includes  formation  of  carbon  monoxide  and  carbon 
dioxide)  indicates  significant  loss  of  mass.  The  implication  is  that  with  the  larger  amounts  of  oxy¬ 
gen,  the  growth  of  PAH’s  is  significantly  accelerated.  These  results  are  considered  preliminary  and 
need  to  be  verified  in  future  experiments. 

C.  Model  of  Soot  Particle  Formation 

1.  Model  Soot  Inception  Kinetics /Nucleation  Calculations 

Soot  inception  kinetics  have  been  modeled  by  considering  a  mechanism  of  stepwise  acetylene 
addition  of  the  form 


Ai  -f-  Hj  Ai  -h  H 

(0) 

Ai  +  CjHz  AiCzH  +  H 

(1) 

AiCjH -t- H  AiCzH  +  H2 

(2) 

AiC2H  -f  C2H2  — ^  Aj+i 

(3) 

Here  Aj  is  a  closed  ring  aromatic  species  (Ai=benzene,  A2=naphthalene,  Ai=phenyl,  etc.);  the 
mass  of  Ai  is  78-f(i-l)50  a.m.u. 

Setting  the  time  derivatives  of  radical,  intermediate  concentrations  equal  to  zero  (steady  state 
approximation)  results  in  the  following  set  of  coupled  rate  equations  for  the  concentrations  of  the 
closed  ring  species. 


d[Ai] 

_  -[All 

dt 

T, 

d(A2] 

_  “(Azl 

dt 

Ti 

dfAs) 

_  -[As] 

dt 

Ti 

dlAi] 

_  -lAil 

dt 

Ti 

+  Q 

,  (Atl 
fz 

+  +  ‘■[Ai] 


J-it 


where  Q  is  the  benzene  source  rate  and 
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[All  =  Q  Ti(l 

|A,|  =  (ll  +  2,)  (A,|  -  r  (li  .  lA,)  -  Q  (^) 

It  is  convenient  to  resort  to  numerical  integration  of  the  rate  equations,  because  the  analytic 
solutions  become  complicated  and  we  want  eventually  to  consider  the  effects  of  agglomeration 
reactions  of  the  form 


Ai  +  Aj-^(Ai)(Aj)  =  Aij  (large  i,j) 

The  Aij  can  participate  in  the  stepwise  acetylene  addition  mechanism  and  agglomerate  to  form 
higher  order  aggregates,  and  the  kinetic  master  equations  have  been  suitably  modified  to  include 
terms  for  these  additional  processes.  The  following  provisional  rate  coefficient  data  were  employed 
in  our  simulations: 

ko  =  5.1  X  10**  g-5020/T  ^nioles/cc/sec) 
k-o  =  2.51  X  10**  e-808j/T 

ki  =  10**  e-*'**^'''*’ 
k-i  =  7.6  X  10*^  e-3030/T 
k2  =  2.51  X  10*^  e-808i/T 
k_2  =  1.51  X  10*® 
ka  =  2  X  10**  e-20»3/T 
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The  following  species  concentration  parameters  were  estimated  from  the  Harris  papers  (Ref.  3-5): 


[C2H2]  =  2.3  X  10  ^  moles/cc 
[H2]  =  5.2  X  10"^ 

[H]  =  7.4  X  10"® 

Q  =  4.9  X  10~^  moles/cc/sec 


For  these  parameters  at  T=1600K,  Tj  ~  1.3  millisec,  and  T1/T2  —  0.3.  A  simulation  of  the  time 
evolution  of  the  closed  ring  concentrations  without  aggregation  is  shown  in  Figure  8.  It  is  seen  that 
the  time  to  achieve  steady  state  for  the  N-th  closed  ring  species  is  approximately  NTi.  Thus  the 
time  required  to  form  a  1000  a.m.u.  particle  for  which  N  ~  20  would  be  20  -  30  milliseconds,  which 
is  about  an  order  of  magnitude  too  slow  for  the  process  of  soot  inception.  The  steady  state  ratio  of 
lAi]/lAi_i](i<2)  is(Ti/T2  -I-  r);  thus,  the  steady  state  concentration  of  the  i-th  closed  ring  is 


[Ai]  ~  QTi 


i-2 


When  Q  is  assumed  to  have  a  reasonable  shut-off  time  of  one  millisecond,  the  time  histories  of  the 
closed  ring  species  are  as  shown  in  Figure  9.  Note  the  qualitative  resemblance  to  the  experimental 
measurements  of  Bittner  and  Howard  (Ref.  6).  While  some  uncertainty  attaches  to  the  rate 
coefficients  and  species  concentrations  which  go  into  this  calculation,  the  indication  is  that  our 
simple  stepwise  acetylene  addition  mechanism  is  too  slow  to  explain  observed  soot  inception  rates 
in  the  absence  of  agglomerative  or  other  mechanisms.  Implicit  in  the  foregoing  analysis  has  been 
the  assumption  that  the  rate  coefficients  ko  -  ks  do  not  scale  with  the  order  of  Ai.  In  fact,  there 
are  qualitative  arguments  one  can  make  that  these  rates  should  increase  monotonically  with  i.  As 
i  increases,  the  number  of  sites  available  for  H-atom  abstraction  and  acetylene  addition  increases, 
and  so  should  the  rate  coefficient.  Preliminary  estimates  are  that  the  scaling  could  be  fairly  strong, 
leading  to  significant  shortening  of  the  time  required  to  form  large  mass  species.  This  question  will 
be  further  addressed  in  the  future;  the  implication  is  that  the  rate-limiting  step  will  prove  to  be  the 
formation  rate  of  the  first  few  closed  ring  species. 

Reactions  of  the  form 

Ai  4-  Aj  ~  AiAj  H 

Ai  -f"  Aj  ■’  ■  AiAj 


could  also  be  important.  The  sequence 


Ai  -f-  A|  5=^  AiAi  -1-  H 
AiAi  -f  H  ^  AiAi  -f  H2 
AiAi  +  C2H2  ^  As  4-  H 


for  example,  could  provide  an  alternative  fast  mechanism  for  the  production  of  A3.  These  reactions 
are  presently  being  incorporated  into  the  steady  state  analysis  to  see  if  they  lead  to  accelerated 
production  of  A3. 

When  physical  agglomeration  processes  are  allowed,  our  preliminary  calculations  show  that  the 
production  of  higher  mass  species  can  be  significantly  accelerated,  but  the  results  are  sensitive  to 
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Evolution  of  Closed  Ring  Aromatics 
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the  model  for  the  sticking  probabilities.  The  question  of  scaling  of  these  probabilities  with  particle 
size  is  one  that  will  have  to  be  addressed  with  some  care  if  the  calculations  are  to  be  meaningful. 
The  results  of  more  extensive  agglomeration  calculations  will  be  presented  in  a  future  report. 

Finally,  we  also  considered  the  stepwise  mechanism 

Ai  +  Hj  A;  -|-  H 
Ai  +  C2H2  ^  AiC2H2 
AiC2H2  +  C2H2  ^  AiC4H4 
AiC4H4  Ai+i  +  H 

Given  reasonable  estimates  for  the  rate  coefficients  in  a  steady  state  analysis,  however,  the  charac¬ 
teristic  formation  times  are  much  too  long  for  this  to  be  the  correct  inception  mechanism. 

2.  Soot  Aerosol  Dynamics  Simulations 

After  the  inception  or  nucleation  stage,  soot  particles  continue  to  grow  as  isolated  spheroids, 
reaching  a  size  of  40-100  nm;  at  later  times,  these  spheroids  fuse  to  form  chainlike  clusters  or 
aggregates.  Prior  to  the  aggregation  phase,  the  growth  of  the  primary  spheroids  occurs  through 
surface  deposition  of  gas  phase  hydrocarbons  and  coalescing  collisions  in  which  colliding  spheroids 
stick  together  and  form  a  new  spheroid. 

In  order  to  model  the  growth  of  primary  soot  spheroids,  the  most  recent  version  of  an  aerosol 
dynamics  code,  MAEROS  (Ref.  7),  heis  been  appropriately  modified.  This  is  a  sectional  model  in 
which  the  particle  size  range  of  interest  is  divided  into  discrete  classes,  or  intervals.  A  set  of  master 
equations  for  the  rate  of  change  of  particle  number  density  in  each  size  class  is  then  integrated 
numerically,  accounting  for  the  processes  of  nucleation,  surface  growth,  and  coalescence.  Each  size 
class  has  a  source  term  representing  nucleation,  but  in  our  simulations  this  is  set  equal  to  zero  for 
all  but  the  smallest  size  class.  Surface  growth  is  assumed  to  occur  primarily  through  acetylene 
deposition,  using  the  Harris- Weiner  growth  rate  expression  (Ref.  3).  In  this  formulation,  the  rate 
of  mass  addition  to  a  particle  is  proportional  to  the  particle  surface  area,  and  is  linear  in  the 
acetylene  vapor  concentration.  The  effective  surface  growth  rate  given  by  Harris  and  Weiner  is 
time  dependent,  and  we  have  chosen  an  early,  or  young,  soot  value.  Provision  is  made  for  acetylene 
vapor  depletion,  and  oxidation  is  accounted  for  by  modifying  the  surface  growth  term  to  allow  for 
oxidative  mass  removal  at  a  rate  given  by  the  well-known  Nagle  and  Strickland-Constable  expression 
(Ref.  8).  The  latter  is  plotted  in  Figure  10.  It  can  be  seen  that  for  oxygen  concentrations  above 
roughly  5%,  and  temperatures  below  approximately  2000K,  the  effective  oxidation  rate  becomes 
essentially  independent  of  oxygen  concentration.  Coalescence  rate  coefficients  are  calculated  in 
the  MAEROS  program  using  the  Fuchs-Sutugin  expression  (Ref.  9),  which  spans  the  Knudsen 
number  range  from  free-molecule  to  continuum.  A  coalescence  sticking  probability  is  input  to  the 
computer  code.  Kennedy  (Ref.  10)  and  Harris  and  Kennedy  (Ref.  11)  have  shown  that  soot  particle 
collision  rates  can  be  significantly  enhanced  by  long-range  Van  der  Waals  forces.  These  forces 
between  electrically  neutral  dielectrics  have  their  origins  in  quantum  fluctuations  as  embodied  in 
the  Heisenberg  uncertainty  principle,  which  create  instantaneous  dipole  moments  in  the  colliding 
bodies.  We  have  written  a  computer  program  to  calculate  the  collision  rate  enhancements  arising 
from  these  forces.  The  program  is  based  on  the  Lifshitz(Ref.  12),  Langbein  (Ref.  13),  and  Marlow 
(Ref.  14)  theory,  and  is  similar  to  the  calculation  outlined  by  Kennedy  (Ref.  10).  All  dielectric 
constant  parameters  are  the  same  as  those  given  in  Reference  10.  Calculated  enhancement  factors 
for  soot  particles  in  the  free-molecule  regime  at  1600K  are  shown  in  Figure  11.  It  can  be  seen 
that  the  enhancements  are  a  function  mainly  of  the  radius  difference  between  particles,  and  for  the 
smallest  particles  reach  a  value  of  about  2.4. 
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Rate  Constant  for  Soot  Oxidation 
Naale  &  Strickland  Constable  Formula 


Figure  10 
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The  most  important  approximation  made  in  our  soot  aerosol  dynamics  simulations  has  been  to 
set  the  lower  limit  of  size  range  at  the  approximate  size  (.5nm  -  5 A)  of  a  benzene  molecule,  and  to 
extrapolate  the  Harris- Weiner  surface  growth  rate  expression  to  this  size  class.  The  nucleation  rate 
for  this  size  range  is  taken  to  be  the  local  benzene  production  rate,  and  the  nucleation  rate  is  set 
equal  to  zero  for  all  other  size  classes.  This  clearly  is  an  approximation,  and  a  more  careful  coupling 
of  soot  inception  kinetics  with  aerosol  dynamics  simulations  might  be  needed  in  the  future,  but,  as 
will  be  seen,  its  use  leads  to  some  encouraging  results. 

The  MAEROS  code  has  the  advantage  of  computational  efficiency;  for  the  size  range  .5  to  1000 
nm,  with  15  -  20  size  classes,  the  computation  times  are  on  the  order  of  VAX-minutes.  There  is, 
however,  some  sensitivity  of  the  results  to  number  of  size  classes;  a  change  from  20  to  15  changes  the 
calculated  numbers  by  about  10%.  Calculations  with  a  larger  number  of  size  intervals  are  planned, 
but  require  non-trivial  code  modifications.  It  may  ultimately  be  desirable  to  use  a  code  employing 
moving  sectional  boundaries,  because  of  reduced  numerical  diffusion  error;  we  are  planning  to  check 
our  results  by  having  F.  Gelbard  run  our  cases  on  an  unreleased  version  of  MAEROS  that  uses  the 
moving  sectional  grid  technique.  All  calculations  to  follow  have  been  performed  with  the  maximum 
20  size  classes  presently  allowed  in  the  program.  In  these  calculations,  the  benzene  source  rate  and 
acetylene  concentration  were  taken  as  before  from  S.  Harris  flame  conditions  (Ref.  3  -  5)  in  order 
to  attempt  to  reproduce  the  measured  soot  volume  fraction  profiles  reported  in  Reference  3. 

Figure  12  shows  the  calculated  evolution  of  the  soot  number  density  distribution  with  time.  It 
is  apparent  that  the  dynamics  model  predicts  the  appearance  of  10  -  100  nm  particles  on  a  time 
scale  of  about  10  milliseconds.  It  is  also  clear  that  after  about  four  milliseconds,  the  shape  of  the 
distribution  has  taken  on  an  invariant  form.  This  is  the  “self-preserving”  size  distribution  (Ref.  15), 
as  can  be  seen  in  Figure  13,  where  the  same  data  are  plotted  as  a  function  of  the  normalized  variables 
n'(v)v/N  and  v/v,  where  v  is  particle  volume,  v  is  average  particle  volume,  and  N  =  n'(v)dv.  A 
self-preserving  form  is  expected  for  coalescence-dominated  kinetics  (Ref.  15);  our  calculations  show 
that  it  can  still  be  expected  even  in  the  presence  of  nucleation  and  surface  growth. 

The  calculated  soot  volume  fraction  is  plotted  as  a  function  of  time  in  Figure  14  for  various 
values  of  the  coalescence  sticking  probability,  which  hM  been  assumed  here  to  be  size-independent. 
For  very  low  values  of  the  sticking  probability,  coalescence  does  not  proceed  rapidly,  and  the  relative 
surface  area  of  the  particles  remains  high.  Surface  deposition  proceeds  rapidly  when  there  are  large 
numbers  of  small  particles,  and  the  total  soot  volume  fraction  increases  to  large  values.  When 
sticking  probability  is  increased,  coalescence  becomes  important  and  a  smaller  number  of  larger 
particles  is  formed.  Coalescence  in  itself  conserves  volume  fraction,  but  the  reduced  relative  surface 
area  of  the  larger  particles  reduces  mass  addition  by  surface  growth.  Thus,  the  calculated  volume 
fraction  decreases  as  seen  with  increasing  sticking  probability.  Only  for  the  very  low  values  of 
sticking  probability  is  significant  depletion  of  the  acetylene  vapor  predicted.  It  is  interesting  to  note 
that  there  is  excellent  agreement  with  experiment  for  sticking  probabilities  in  the  expected  range 
1-2.  Both  the  magnitude  and  time  dependence  of  the  Reference  3  results,  which  were  obtained 
from  absorption  measurements,  are  seen  to  be  reproduced  well.  At  the  present  time,  the  success  of 
a  model  which  omits  the  fine  details  of  the  inception  kinetics  is  possibly  fortuitous,  but  certainly 
encouraging.  These  results  support  previous  work  indicating  that  the  rate-limiting  process  of  soot 
formation  is  the  formation  of  the  first  few  rings.  The  extrapolation  of  the  Harris- Weiner  surface 
growth  rate  to  the  smallest  molecule/particle  size  classes  appears  to  effectively  reproduce  the  effects 
of  more  detailed  inception  kinetics;  this  result  (supported  by  recent  calculations  of  Howard,  Harris, 
and  Frenklach)  might  form  the  basis  of  a  very  simple  and  powerful  model  for  soot  formation.  Further 
examination  of  these  questions  is  planned.  Another  point  of  interest  is  whether  the  aging  of  the 
acetylene  surface  deposition  rate  reported  in  Ref.  3  is  really  due  to  temperature  effects  and  an 
activation  energy  for  the  process. 
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Calculated  Evolution  of  Soot  Size  Distribution 
Coagulation  &  Surface  Growth  from  Maeros  Code 
Harris  Surface  Growth  Rate  and  Flame  Conditions 


Figure  12 
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Log(lO)  Particle  Diameter  (Angstroms) 


i  Class  Distributions  from  MAEROS  Code 
Surface  Growth  Rate  and  Flame  Conditions 


Normalized  Volume 


Calculated  Soot  Volume  Fraction 
Harris  Flame  Conditions 
Acetylene  Vapor  Depletion 


Figure  14 
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Calculated  Soot  Volume  Fraction 


Figure  15 
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Performing  the  soot  growth  calculations  with  oxygen  present  shows  that  even  modest  amounts 
of  oxygen  can  have  rather  drastic  effects  on  the  predicted  amount  of  soot  formed.  Figure  15  displays 
the  calculated  soot  volume  fraction  versus  oxygen  concentrations.  Given  the  Harris- Weiner  acetylene 
surface  growth  rate  and  the  Nagle  and  Strickland-Constable  surface  oxidation  rates,  it  appears  not  to 
be  possible  to  form  soot  for  oxygen  concentrations  beyond  roughly  0.2%  at  a  temperature  of  1600K. 
No  depletion  of  the  oxygen  by  the  oxidation  process  has  been  allowed  for  in  these  calculations, 
however. 


IV.  List  of  Publications 

A  paper  entitled  “The  Pyrolysis  of  Acetylene  Initiated  by  Acetone”  by  M.  B.  Colket,  H.  B. 
Palmer  and  D.  J.  Seery  has  been  published  in  Combustion  and  Flame,  Vol.  75,  pp.  343-366,  1989. 
This  work  was  initiated  under  contract  F49620-85-C-0012,  and  revisions  to  the  manuscript  were 
performed  under  this  contract.  See  Appendix  C  for  reprint. 

A  paper  entitled  “The  Pyrolysis  of  Acetylene  and  Vinylacetylene  in  a  Single-Pulse  Shock  Tube” 
by  M.  B.  Colket  was  published  in  the  Twenty- First  Symposium  (International)  on  Combustion,  The 
Combustion  Institute,  pp.  851-864,  1986.  See  Appendix  D  for  reprint.  Research  for  this  publication 
was  performed  under  contract  F49620-85-C-0012. 

A  manuscript  entitled  “Oxidative  Pyrolysis  of  Ethene  in  a  Single-Pulse  Shock  Tube”  by  M.  B. 
Colket  and  D.  J.  Seery  will  be  written  late  in  1989  and  submitted  to  the  Twenty-Third  Symposium 
(International)  on  Combustion,  Jan.  1990. 

A  manuscript  entitled  “A  Simplified  Model  for  the  Production  of  Soot  in  a  Premixed  Flame” 
by  R.  J.  Hall  and  M.  B.  Colket  will  be  written  in  1989  and  submitted  to  Combustion  Science  and 
Technology. 

V.  Meeting  Interactions  and  Presentations 

1.  A  round  table  discussion  on  “Current  Problems  in  Soot  Formation  During  Combustion, 
Especially  the  Mechanism  of  Soot  Formation”  was  held  in  Gottingen,  West  Germany  on 
March  29-30,  1989.  The  meeting  was  organized  by  Professor  H.  Gg.  Wagner  and  was 
attended  by  about  twenty  engineers/scientists  currently  involved  with  understanding  soot 
formation  phenomena.  Financial  support  was  provided  by  the  Commission  for  Condensa¬ 
tion  Phenomena  of  the  Academy  of  Sciences  in  Gottingen,  West  Germany. 

2.  Eastern  Section  of  the  Combustion  Institute,  Clearwater  Beach,  Florida,  Dec.  5-7,  1988. 
M.  Colket  presented  an  invited  talk  entitled  “The  Role  of  Oxidative  Pyrolysis  in  Preparticle 
Chemistry” .  See  Appendix  A. 

3.  Brookhaven  National  Laboratory,  Upton,  New  York,  November  9,  1988.  M.  Colket  pre¬ 
sented  an  invited  seminar  entitled  “Kinetic  Mechanisms  for  the  Pyrolysis  of  Unsaturated 
Hydrocarbons”.  Financial  support  provided  by  BNL. 

4.  Department  of  Energy  -  Office  of  Basic  Energy  Sciences  Combustion  Research  Meeting 
held  at  Lake  Geneva,  Michigan,  June  1-3,  1988.  M.  Colket  was  invited  by  W.  H.  Kirchhoff 
to  be  an  observer  and  participant  at  this  D.O.E.  contractor’s  meeting.  The  meeting  was 
attended  under  corporate  sponsorship. 
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The  Role  of  Oxidative  Pyrolysis  in  Preparticle  Chemistry 
M.  B.  Colket  III 

United  Technologies  Research  Center 

Presentation  to  the  Eastern  Section  of  the  Combustion  Institute 
Clearwater  Beach,  Fla.,  Dec.  5-7,  1988 

The  specific  steps  leading  to  soot  formation  have  been  described  recently  by  Classman*.  In  his 
review  article.  Classman  unites  the  individual  efforts  and  results  of  a  great  many  dedicated  re¬ 
searchers  and  presents  a  cohesive  story  of  the  processes  leading  to  soot  formation.  His  analysis 
focuses  on  flame  structure  and  temperature  as  two  of  the  most  important  parameters  in  controlling 
soot  formation:  increasing  temperature  in  a  diffusion  flame  enhances  soot  formation  since  fuel  py¬ 
rolysis  rates  are  increased;  yet,  in  premixed  flames,  increasing  temperature  decreases  soot  formation 
since  competitive  oxidation  processes  are  enhanced.  In  the  present  article,  oxidation  processes  are 
considered.  First,  global  effects  due  to  the  presence  of  oxygen  are  reviewed.  Secondly,  the  oxidation 
of  radical  intermediates  critical  to  soot  formation,  i.e.,  vinyl  and  phenyl,  are  examined,  and  then 
the  phenomena  of  soot-oxidation  are  briefly  discussed. 

Oxygen  is  often  found  in  the  fuel  side  of  diffusion  flames.  For  example,  probe  sampling^’®  has 
demonstrated  the  existence  of  oxygen  at  concentrations  near  1%  in  the  interior  fuel  jet  of  a  co- 
annular  flame.  The  oxygen  may  be  from  finite  rate  chemistry  (diffusion  of  Oj  through  the  flame 
front),  leakage  at  the  burner  lip  or  partial  mixing  (turbulent,  diffusion  flame).  The  effect  of  oxy¬ 
gen  has  been  shown  to  vary  depending  on  fuel-type^:  for  ethene,  oxygen  addition  enhances  soot 
formation,  but  for  alkanes,  soot  formation  remains  the  same  or  actually  decreases.  This  efl'ect  has 
been  explained^  by  the  fact  that,  at  1400  K,  oxygen  addition  noticeably  enhances  the  pyrolysis 
rate  and  radical  production  for  ethene  but  barely  changes  these  rates  in  the  case  of  the  pyrolysis 
of  an  alkane.  Thus,  oxygen  addition  to  ethene  accelerates  the  chemistry  but  has  little  (or  only 
a  dilution)  effect  for  alkanes,  which  have  substantially  weaker  C-C  and  C-H  bonds.  This  result 
is  qualitatively  consistent  with  shock  tube  results  on  soot  scattering®  and  ring  formation®  which 
showed  that  both  soot  and  aromatic  production  are  enhanced  at  lower  temperatures  when  small 
amounts  of  oxygen  are  added  to  the  pyrolysis  of  unsaturated  hydrocarbons.  Numerical  modeling^ 
confirms  the  acceleration  in  chemistry  due  to  the  presence  of  oxygen. 

Quantitative  prediction  of  soot  production,  should  it  become  a  reality,  will  require  knowledge  of 
competition  between  processes  associated  with  ring  formation  and  fragmentation.  As  an  example 
of  the  uncertainty  associated  with  these  processes,  consider  reactions  involving  the  vinyl  radical 
which  has  been  suggested  to  be  a  primary  radical  leading  to  the  formation  of  aromatic  rings. 


CjHs  +  CjHj -(- H  (1) 

CjH3-(-C3H4  — C4H7  -  C,Hr.  +  H  (2) 

C:H3(-(-M)  C7H7 -t- H(*.M)  (3) 

CjHs +  0,  -  HCO+ H7CO  (4) 

C7H3  +  Oj -- HO7  +  C,H7  (5) 


The  relative  rates  of  these  reactions  will  define  the  fate  of  the  vjnyl  radical,  i.e.,  does  it  undergo 
growth,  does  it  decompose  to  acetylene,  or  is  it  oxidized.  Unfortunately,  there  is  a  relatively  large 
amount  of  uncertainty  in  these  rate  constants  and  the  accuracy  of  models  may  therefore  be  in 
question. 


A-2 


The  rate  constant  for  Reaction  1  has  recently  been  measured*  to  be  10*^  *  exp(-5.0  kcal/moIe/RT) 
cm*/mole/sec  and  recent  predictions*  and  most  literature  values  are  reasonably  close  to  this  deter¬ 
mination  near  1200  K.  Two  values  for  kj®  differ  by  an  order  of  magnitude  although  the  direct 
measurement  by  Fahr  and  Stein®  is  probably  the  more  reliable  value. 

Values  for  ks  vary  over  three  orders  of  magnitude.  Very  few  experimental  rate  constants  have 
been  determined  and  these  experiments  have  been  subject  to  alternative  interpretation  or  large 
uncertainty.  A  reasonable  estimate  for  ka  should  be  derivable  from  a  measurement  of  the  reverse 
process  and  detailed  balancing.  Unfortunately,  the  heat  of  formation  of  vinyl  radicals  has  been  under 
much  discussion  recently  and  values  range  from  61  to  72  kcal/mole  (see  Ref.  9).  This  uncertainty 
leads  to  a  range  in  the  calculation  of  ka  of  more  than  two  orders  of  magnitude  at  1200  K. 

Information  regarding  Reaction  4  has  been  obtained  only  at  low  temperatures*^  and  interpreta¬ 
tion  of  these  results  has  recently  been  questioned®.  There  is  no  experimental  confirmation  of  this 
reaction  or  a  rate  determination  at  combustion  temperatures.  Inclusion  of  this  reaction  and  its  low 
temperature  rate  constant**  in  detailed  chemical  models  indicates  that  this  reaction  and  its  rate 
constant  are  critical  in  the  modeling  of  hydrocarbon  oxidation.  This  and  similar  reactions  could 
also  play  a  very  important  role  in  arresting  growth  to  higher  molecular  weight  species.  Reaction  5 
is  perhaps  of  lesser  importance,  and  its  rate  constant  is  also  poorly  known.  Published  values  (see 
Ref.  12)  range  two  orders  of  magnitude.  Evidence  for  its  existence  is  minimal  and  rate  constants 
have  been  inferred  from  detailed  modeling. 

Considering  these  large  uncertainties  and  the  importance  of  these  vinyl  radical  reactions  to  the 
rich  combustion  of  hydrocarabons  and  pre-particle  kinetics,  a  series  of  experiments  were  performed 
in  a  single-pulse  shock  tube  (SPST)  coupled  with  gas  chromatographic  analysis**.  Experimental 
results  for  a  mixture  of  3.5%  ethene,  0.52%  oxygen  and  the  balance  argon  are  shown  in  Figs,  la-ld 
for  total  pressures  near  10  atm  and  average  dwell  times  of  550  microseconds.  Profiles  were  also 
obtained  for  cis-  and  trans-2-butene,  as  well  as  Cs-  and  Cs-  species.  Although  not  shown  in  Fig.  1, 
they  were  included  in  the  mass  balances  shown  in  Figs,  la  and  lb.  Similar  experiments  were  also 
performed  with  no  added  oxygen  (but  there  was  an  oxygen  impurity  of  about  1000  ppm)  and  for 
2.0%  oxygen.  The  pyrolysis  products  were  virtually  the  same  as  those  obtained  by  Skinner  and 
Sokoloski*^  although  decomposition  rates  were  faster.  The  principle  effect  of  oxygen  addition  is 
to  accelerate  the  kinetics  with  little  overall  change  to  the  product  distribution,  although  carbon 
oxides,  ketene,  and  acetaldehyde  were  also  detected.  An  increased  rate  of  aromatic  production  was 
observed  at  lower  temperatures  when  oxygen  was  added  to  the  ethene,  which  is  consistent  with 
calculations^.  The  production  of  aromatics  has  in  turn  been  correlated  with  soot  production**. 

An  analysis  of  the  data  indicated  that  information  on  the  vinyl  radical  could  be  obtained  from 
the  yields  of  1,3-butadiene  (via  Reaction  2).  If  one  assumes  that  (1)  butadiene  decomposition  is 
negligible  at  low  extents  of  reaction  (<1200  K),  (2)  all  oxidized  products  arise  principally  from 
Reaction  (4),  and  (3)  these  products  are  rapidly  converted  to  carbon  monoxide,  then  the  relative 
rates  of  Reactions  2  and  4  can  be  written 


R,  ^  k,  (0,|  _  1  d[CO)/dt 
Rj  kj(C,H<l  2d(C,Hr.|/dt 

where  the  factor  of  one-half  arises  from  the  production  of  both  formyl  and  formaldehyde  (or  two 
CO’s)  in  Reaction  4.  Integrating  this  expression  over  the  length  of  the  experiment  and  rearranging 
gives; 

k,  ^  I  jeoir 

k:  2(C,H6if  (Oil.  ^  ^ 

where  the  subscripts  i  and  f  refer  to  the  initial  and  final  concentrations,  respectively.  If  assumpn 
tions  (1)  or  (2)  are  invalid,  then  Eq.  (B)  represents  an  upper  limit  ratio.  Alternatively,  if  the 
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formaldehyde  produced  in  Reaction  4  does  not  rapidly  decompose,  Eq.  (B)  could  be  a  factor  of 
two  too  high.  A  preliminary  analysis  using^detailed  chemical  kinetics  indicates  that  in  the  temper¬ 
ature  range  1100-1200  K,  about  half  of  the  carbon  monoxide  produced  could  be  formed  (directly 
or  indirectly)  from  other  reactions  such  as  O  or  OH  attack  on  ethene.  In  addition,  not  all  of  the 
formaldehyde  decomposes  to  carbon  monoxide  during  the  SPST  experiments.  Consequently,  we 
make  the  preliminary  assessment  that  rate  constants  determined  from  Elq.  (B)  represent  an  upper 
limit  and  values  reduced  by  a  factor  of  one-third  are  reasonable  estimates  for  this  rate  constant.  The 
upper  limit  and  the  reduced  values  are  compared  to  the  value  obtained  by  Slagle  and  co-workers" 
in  Fig.  2.  In  addition,  a  rate  constant  was  determined  from  flow  reactor  data^®  using  Eq.  (A) 
(for  an  experiment  of  3.1%  ethene,  0.5%  oxygen  at  one  atmosphere).  The  value  (also  reduced  by 
one-third)  is  shown  in  Fig.  2  and  agrees  well  with  the  SPST  data.  Overall,  these  new  values  are 
consistent  with  an  extrapolation  of  the  low  temperature  results,  although  a  stronger  decay  with 
increasing  temperature  is  indicated.  This  negative  temperature  dependence  is  not  unexpected  for 
addition/rearrangement  reactions.  As  a  test  for  this  technique,  the  same  procedure  was  used  in  the 
case  when  2.2%  ethene  and  1.1%  acetylene  were  co-pyrolyzed.  Thus,  a  comparison  of  rate  constants 
for  reactions  1  and  2  could  be  obtained  for  comparison  to  other  data.  By  analogy  to  the  previous 
derivation, 


kj  (C^Hejf  [CjHili 

This  ratio  wais  found  to  be  1.25±0.3  over  the  temperature  range  1240  to  1360  K,  and  compares  well 
to  the  ratio  of  1.42  obtained  by  Fahr  and  Stein  at  slightly  lower  temperatures. 

Addition  of  a  little  oxygen  to  ethene  pyrolysis  (as  for  many  other  hydrocarbons)  results  principally 
in  the  acceleration  of  the  reaction:  some  carbon  monoxide,  formaldehyde,  ketene  and  other  oxidative 
products  are  formed,  but  the  identity  of  the  hydrocarbon  species  and  their  relative  concentrations 
are  similar  in  the  case  of  pure  pyrolysis  and  oxidative  pyrolysis.  For  SPST  experiments  at  UTRC 
this  situation  does  not  hold  in  the  case  of  benzene.  Based  on  the  following  analysis,  the  difference 
is  due  totally  to  the  different  fate  of  the  phenyl  radical  in  the  pyrolysis  versus  that  in  the  oxidative 
pyrolysis.  In  the  pyrolysis, 


Comparison  of  Rote  Consfonts  for 
C2H3  +  02-^HC0+  H2CO 


10A(VK) 


CsHs  -  C^Ha  +  CjHj  (overall)  (6) 

-l-H  (7) 

♦-*  CjH  -f  CjHj  (8) 

dominate,  leading  to  the  production  of  CjHj 
and  as  the  principal  species"  *®.  In 

the  oxidative  pyrolysis,  however,  a  sequence 
similar  to 


CijHs  -4-  Oj  — '  C'r.HsO  -f  O  (9) 

C«HsO  -  CjHs -4- CO  (10) 

CsHs  +  Oi -CsHsO  +  O  (11) 

CsHsO  -  C4H5 +CO  (12) 

C%Hs-C,H<  +  H  (13) 

C^Hs  -  CjHj -I-C2H3  (14) 

C2H3 -- CjH: -t- H  (15) 
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overrides  the  pyrolytic  ring  fracture  and  phenol,  cyclopentadiene,  vinylacetylene,  and  acetylene 
are  principle  products.  These  sequences  ar^  qualitatively  consistent  with  observed  SPST  product 
distributions.  This  mechanism  for  phenyl  oxidation  is  also  consistent  with  and  derived  principally 
from  the  work  of  Venkat,  et  aP®.  The  explanation  for  this  major  difference  between  the  pyrolysis 
and  the  oxidative  pyrolysis  is  that,  at  1400  K,  the  lifetime  of  phenyl  based  on  Reaction  6  is  3500 
microseconds,  but  for  decomposition  via  Reaction  9,  the  lifetime  is  just  40  microseconds  with  only 
17  torr  Oj.  Biphenyl  was  also  observed  in  the  oxidative  pyrolysis  runs  since,  for  these  series  of 
experiments,  the  sampling  system  had  been  modified  to  collect  polyaromatic  species.  Since  biphenyl 
is  formed  principally  via 


CeHe  +  CsHs-CiiHio  +  H  (16) 

The  relative  rates  of  Reactions  9  and  16  can  be  written 

Rg  _  [Qi] 

Ri6  ki6  (CeHe] 

Assuming  negligible  back  reaction,  the  integrated  value  of  Rir  is  simply  the  yield  of  biphenyl, 
while  that  of  Rg  can  be  estimated  by  either  one-half  (1)  the  yields  of  phenol,  cyclopentadiene,  and 
vinylacetylene  plus  one-half  the  yield  of  acetylene  or  one-half  (2)  the  yield  of  phenol  plus  that  of 
carbon  monoxide  minus  the  CO  produced  from  Reaction  12  (the  yield  of  vinylacetylene  plus  one-half 
of  that  from  acetylene).  The  correction  of  one-half  is  to  account  for  the  fact  that  0-atoms  produced 
in  Reaction  6  add  to  benzene  to  form  phenol.  Method  (2)  produced  values  for  the  integrated  rate  of 
Reaction  9  about  40%  greater  than  those  of  method  (1).  Since  additional  phenol  and  hence  other 
products  can  also  be  formed  from  O  and  OH  via 

H  +  Oj  — OH  +  0 
O  -I-  CsHs  CsHsOH 
OH -1- CsHs  -  CsHsOH  +  H, 

the  integrated  rate  of  Reaction  9  calculated  using  these  methods  is  an  upper  limit.  Using  values  av¬ 
eraged  from  methods  (1)  and  (2),  the  ratio  kg/kie  could  be  calculated.  Using  a  recent  measurement® 
of  kjc,  values  of  kg  were  determined  (from  a  run  of  1.1%  benzene  and  0.2%  oxygen)  and  are  plot¬ 
ted  in  Fig.  3.  Also  shown  on  this  figure  are  literature  values  based  on  detailed  modeling*®  and 
thermochemical  estimates*'.  The  SPST  upper  limits  are  consistent  with  the  previous  results. 

The  drastic  difference  between  products 
of  the  pyrolysis  and  the  oxidative  pyrolysis 
of  benzene  (and  quite  likely  other  aromatic 
hydrocarbons)  could  have  important  im¬ 
plications  for  the  mechanisms  of  soot  for¬ 
mation.  Cyclopentadienyl,  whose  stabil¬ 
ity  is  resonantly  enhanced,  is  an  impor¬ 
tant  and  relatively  long-lived  intermediate 
during  the  oxidative  pyrolysis  of  benzene, 
but  is  unimportant  in  the  pyrolysis.  It  is 
quite  reasonable  to  expect  that  collisions 
involving  cyclopentadienyl  and  aromatics 
(or  fused  six-membered  rings)  may  lead  to 
growth  of  aromatic  rings.  In  fact,  an  un¬ 
usual  result  from  the  present  experiments 
is  the  production  of  1-methylnaphthalene 
at  the  early  (low  temperature)  stages  of  ox- 


CeHs+Oz 
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idative  pyrolysis  of  benzene.  Concentrations  of  this  product  are  nearly  a  factor  of  three  higher  than 
that  of  naphthalene  and  100  times  higher  than  2-methylnaphthalene.  A  possible  explanation  is  the 
reaction 


c-CsHs  +  CeHfl  — *  H  +  I-C10H7CH3  (1-methylnaphthalene). 

Although  the  rate  constant  for  this  reaction  might  be  low  due  to  a  complex  rearrangement,  the  net 
rate  could  be  high  due  to  high  concentrations  of  cyclopentadienyl.  In  fact,  the  naphthalene  that  is 
observed  could  be  a  result  of  this  reaction  followed  by  H-atom  substitution  of  the  methyl  group  rather 
than  acetylene  addition  to  phenyl  and  then  to  the  benzyne-acetylene  radical**.  Another  interesting 
aspect  of  oxidation  phenomena  relates  to  the  question  of  six-membered  vs.  five-membered  rings 
in  pre-particle  kinetics.  Thermodynamic  and  other  arguments**  have  been  used  to  suggest  that 
principally  six-membered  rings  exist  in  pre-particle  molecules.  These  arguments  contrast  with 
interpretation  of  experimental  data*^'*®  supporting  the  existence  of  buckminsterfullerene  and  similar 
species  which  contain  a  mixture  of  six-  and  five-membered  rings.  Based  on  (l)  the  interpretation 
of  benzene  (and  phenyl)  oxidation,  which  includes  formation  of  either  phenoxy  and  phenol  due  to 
attack  by  O2,  O,  or  OH  followed  by  loss  of  a  single  carbon  (to  carbon  monoxide)  and  the  formation  of 
a  five-membered  ring  and  (2)  the  interpretation  by  Harris*®  that  oxidation  of  nascent  soot  particles 
is  a  detectable  process  in  fuel-rich  premixed  flames,  we  reach  the  following  conclusion;  during  the 
growth  of  PAH’s,  partial  oxidation  of  these  species  will  occur  leading  to  the  production  of  carbon 
monoxide  and  five-membered  rings.  The  resultant  compounds  will  contain  mixtures  of  five-  and 
six-membered  rings,  some  of  which  will  fragment,  and  others  may  continue  to  grow  via  addition 
processes. 

Soot  Oxidation 

The  focus  of  this  paper  thus  far  has  been  the  oxidation  of  small  hydrocarbon  radicals.  Oxidation 
of  larger  species  and  soot  also  play  an  important  role  in  sooting  flames.  Harris*®  recently  interpreted 
experimental  data  to  conclude  that  oxidation  competes  with  surface  growth  and  coagulation  during 
particle  inception,  and  it  is  generally  acknowledged  that  oxidation  can  retard  growth  throughout 
the  life  of  a  particle.  The  parameter  controlling  complete  burn-out  of  soot  has  been  determined 
to  be  temperature.  Kent  and  Wagner**  have  shown  that  once  the  flame  temperature  decreases 
below  1400  K,  soot  oxidation  ceases,  and  a  soot  trail  emerges  from  the  tip  of  the  flame.  A  physical 
interpretation  of  this  process  can  be  found  by  considering  a  heat  balance  on  an  isolated  soot  particle. 
Assuming  the  dominant  factors  to  be  radiation  loss  (Qr>d)  and  exothermic  heat  rele2ise  due  to 
oxidation  (Qox)i  the  following  simplified  analysis  can  be  used  to  explain  quantitatively  the  cessation 
of  soot  burn-out. 

Radiative  loss  can  be  calculated  based  on  an  expression  from  Hall*®: 

=  2.97  X  10“‘°  watts/cmVK^fvT® 
for  radiation  from  unit  gas  volume,  or 

Q„d  =  2.97  X  I0~‘°  walts;cm*/K"T^ 

=  7.1  X  10”’^  kcal/sec/cm*/K*T* 

for  radiation  from  unit  particle  volume,  fv  represents  the  particle  volume  fraction  and  the  fifth 
power  dependency  on  temperature  arises  from  a  1/A  dependency  on  the  particle  emissiviiy,  A 
comparison  of  oxidation  rates  indicates  (see  e.g..  Ref.  27)  that  the  oxidation  of  soot  particles  is 
nearly  independent  of  the  concentration  of  molecular  oxygen  below  2000  K.  A  fit  to  measured  rates 
has  been  reported**  and  the  rate  of  soot  oxidation  is 
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A 

1.  — 

fit 

with  kox  =  91  X  g/cm^/sec 

and  Sn,  is  the  soot  mass  and  A  is  the  surface  area  of  the  particle.  The  rate  expression  is  valid  only 
for  temperatures  below  2000  K.  Assuming  the  reaction  sequence  to  be  equivalent  to 

C(s)  +  1/2  Oj  CO. 

the  reaction  is  21  kcal/mole  (or  1.75  kcal/gm)  exothermic.  The  assumption  that  carbon  monoxide 
is  the  principle  product  is  consistent  with  recent  results^®  on  carbon  oxidation  using  a  diode  laser 
for  product  detection.  The  heat  release  due  to  oxidation  is  therefore 


=  (1.75  kcal/gm)  x 


40  rm  partide 


Fig.  4 


=  160  X  A  X  kcal/cmVsec. 

CompOTSOn  of  Rotes  of  Heat  Gocytoss  .  Assuminj  »  particle  diameter  of  40  nm, 
For  on  OxidMng/Racliating  Soot  Partide  ,  ■“  ” 

f _  temperature  for  comparison.  As  long  as 

40  nm  partide  ^  particle  temperature  exceeds  approxi- 

mately  1400  K,  the  oxidation  process  gen- 
crates  sufficient  heat  to  overcome  that  loss 
_io-  by  radiation.  Temperatures  can  rise  in  this 

'  . . region  (see  Ref.  30)  but  usually  decay  due 

-ti-  Qrod  _  to  other  (conductive  and  convective)  losses. 

Furthermore,  the  temperature  may  de- 
-t2-  creaise  since  heat  release  due  to  oxidation 

may  be  inhibited  by  the  inability  of  oxygen 
“'2-  to  diffuse  to  the  center  of  the  flame.  Once 

^1800K^  ^  ^  U(^  ^  COOK  TKWK  particle  temperature  reaches  about 

I  r  1400  K,  however,  the  oxidation  quenches, 

nnn/T  flAO  “  radiation,  conductive  and  convective 

'  ^  processes  rapidly  cool  the  particle. 

Characteristic  times  for  oxidation  can  also  be  estimated,  assuming  the  particle  is  a  perfect  sphere. 

^  =  -kA 
dt 

substituting  for  the  surface  area  and  integrating 


0.7  0.6 

xxx/rW 


_  /■*"'  4>rkdt 


t„x  =  ^ 


=  4.0  X  sec 

=  62  msec  at  1400  K 
=  24  msec  at  1500  K 
=  10  msec  at  1600  K 
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for  a  particle  40  nm  in  diameter  and  with  a  density  of  1.8  gm/cm^.  Considering  dwell  times  in  the 
luminous  zone  of  a  co-annular  diffusion  flame,  as  well  as  measured  temperatures^*,  one  can  easily 
argue  that  once  the  particle  temperature  drips  below  1500  K,  there  is  insufficient  time  for  complete 
oxidation.  Conductive,  convective,  and  radiation  losses  slightly  dominate  over  the  chemical  heat 
release  and  the  particle  cools  slowly  until  1400  K,  at  which  point  heat  release  from  oxidation  is  too 
slow,  the  particle  cools  very  rapidly,  and  a  soot  trail  is  emitted. 
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The  Rich  Oxidation  of  Ethylene  in  a  Single-Pulse  Shock  Tube 


M.  B.  Colket  III 

United  Technologies  Research  Center 

Second  International  Conference  on  Chemical  Kinetics 
Gaithersburg,  MD.,  July  23-27,  1989 

The  fate  of  vinyl  radicals  during  pyrolysis  and  combustion  of  hydrocarbons  has  been  a  subject  of 
much  discussion.  Possible  reactions  include 


C2H3  +  C2H2  ^(C^Hs)'  ^  C4H4 -I- H  (1) 

C2H3  +  C2H4--(C4H7)*  --C4H6  +  H  (2) 

C2H3(+M)  ^  C2H2  +  H(+M)  (3) 

C2H3  +  O2  ->  HCO  -h  H2CO  (4) 

C2H3  +  O2 -- HO2  +  C2H2  (5) 


The  relative  rates  of  these  reactions  will  define  the  fate  of  vinyl;  yet  except  for  the  first  two, 
the  rate  constants  for  these  processes  are  poorly  known.  To  shed  some  light  on  these  steps,  a 
series  of  experiments  using  a  single-pulse  shock  tube  have  been  completed  on  the  rich  oxidation  of 
ethylene.  Post-shock  temperatures  ranged  from  about  1000  to  1800K  and  total  pressures  were  9  to 
12  atmospheres.  Ethylene  and  oxygen  were  mixed  in  a  bath  of  argon.  Concentrations  of  ethylene 
were  0.3  to  3.5%  and  overall  stoichiometries  varied  from  5.25  to  about  100.  Final  products  were 
collected  and  analyzed  using  gas  chromatography.  Measured  products  include  Cj-Cio  hydrocarbons, 
carbon  oxides,  and  hydrocarbon  oxygenates  such  as  acetaldehyde  and  ketene.  Major  hydrocarbon 
products  include  acetylene,  methane  and  ethane.  Detailed  chemical  kinetic  modeling  indicates  that 
the  latter  is  formed  principally  by  the  sequence 

C2H4  +  H(  +  M)^C2H5(  +  M)  (6) 

C2H5  +  C2H4  ^C2H6  +  C2H3  (7) 

Direct  evidence  for  Reaction  4  (oxidation  of  vinyl)  was  not  obtained  but  an  upper  limit  for  its  rate 
constant  was  found  to  be  3  x  10‘*  cm^/mole/sec  at  HOOK  which  is  in  good  agreement  with  the 
determination  by  Slagle  and  coworkers  (1984).  ks  was  estimated  to  be  5  x  10*°  cm®/mole/sec  at 
HOOK.  Large  amounts  of  acetylene  are  observed  and  are  indicative  of  either  an  unusually  high  rate 
constant  for  k3  or  very  significant  contributions  due  to  reactions,  such  as: 


C2H4(  +  M)->C2H2+H2(  +  M)  (8) 

C2H3  -l-  C2H4  ‘-+(04117)  <-*  C2H5  -f  C2H2  (9) 

By  copyrolyzing  ethylene  and  acetylene,  the  ratio  ki/k2  was  determined  to  be  1.25±0.3  from  1240  to 
1360K  and  agrees  well  with  the  ratio  of  1.42  obtained  from  the  Fahr  and  Stein  (1989)  measurements. 

ACKNOWLEDGEMENTS 

This  work  has  been  supported  in  part  by  the  Air  Force  Office  of  Scientific  Research  under  Contract 
Nos.  F49620-85-C-0012  and  F49620-88-C-0051 . 


B-2 


Appendix  C 

The  Pyrolysis  of  Acetylene 
Initiated  by  Acetone 


C-1 


COMBUSTION  AND  FLAME  75:  343-366  (1989) 


The  Pyrolysis  of  Acetylene  Initiated  by  Acetone 

M.  B.  COLKET  III  and  D.  J.  SEERY 

United  Technologies  Research  Center,  East  Hartford,  CT  06108 

and 

H.  B.  PALMER 


Pennsylvania  State  University,  University  Park,  PA  16802 


A  detailed,  radical  chain  mechanism  is  used  to  model  the  pyrolysis  of  acetylene  near  1000  K.  The  initiation  process. 
CjHj  +  C;H;  —  C4H1  +  H.  appears  to  be  inconsistent  with  thermochemistry.  Since  experimental  evidence  indicates 
the  presence  of  a  chain  mechanism,  alternative  sources  of  initiation  are  considered.  Acetone,  a  common  impurity  in 
■■purified"  acetylene,  was  found  to  dominate  radical  initiation  during  the  pyrolysis  of  acetylene  near  1000  K  despite 
concentration  levels  only  0. 1  %  that  of  acetylene.  Modeling  results  compare  favorably  with  the  experimental  results  of 
Munson  and  Anderson  for  acetylene  decay  and  the  formation  of  products  vinylacetylene,  benzene,  and  ethylene.  Rate 
constants  for  the  following  reactions  were  adjusted  to  be 

CjHi  +  C:Hj  ~  k,  =  2.Qx  10”  cm’  mol''  sec"' 

n-C4H5  C4H4  +  H  *j2=  l.6x  sec' 

n-C4H5  +  C;H2-’CtH*+  H  A'jj  =  4.5  x  10'’ xe'"®'  '^■'''rr)  cm’  mo!"'  sec' ' 

n-C4H,  +  C4H4  -C,H,  +  H  =  1 .2  X  1 0  ”  X  e  -  cm  ’  mol ' '  sec  ' ' 

to  optimize  the  fit  to  the  experimental  data.  A  sensitivity  analysis  shows  that  the  computed  results  were  most  sensitive 
to  rate  constants  for  these  and  a  few  other  reactions.  By  addition  of  high-temperature  reactions  (such  as  those 
involving  ethynyl  radicals),  the  mechanism  was  found  to  be  consistent  with  experimental  results  using  single-pulse 
shock  tubes  and  the  very-high-temperature  studies  of  Frank  and  Just  and  Wu,  Singh,  and  Kem.  Calculations  are 
quantitatively  consistent  with  the  experimental  finding  that  the  effect  of  acetone  is  negligible  at  high  temperatures 
(>2000  K)  and  low  concentrations  of  acetylene  (<5  x  10''®  mol  cm*’).  The  model  indicates,  however,  that  at 
either  higher  concentrations  or  lower  temperatures,  impurities  can  play  a  significant  role. 


INTRODUCTION 

Speculation  on  the  mechanism  of  acetylene  pyroly¬ 
sis  has  continued  for  many  years.  Experiments  at 
low  temperature  (6(X)-1200  K)  during  the  1950s 
and  1960s  resulted  in  a  variety  of  suggestions, 
including  radical  chain  processes  and  nonchain 
dimerizations  involving  diradical  or  electronically 
excited  intermediates.  During  the  1970s,  a  radical 
chain  mechanism  was  preferred  [1-3]  and  initia- 
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tion  was  assumed  to  be  caused  by 

C2H2  +  C2H2-C4H,-|-H.  (1) 

The  structure  of  was  not  identified  [1]  but 
the  overall  endothermicity  45.9  kcal/mol  was 
consistent  with  an  early  calculation  [4]  for  the  heat 
of  formation  of  the  butatrienyl  radical  (102  kcal/ 
mol).  FrenkJach  and  coworkers  [5],  in  a  recent 
modeling  study,  tracked  two  of  the  isomers, 
namely  l-buten-3-yn-l-yl  (n-CaHj)  and  1-buten- 
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3-yn-2-yI  (i-C4H3),  and  used  Tanzawa  and  Gar¬ 
diner’s  rate  constant  [1],  but  assumed  that  Reac¬ 
tion  la  was  the  initiation  reaction. 


C2H2  +  C2H2-n-C4H3  +  H, 

(la) 

C2H2  +  C2H2“^i”C4H3  +  H. 

(lb) 

Recent  estimates  [6,  7]  for  the  heat  of  formation 
of  the  product  of  Reaction  la  (n-C4H3)  of  124 
kcal/mol  (see  later  discussion)  indicate  that  this 
reaction  is  about  68  kcal/mol  endothermic.  Wu  et 
al.  [8]  noted  the  descrepancy  between  the  en- 
dothermicity  of  Reaction  la  and  the  observed 
activation  energy,  and  suggested  that  i-C4H3  is  the 
product  (Reaction  lb).  Because  isomerization  of 
the  (C4H4)  intermediate  is  required  prior  to  de¬ 
composition,  it  is  probable  that  the  reaction  will 
have  a  low  A-factor  and/or  an  activation  barrier  in 
excess  of  the  reaction  endothermicity.  Their  pro¬ 
posed  Arrhenius’ parameters  (A  =  2  X  10'^  cm^ 
mol"‘  sec"',  E  ~  44.5  kcal  mol"')  are  inconsis¬ 
tent  with  this  constraint. 

Discounting  the  possibility  of  isomerization  and 
invoking  an  activation  energy  of  68  kcal/mol  for 
Reaction  la,  Duran  et  al.  [9]  were  unable  to 
describe  experimental  data  with  detailed  chemical 
kinetic  modeling.  In  addition,  Duran  and  cowor¬ 
kers  [10]  found  their  data  to  be  consistent  with  the 
assumption  that  acetylene  pyrolyzes  by  isomeriz- 
ing  to  vinylidene  which  then  inserts  into  a  C-H 
bond  of  acetylene  to  form  vinylacetylene.  The 
mechanism 

C2H2  +  M  «  HjCCi-t-M,  (Ic) 

H2CC:-)-C2H2  -  (C4H4)*,  (Id) 

(C4H4)*  +  M  «  C4H4-t-M  (le) 

requires  that  the  reverse  of  the  first  step  be  very 
fast  so  that  vinylidene  is  always  in  equilibrium 
with  acetylene.  Alternatively,  Kiefer  et  al.  [11] 
find  support  for  the  vinylidene  mechanism  from 
detailed  balancing  and  interpretation  of  experi¬ 
mental  data  on  the  reverse  process,  i.e.,  the 
decomposition  of  vinylacetylene. 

We  cannot  disprove  a  contribution  from  this 
nonchain  dimerization;  however,  upon  review  of 
our  own  experimental  data  as  well  as  that  obtained 
by  others  we  reach  the  unavoidable  conclusion  that 
a  radical  chain  process  occurs.  In  attempts  to 
model  experimental  data,  we  have  identified  the 


probable  role  of  acetone,  a  common  impurity  in 
acetylene,  in  initiating  a  radical  chain  during 
laboratory  studies  of  the  pyrolysis  of  acetylene. 

In  this  paper,  evidence  supporting  a  radical 
chain  is  presented,  issues  related  to  acetone- 
initiation  are  discussed,  and  modeling  results  are 
compared  to  existing  experimental  data  near 
1000-1500  K.  Finally,  the  proposed  model  is 
extended  to  examine  the  possible  role  of  acetone 
under  higher  temperature  conditions  (<2000  K). 

EVIDENCE  FOR  CHAIN  MECHANISM 

Palmer  and  Cullis’s  review  article  [12]  cites 
results  such  as  inhibition  by  radical  scavengers  and 
wall  effects  as  evidence  that  the  thermal  decompo¬ 
sition  of  acetylene  is  a  radical  process.  Because  it 
is  not  known  whether  radical  scavengers  and  walls 
might  effectively  remove  singlet  vinylidene  as 
well  as  monoradicals,  other  evidence  is  needed 
regarding  the  possible  role  of  vinylidene.  The 
arguments  [10,  13]  for  an  extremely  short  lifetime 
(with  respect  to  isomerization  to  HCCH)  of 
vinylidene  do  strongly  suggest  that  a  dominant  role 
for  vinylidene  is  not  compatible  with  the  experi¬ 
mental  observation  of  induction  periods  (see  ref. 
12)  in  the  pyrolysis;  Reaction  (Ic,  -  Ic)  will 
essentially  equilibrate  in  a  time  equal  to  (or  less 
than  [13])  the  time  required  for  vinylidene  to 
undergo  one  collision. 

Therefore,  the  observation  of  induction  periods 
supports  the  participation  of  a  radical  chain  in 
which  vinylidene  does  not  play  a  significant  role. 
[Although  a  hypothetical  chain  involving  vinyli¬ 
dene  can  be  constructed  by  adding  a  fourth 
reaction, 

(C4H4)*  +  C2H2  H2CC:  -I- C4H4  (If) 

to  the  above  three,  the  “chain  length"  in  this 
sequence  will  be  insignificant  because  of  the 
extremely  rapid  isomerization  of  H2CC:  to  acety¬ 
lene,  i.e.,  R(10  cannot  compete  effectively  with 
R(-lc).] 

Evidence  for  a  radical  chain  also  arises  from  the 
work  of  Callear  and  Smith  [14],  who  investigated 
the  reaction  between  acetylene  and  atomic  hydro¬ 
gen  in  a  quartz  vessel  at  293  K.  Atomic  hydrogen 
was  generated  by  decomposing  H2  using  a  low- 
pressure  mercury  discharge  lamp.  Major  products 
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included  ethene,  1 ,3-buiadiene,  benzene,  and 
trans-l,3,5-hexatriene.  To  explain  these  products, 
Callear  and  Smith  suggested  a  reaction  mechanism 
similar  to  that  used  in  the  present  study  except  that 
in  the  low  temperature,  hydrogen-rich  investiga¬ 
tion,  C2H3,  n-C4Hj,  and  I-C6H7  (1,3,5-hexa- 
trienl-yl)  radicals  formed  principally  C2H4,  1,3- 
C4H6,  and  l,3,5-C6Hg  rather  than  losing  H-atoms 
to  form  C2H2,  C4H4,  and  benzene.  Some  benzene 
was  still  formed  at  low  temperatures,  presumably 
because  of  the  high  stability  of  the  aromatic.  A 
main  conclusion  from  this  low-temperature  study 
was  that  chain  processes  exist.  If  these  chain 
processes  exist  at  room  temperature,  they  are  not 
likely  to  vanish  at  elevated  temperatures.  The 
principal  concerns  in  extending  these  results  to 
pyrolysis  conditions  are  the  stability  of  the  radical 
adducts  (relative  to  the  reactants)  and  whether  or 
not  sufficient  radicals  can  be  generated  to  initiate 
the  chain  and  account  for  the  measured  rate  of 
acetylene  pyrolysis.  These  aspects  can  only  be 
addressed  properly  through  detailed  chemical  ki¬ 
netic  modeling. 

At  elevated  temperatures,  there  is  a  multiplicity 
of  products  from  the  thermal  decomposition  of 
acetylene.  Near  1000  K,  using  20%  C2H2  in 
helium,  Munson  and  Anderson  [15]  found  vinyla- 
cetylene  and  also  benzene,  ethene,  methane,  and 
hydrogen.  Working  at  even  higher  temperatures, 
Ogura  [3]  detected  the  same  products  during 
pyrolysis  in  a  single  pulse  shock  tube  and  also 


diacetylene  and  traces  of  1,3-butadiene,  and  C3 
and  C5  species.  Under  similar  conditions  but  with 
a  heated  sampling  system,  Colket  [16]  detected  the 
same  species  and  benzene,  phenylacetylene,  and 
traces  of  naphthalene  as  well.  Although  the  vinyli- 
dene  mechanism  apparently  explains  the  overall 
decomposition  rates  of  the  reactant,  it  does  not 
describe  the  variety  of  products  that  have  been 
observed  experimentally.  Thus,  contributions 
from  a  chain  sequence  are  expected. 

Further  indications  of  a  chain  mechanism  are 
provided  by  additional  single-pulse  shock  tube 
data  in  Fig.  1,  which  compares  vinylacetylene  and 
benzene  production  from  three  different  gases;  (a) 
2040  ppm  biacetyl  in  argon,  (b)  3.7%  acetylene  in 
argon,  and  (c)  a  mixture  of  3.55%  acetylene  and 
1500  ppm  biacetyl  in  argon.  It  also  reveals  that  the 
production  of  C4H4  and  is  substantially 

larger  from  the  mixture  than  would  be  expected 
from  purely  additive  considerations.  Near  1 100  K, 
the  mixture  produces  an  order  of  magnitude  more 
C4H4  and  C6H6  than  is  produced  from  acetylene  or 
biacetyl.  Because  there  is  no  reason  to  expect  that 
biacetyl  enhances  the  rate  of  the  vinylidene  mecha¬ 
nism,  this  enhancement  must  be  a  result  of  a  chain 
mechanism  initiated  by  biacetyl  decomposition 
first  into  acetyl  and  then  into  methyl  radicals. 
Mechanisms  by  which  methyl  radicals  are  con¬ 
verted  to  H-atoms  (required  for  the  acetylene 
chain)  are  discussed  elsewhere  in  this  article. 
Methyl  radical  addition  to  acetylene  and  to  subse- 
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Fig.  I .  Comparison  of  production  of  vinylacetylene  and  benzene  in  a  single-pulse  shock  ntbe. 

- ,  3.7%  acetylene;  — .  2040  ppm  biacetyl;  •  ■  3.55%  acetylene/ 1 500  ppm  biacetyl:  □, 

vinylacetylene;  A,  benzene.  Total  pressure  -  8  atmo,  dwell  times  ~  500-700  jisec. 
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quent  adducts  forms  a  variety  of  products,  e.g., 
C3H4,  C4H6,  C4Hg,  etc.,  but  the  concentrations  of 
these  species  are  small  relative  to  those  of  C4H4 
and  C6H6.  There  appears  to  be  no  facile  method  to 
produce  vinylacetylene  and  benzene  at  the  ob¬ 
served  rates  other  than  from  the  chain  mechanism 
related  to  acetylene  pyrolysis. 

Finally,  Ogura  [3]  shock-heated  equimolar  mix¬ 
tures  of  C2H2/C2D2  diluted  in  argon.  He  found 
more  than  40%  of  the  vinylacetylene  in  the  form  of 
the  C4H3D  and  C4HD3  isotopes  at  temperatures  as 
low  as  1070  K.  It  is  highly  unlikely  that  such 
significant  scrambling  can  be  accounted  for  by  a 
vinylidene  mechanism.  Again  a  radical  chain 
sequence  is  favored. 

In  summary,  (1)  an  induction  period  has  been 
observed  during  low-temperature  acetylene  pyro¬ 
lysis,  (2)  a  radical  chain  mechanism  has  been 
observed  at  room  temperature,  (3)  acetylene  pyro¬ 
lysis  produces  a  wide  variety  of  products,  (4)  a 
large  enhancement  in  the  production  of  vinylacety¬ 
lene  and  benzene  is  observed  when  acetylene  is 
copyrolyzed  with  biacetyl,  and  (5)  significant 
scrambling  was  observed  during  pyrolysis  of 
C2H2/C2D2  mixtures.  Because  these  phenomena 
are  either  indicative  of  a  chain  mechanism  or  can 
only  be  explained  by  a  chain  mechanism,  they 
support  the  existence  of  a  chain  mechanism  during 
acetylene  pyrolysis  at  tempieratures  of  900-1400 
K.  They  do  not  exclude  possible  contributions^  a 
vinylidene  or  other  nonchain  mechanism,  nor  do 
they  provide  information  on  the  relative  contribu¬ 
tions  of  the  two  mechanisms  if  two  coexist. 
Nevertheless,  as  recent  modeling  efforts  have 
focused  on  a  nonchain  mechanism,  this  work 
focuses  solely  on  chain  processes  to  examine 
whether  or  not  they  describe  acetylene  decomposi¬ 
tion  and  product  formation  adequately. 

HEATS  OF  FORMATION  OF  C4H3 
ISOMERS 

An  interesting,  although  not  critical,  issue  of  this 
analysis  is  the  heats  of  formation  of  the  isomers  of 
C4H3.  Most  modelers  have  relied  on  the  early 
calculations  of  Cowperthaiie  and  Bauer  [4],  which 
were  performed  for  the  butatrienyl  isomer.  Later 
estimates  for  the  C4H3  isomers  differ  substantially 
from  these  calculations.  Nevertheless,  recent  com¬ 


pilations  [17,  18]  have  retained  the  early  results 
and  differences  between  the  isomers  are  often 
ignored.  For  the  purpose  of  this  study,  we 
hypothesize  two  separate  isomers  of  C4H3:  HCC- 
CCH2  (i-C4H3)  and  HCCCHCH  (n-C4H3),  We 
assume  the  first  of  these  to  be  a  resonant-hybrid 
with  the  two  Kekule  structures:  HC  =  C  =  C  =  CH2 
and  HC  =  C-C  =  CH2.  This  assumption  is  consist¬ 
ent  with  interpretation  of  electron  spin  resonance 
measurements  [19]. 

Recent  estimates  by  Stein  [6]  and  Bittner  [7] 
place  the  heat  of  formation  (A///.298°)  of  n-C4H3 
near  124  kcal/mol.  These  estimates  are  based  on 
the  heat  of  formation  of  C4H4  (about  69  kcal/mol) 
and  the  assumption  that  D(n-C4H3-H)  is  equiva¬ 
lent  by  analogy  to  D(C2H3-H),  i.e.,  about  108 
kcal/mol.  Accepting  the  analogous  reasoning, 
there  are  remaining  uncertanties  in  both  the  heat 
formation  of  vinylacetylene  (C4H4)  and  the  C2H3- 
H  bond  strength.  Recent  proposals  [17,  20]  for  the 
heat  of  formation  of  C4H4  vary  between  68  and  74 
kcal/  mol.  We  chose  to  accept  69  kcal/mol  based 
on  the  agreement  between  group  additivity  calcu¬ 
lations  and  Bittner’s  analysis  [7],  which  compared 
experimental  heats  of  hydrogenation  for  several  3- 
en-l-yne  hydrocarbons.  The  heat  of  formation  of 
vinyl  and  consequently  the  vinyl-H  bond  strength 
has  been  under  much  discussion  in  recent  years 
with  most  suggested  values  ranging  from  62  [17] 
to  72  [21]  kcal/mol.  A  midrange  value  of  67.5 
kcal/mol  leads  to  124(-i- 10,  -7)  kcal/mol  for  the 
heat  of  formation  of  n-C4H3,  where  the  assigned 
uncertainties  reflect  the  range  in  the  published 
values  for  the  heats  of  formation  of  vinyl  radical 
and  vinylacetylene. 

Thermodynamics  for  i-C4H3  have  also  been 
estimated  by  both  Stein  [6]  and  Bittner  [7].  Stein 
argues  that  the  difference  between  n-C4H3  and  i- 
C4H3  should  be  approximately  equal  to  the  reso¬ 
nant  energy  in  HCCCCH2.  An  upper  limit  of  this 
resonant  energy  is  that  in  the  propargyl  radical 
(  —  8  kcal/mol  [58]).  He  then  estimates  a  lower 
limit  heat  of  formation  of  i-C4H3  to  be  116  kcal/ 
mol.  Alternatively,  Bittner  used  the  procedure 
described  by  Cowperthwaite  and  Bauer  [4]  to 
determine  the  delocalized  electron  interactions. 
Using  parameters  updated  from  those  used  previ¬ 
ously  [4],  Bittner  calculates  A///,298°  =  110.6 
kcal/mol.  It  is  important  to  note  that  as  in  the  case 
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of  n-C4H3,  both  estimates  of  i-C4H3  are  depen¬ 
dent  on  the  uncertainties  of  the  heats  of  formation 
of  C4H4  and  C2H3.  Both  the  Stein  and  Bittner 
values  are  substantially  higher  than  the  previous 
calculation  [4]  of  102  kcal,  although  the  large 
estimated  uncertainties  (-^  10,  -7  kcal)  place  all 
three  values  in  “agreement.”  In  the  remainder  of 
the  present  study  we  have  elected  to  use  the  lower 
limit  value  suggested  by  Stein.  Despite  the  com¬ 
mon  use  of  the  early  calculation  [4]  for  the  heat  of 
formation  of  C4H3,  we  believe  this  value  is 
significantly  in  error  due  to  the  updated  calculation 
of  Bittner  and  the  recent  estimate  [20]  of  81  kcal/ 
mol  for  the  heat  of  formation  of  butatriene,  which 
coupled  with  102  kcal/mol  for  C4H3  would  predict 
a  rather  low  (and  unlikely)  value  of  73kcaI/mol  for 
the  C-H  bond  strength  in  butatriene. 

PRELIMINARY  MODELING 

For  low-temperature  pyrolysis  (900-1400  K), 
Tanzawa  and  Gardiner  [1]  invoked  bimolecular 
initiation  followed  by  addition  of  H-atoms  to  C2HJ 
to  form  vinyl  radicals.  The  proposed  chain  mecha¬ 
nism  accounted  quite  well  for  the  experimental 
results  for  acetylene  decomposition.  However,  it 
did  not  describe  the  known  production  of  benzene 
and  higher  molecular  weight  species,  which  to¬ 
gether  accounted  for  more  than  50%  of  the 
decomposed  acetylene  (near  1000  K).  Frenklach 
and  coworkers  [5]  extended  this  mechanism  to 
describe  formation  of  heavier  species  including 
PAH  and  soot  formation  during  pyrolysis  of 
acetylene  at  higher  temperatures.  Colket  [16], 
subsequently,  used  a  similar  mechanism  for  com¬ 
parison  to  experimental  profiles  of  species  with 
molecular  weights  up  to  102.  Each  of  these  studies 
used  Reaction  1  or  la  as  the  initiation  step  with 
Arrhenius  parameters  as  defined  by  Tanzawa  and 
Gardiner  {E  =  45.9  kcal,  A  =  2  x  10'^  cm’ 
mol"'  sec"') 

Based  on  these  studies,  a  simplified,  low- 
temperature  sequence  can  be  written  (Table  I). 
This  sequence  differs  from  that  proposed  by 
Tanzawa  and  Gardiner  [1]  in  that  Reaction  la  is 
assumed  to  initiate  the  reaction  and  vinyl  recombi¬ 
nation  was  found  to  dominate  radical  termination. 
In  the  earlier  modeling  [1],  the  high  stability 
assumed  for  C4H3  led  to  its  important  role  in 
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TABLE  1 


SimpUncd  Reaction  Sequence  For  Low  Temperature  Pyrolysis 
of  Pure  Acetylene 


CjHj  +  CjHj  «  n-C4H,  +  H 

(la) 

H  +  CjH.  «  CzH, 

(2) 

C,Hj  +  CjHi  ~  C4H4  +  H 

(3a) 

CjHj  +  CjH,  «  C4H4 

(4) 

termination  and  C4H3  +  H  dominated.  Termina¬ 
tion  by  vinyl  recombination  has  also  been  sug¬ 
gested  by  Ogura  [3].  Duran  et  al.  [9]  and 
Frenklach  et  al.  [5]  also  included  vinyl-vinyl 
recombination  in  their  models.  In  the  former 
study,  an  unusually  high  rate  constant  (10'^  cm’ 
mol"'  sec"')  was  used  which  partially  explains 
why  their  attempts  to  model  acetylene  pyrolysis 
failed.  In  the  latter  case,  the  modeling  was 
performed  at  temperatures  above  1400  K  where 
thermal  decomposition  of  vinyl  is  fast  so  alterna¬ 
tive  recombination  processes  can  compete.  A 
dominant  termination  step  was  not  identified  by 
Frenklach,  et  al. 

Performing  a  steady-state  analysis  on  the  reac¬ 
tions  in  Table  1  leads  to 


£/[C2H2] 

dt 


=  2A:3a 


1/2 

[C2H2]' 


+  2/:u[C2H2]=. 


(Eq.  A) 


in  agreement  with  the  often  quoted  experimental 
observation  of  second-order  decomposition  kinet¬ 
ics.  Assuming  a  long  chain,  the  second  term  in  Eq. 
A  can  be  ignored  and  the  overall  activation  energy 
is  Eov  =  £'3a  +  (£1,  -  Et)/2.  Assuming  £3^  = 
4,  £4  =  0  kcal/mol  and  taking  £|,  =  68  (in  accord 
with  the  reaction  endothermicity),  then  Eov  =  38 
kcal/mol,  which  agrees  with  measured  values. 
Substitution  of  46  kcal/mol  (the  previously  used 
value  of  £1)  for  £,,  gives  a  steady-state  £ov  of  27 
kcal/mol,  substantially  lower  than  experimental 
results. 

Despite  the  attractiveness  of  the  “updated”  £|„ 
our  attempts  to  model  experimental  results  of 
Munson  and  Anderson  [15]  were  not  fully  satisfy¬ 
ing.  The  mechanism  and  rate  constants  adopted 
were  similar  to  those  used  by  Colket  [16].  Princi¬ 
pal  differences  included  values  of  two  rate  con- 


M.  B.  COLKET  ET  AL. 


slants 

=  mol"'  sec"' 

and 


k4  =  2x  lO'^cm^  mol" '  sec”'. 

The  first  expression  is  denoted  /Tu'”*'’  because  it 
is  higher  than  the  maximum  value  one  would 
expect.  The  A-factor  is  obviously  much  too  high 
for  a  bimolecular  initiation  of  this  type.  The  value 
of  k4,  for  vinyl-vinyl  recombination,  is  more  than 
five  times  lower  than  that  typically  used  for  this 
reaction.  Thus,  despite  reasonable  agreement  of 
experimental  acetylene  decay  and  product  forma¬ 
tion  with  the  model,  these  results  must  be  viewed 
with  caution. 


INDICATIONS  OF  ACETONE 
INITIATION 

With  concern  about  the  above  results  and  yet  with 
the  evidence  for  the  existence  of  a  chain,  possible 
initiation  steps  were  reexamined.  Of  the  likely 
possibilities 

CzHz-CjH-t-H- 128  kcal  mol”', 

CzHz-t-CzHz-^CjH-l-CjHj-S?  kcal  mol”', 

C2H2  +  C2H2  -*  i-C4H j  +  H  —  60  kcal  mol "  * , 

only  the  third  has  a  sufficiently  low  endothermic- 
ity,  and  yet  it  also  would  be  expected  to  have  a  low 
A-factor  in  order  to  form  the  H2CCCCH  radical. 
This  reaction  remains  a  possibility  and  cannot  be 
ruled  out.  Further  consideration  of  the  initiation 
problem  led  us  to  consider  impurities  as  potential 
initiators.  Acetone,  stored  in  acetylene  bottles  as 
an  inhibitor  of  self-detonation,  is  a  well-known 
impurity  of  acetylene  and  coelutes  from  tanks  in 
concentrations  of  l%-20%  that  of  acetylene  (23). 
Because  the  vapor  pressure  of  acetone  is  about  one 
fourth  an  atmosphere  at  room  temperature  and 
peak  cylinder  pressures  (Matheson)  are  about  18 
atm,  nominal  acetone  concentrations  can  be  esti¬ 
mated  to  be  1.3%  of  the  value  of  acetylene.  When 
precautions  are  taken  to  reduce  the  concentration 
of  impurities,  acetone  concentrations  are  typically 
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(see,  for  example,  ref.  23)  decreased  to  approxi¬ 
mately  0.1%  of  the  acetylene  concentration.  As¬ 
suming  that  a  radical  initiation  rate  of 
[C2H2]^  is  needed  to  produce  sufficient  radicals 
for  acetylene  pyrolysis,  then  this  rate  can  be 
compared  to  the  measured  [24]  rates  of  chain 
initiation  by  acetone 

CHjCOCHj-CHj  +  CHjCO,  (5) 

CH3CO-CH3  +  CO.  (6) 


Assuming  [acetone]/[C2H2]  =  0.001,  one  ob¬ 
tains 


ki  [acetone] 


=  0.004  X  cal/mol/RT)  [CjHj  ]  "  ' , 


with  Szwarc  and  Taylor’s  expression  [24]  for  k^ 
(  =  2.4  X  10''*  exp(- 72  kcal/mol/RT)  sec” ').  At 
20%  C2H2  and  1000  K  (Munson  and  Anderson’s 
conditions),  this  ratio  is  22  while  at  [C2H2]  = 
2.35  X  10'^  mol/cm’  and  1200  K  (Ogura),  the 
ratio  is  3.2.  Because  this  ratio  is  greater  than  1  at 
both  experimental  conditions,  radical  initiation  by 
thermal  decomposition  of  acetone  dominates  over 
that  of  acetylene. 

Considering  that  Ati,'"*'’  is  substantially  higher 
than  can  be  expected  from  thermodynamic  argu¬ 
ments  and  that  the  radical  initiation  rate  from 
acetone  is  higher  than  the  required  rate  of  initia¬ 
tion,  it  seems  quite  reasonable  to  expect  that 
acetone  plays  a  role  in  the  initiation  process.  For 
reference,  it  is  instructive  to  examine  the  ratio 
Ats  [acetone]/Ar|[C2H2j^  when  AT)  is  assigned  the 
value  determined  by  Tanzawa  and  Gardiner  [1]. 
This  ratio  is  then  0.13  and  0.12  for  the  same 
conditions  cited  at  the  end  of  the  previous  para¬ 
graph.  These  values,  in  all  probability,  led  to  early 
conclusions  that  acetone  initiation  contributed 
negligibly  to  acetylene  pyrolysis.  This  difference 
in  conclusions  suggests  a  major  difference  in  the 
values  of  ki  and  A:,,  so  reasons  for  this  difference 
need  to  be  examined. 

The  absolute  value  of  k\  determined  by  Tan¬ 
zawa  and  Gardiner  is  nearly  two  orders  of  magni¬ 
tude  higher  than  the  A:i,  found  in  the  present 
analysis.  The  principal  reason  for  their  higher 
value  is  believed  to  be  caused  by  the  low  AH/  for 
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the  C4H3  radical  (-20  kcal/mol  lower  than  that 
used  for  n-C4H3  in  the  present  analysis).  The 
resultant  high  stability  of  C4H3  radical  as  used  in 
their  model  led  to  high  radical  concentrations  and 
high  rates  of  radical  termination  via  C4H3  +  H, 
the  reverse  of  Reaction  1 .  This  step  played  only  a 
minor  role  in  the  present  study;  rather  vinyl-vinyl 
recombination  dominated  termination.  Consistent 
with  the  interpretation  of  high  radical  concentra¬ 
tions  (predicted  by  the  Tanzawa  and  Gardiner 
model  near  1000  K)  is  they  very  low  rate  constant 
that  was  used  for  the  principal  chain  step  (Reaction 
3a).  At  1000  K,  their  rate  constant,  =  1.6  x 
is  three  order  of  magnitude  below 
that  (for  Reaction  3)  used  in  the  present  study. 
Although  Reaction  3  was  found  to  be  nearly 
equilibrated  in  the  present  study,  the  net  rate  of 
Reaction  3  was  still  more  than  ten  time  faster  than 
Reaction  3a. 

The  important  conclusion  from  this  cursory 
analysis  is  that  acetone  decomposition  is  fast 
enough  to  influence  the  radical  concentrations 
during  acetylene  pyrolysis. 

ISSUES  RAISED  BY  ACETONE 
INITIATION 

Several  critical  issues  must  be  raised  if  initiation 
by  acetone  contributes.  First  of  all,  a  rapid 
mechanism  for  conversion  of  methyl  radicals  to  H- 
atoms  must  exist;  secondly,  the  “well-known” 
second-order  behavior  of  acetylene  must  be  ad¬ 
dressed,  specifically  to  see  if  this  reaction  order 
can  be  predicted  and  if  not,  why  not;  and  thirdly,  if 
acetone  initiation  dominates,  then  why  is  there 
such  good  experimental  agreement  (see  ref.  12)  on 
acetylene  decomposition  from  study  to  study  when 
initial  concentrations  of  acetone  can  be  expected  to 
vary? 

It  is  readily  shown  that  methyl  radicals  can  be 
rapidly  converted  to  H-atoms.  A  mechanism  for 
conversion  is  shown  in  Table  2.  Rate  constants 
were  obtained  from  the  literature  where  possible. 
The  principal  fate  of  methyl  near  1000  K  will  be  its 
addition  to  acetylene  (R7).  According  to  a  recent 
QRRK  analysis  [25]  the  resultant  adduct  should 
stabilize  under  these  conditions.  Using  rate  con¬ 
stant  (in  Tables  2  and  4)  and  thermodynamic 


parameters  nearly  identical  to  those  given  by  Dean 
and  Westmoreland  [25]  (see  Table  6)  the  rate  of 
CH3CHCH  isomerization  to  allyl  radical  relative 
to  its  rate  of  decomposition  to  methyl  acetylene 
and  H  is  approximately  5.  As  a  consequence  of  the 
stability  of  allyl,  its  decomposition  is  slow  and 
principally  adds  to  acetylene  to  form  C3H7  (R9). 
This  linear  radical  can  then  either  decompose  to 
reactants  (R(-8))  or  form  cyclopentadiene  plus 
atomic  hydrogen  (RIO).  Thus,  methyl  radicals 
produced  from  acetone  decomposition  will  be 
rapidly  converted  to  H-atoms.  Since  the  k5[ace- 
tone]/A:i,[C2H2]^  ratios  as  previously  calculated 
are  significantly  greater  than  1 ,  production  of  H- 
atoms  via  acetone  decomposition  appears  to  be 
sufficiently  fast  to  initiate  the  thermal  decomposi¬ 
tion  of  acetylene. 

At  higher  temperatures,  methyl  radical  addition 
to  acetylene  must  compete  with  methyl  recoribina- 
tion  (R17).  The  latter  process  can  become  impor¬ 
tant  at  high  temperatures  because  its  rate  is 
proportional  to  the  square  of  the  concentration  of 
methyl  radicals  and  the  rapid  decomposition  of 
acetone  (at  elevated  temperatures)  quickly  pro¬ 
duces  high  concentrations  of  methyl.  Also,  methyl 
can  extract  hydrogen  from  acetylene  (R16),  be¬ 
cause  at  high  temperatures  thermal  energy  is 
available  to  overcome  the  substantial  endothermic- 
ity  ( —  25  kcal/mol)  of  this  reaction.  In  either  high- 
temperature  case,  a  variety  of  alternative  se¬ 
quences  will  again  provide  H-atoms. 

The  order  n  for  the  overall  reaction  as  defined 


by 


t^[C2H2] 

dt 


=  A:ov[C2H2]" 


is  often  cited  at  2.  A  steady-state  analysis  using 
reactions  in  Table  1  and  Reactions  5-10  in  Table  2 
leads  to 


C/[C2H2]  flk, 

- =  2Ar3A  — 

dt  \  k. 


1/2 

[C2H2]  [acetone] 


or  assuming  [acetone]/[C2H2]  =  0.001,  then 


tf[C2H2] 

dt 


0.09k, ^  ^ 

^4 


1/2 


(Eq.  B) 


C-8 


M.  B.  COLKET  ET  AL. 


TABLE  2 

Proposed  Set  of  Reactions  and  Rate  Coefficients  for  Initiation  by  Acetone 
log  It  =  log  /J  +  /t  log  r  -  £/A’/r/2.303 


Forward 

Reverse 

rate  constant  ‘ 

rate  constant* 

Reactions* 

log  /t  n  E 

log  /I  n  E 

Ref 

5 

CHjCOCH,  =  CHjCO  +  CHj 

14.38 

0.0 

72.0 

11.55 

0.0 

-5.6 

24 

6 

CHjCO  =  CHj  +  CO 

12.48 

0.0 

16.7 

10.68 

0.0 

4.6 

e 

7 

CHj  +  C2H2  =  CHjCHCH 

11.79 

0.0 

7.7 

12.87 

0.0 

33.9 

48 

8 

CHjCHCH  =  CjH, 

13.15 

0.0 

36.0 

14.10 

0.0 

57.4 

25 

9 

CjHj  +  C2H2  =  CjH, 

12.00 

0.0 

8.0 

13.11 

0.0 

20.1 

est* 

10 

CjH,  -  c-CjtL  +  H 

10.30 

0.0 

5.0 

0.00 

0.0 

0.0 

CSl 

11 

H  +  CjR,  =  CHjCHCH 

12.76 

0.0 

3.1 

12.52 

0.0 

38.2 

49 

12 

H  +  ALLENE  =  CjH, 

12.60 

0.0 

2.7 

13.13 

0.0 

60.9 

49 

13 

2C,H5  =  CjH*  +  CjH4 

12.70 

0.0 

0.0 

13.20 

0.0 

30.5 

est 

14 

CjHj  +  C2H,  =  CjH. 

12.70 

0.0 

0.0 

15.62 

O.n 

82.5 

est 

15 

CjHj  +  H  =  CjHs 

13.60 

0.0 

0.0 

14.81 

0.0 

87.0 

est 

16 

CHj  +  CjHj  =  CR,  +  CjH 

12.40 

0.0 

35.0 

13.11 

0.0 

15.1 

est 

17 

2CHj  =  CjH* 

14.38 

-0.4 

0.0 

16.31 

0.0 

87.1 

e 

18 

2CHj  =  CjHj  +  H 

14.90 

0.0 

26.5 

16.92 

0.0 

16.3 

c 

19 

CjH«  +  CHj  =  CjH,  +  CR, 

-0.26 

4.0 

8.3 

14.71 

0.0 

25.1 

c 

20 

CjHj  =  CjH4  +  H 

13.30 

0.0 

39.7 

12.50 

0.0 

1.9 

e 

21 

CHj  +  CHjCOCHj  -  CH.  +  CHjCO  +  CH, 

11.60 

0.0 

9.7 

0.00 

0.0 

0.0 

50 

22 

H  +  CHjCOCHj  -  Hj  +  CHjCO  +  CHj 

13.28 

0.0 

6.4 

0.00 

0.0 

0.0 

49 

23 

CjHj  +  CHjCOCHj  -  CjH4  +  CHjCO  +  CHj 

12.48 

0.0 

6.4 

0.00 

0.0 

0.0 

est 

24 

CHjCO  =  CH;  +  CO 

14.00 

0.0 

71.0 

11.65 

0.0 

-4,4 

e 

25 

CHj  +  H  =  CHj  +  Hj 

14.86 

0.0 

15.1 

14.09 

0.0 

10.1 

51 

26 

CHj  +  CHj  =  CjHj  +  H 

13.30 

0.0 

0.0 

15.89 

0.0 

61.3 

51 

27 

2CHj  =  CjHj  +  Hj 

13.30 

0.0 

0.0 

15.65 

0.0 

128.7 

51 

28 

CHj  +  CjHj  =  CjH,  +  H 

12.26 

0.0 

0.0 

13.36 

0.0 

11.5 

C 

29 

2CjHj  -  C»H* 

12.70 

0.0 

0.0 

0.00 

0.0 

0.0 

est 

30 

CHjCO  +  H  =  CH,  +  CO 

12.85 

0.0 

3.0 

11.87 

0.0 

37.0 

r 

31 

CHjCO  +  CHj  =  CjHj  +  CO 

12.30 

0.0 

3.0 

13.34 

0.0 

26.7 

est 

*  represents  forward  and  reverse  directions  included  in  model.  ”  -  ”  represents  forward  direction  only  included  in  model. 

*  Units  for  A;  cm’/mol/sec.  Units  for  E:  kcal/mol. 

^  Rate  constant  from  ref.  29. 

''  Estimate  based  on  thermodynamics  and/or  analogous  reactions. 

'  Rate  constant  a  factor  of  3  lower  than  that  reported  in  ref.  29. 


This  Steady-State  result  that  n  =  3/2  appears  to  be 
a  major  drawback  of  the  present  proposal  regard¬ 
ing  initiation  by  acetone.  A  review  of  the  literature 
was  performed  in  an  attempt  to  find  references  in 
which  the  reaction  order  was  determined  or  which 
provide  data  from  which  the  order  could  be 
computed.  A  summary  of  reaction  orders  is 
presented  in  Table  3.  Although  there  is  substantial 
variation  in  the  reaction  order,  most  values  fall  in 
the  range  of  1. 5-2.0.  The  low-temperature  data 

c-n 


from  Silcocks  [26]  as  well  as  the  flow  reactor  data 
from  Palmer  and  Dormish  [27]  were  originally 
interpreted  assuming  simultaneous  homogeneous 
and  heterogeneous  processes.  Silcocks  found  the 
reaction  orders  to  be  2  and  1,  respectively, 
whereas  Palmer  and  Dormish  found  that  this 
assumption  was  consistent  with  their  data.  An 
attempt  has  been  made  to  reinterpret  the  latter 
experiments,  with  a  resultant  overall  order  ranging 
from  1.2  to  1.7.  Munson  and  Anderson’s  low- 


ACETONE-INITIATED  PYROLYSIS  OF 

TABLE  3 

Experimental  Overall  Reaction  Orders 


Apparatus 

Temp  (K) 

Order  (n) 

Silcocks  (26) 

static  reactor 

625-745 

l(het.).  2(hom.) 

Cullis  and  Franklin  [28] 

static  reactor 

983 

2 

Munson  and  Anderson  [IS] 

flow  reactor 

873,  923 
973-1073 

1.7  ±  0.2 

1.5  ±  0.1 

Palmer  and  Donnish  [27] 

flow  reactor 

1333-1528 

l(het.).  2(hom.) 

[1.2  to  1.7]'’ 

Ogura  [3] 

SPST 

1000-1670 

(-2,3) 

Colket  [16],  this  work 

SPST 

1100-1400 

1.6  ±  0.2 

Wu,  Singh,  and  Kern  (8) 

shock  tube 

1900-2500 

1.75  ±  0.2“ 

Aten  and  Greene  [52] 

shock  mbe 

1400-2500 

1.72  ±  0.1 

Cundall  et  al.  [53] 

shock  tube 

1500-2500 

I* 

Gay  et  al.  [54] 

shock  tube 

1600-2400 

not  2 

Towell  and  Martin  (55) 

flow  reactor 

1220-1400 

1.5“ 

*  Evaluated  in  present  study. 

*  Overall  order  for  formation  of  products. 


temperature  experiments  [15]  exhibit  significant 
scatter  in  plots  of  log(rate)  vs.  log[C2H2]  but  the 
order  appears  to  be  1.7  ±  0.2.  Their  higher 
temperature  experiments  clearly  give  n  =  1.5  ± 
0.1.  The  results  from  Colket  [16]  of  n  =  1.6  ± 
0.2  are  from  comparison  of  single-pulse  shock 
tube  experiments  at  4.9%  (unpublished).  3.5% 
(ref.  11),  and  0.2%  acetylene  (unpublished).  The 
order  obtained  by  Ogura  [3]  is  substantially  higher 
than  from  other  experiments.  Results  from  Cullis 
and  Franklin  [28]  at  983  K  give  beautiful  agree¬ 
ment  with  the  assumption  of  a  second-order 
reaction.  If  indeed  a  chain  mechanism  controls  the 
reaction,  however,  then  Cullis  and  Franklin’s  data 
need  reinterpretation.  They  pyrolized  pure  acety¬ 
lene  at  total  pressures  ranging  from  96  to  402  torr. 
Under  these  conditions,  some  critical  reactions, 
particularly  H  +  C2H2  -*  C2HJ,  can  be  expected 
to  exhibit  pressure  dependence  (see  Wamatz  [29]) 
and  will  lead  to  a  net  decrease  in  overall  reaction 
rate  with  decreasing  pressure.  If  pressure  depen¬ 
dence  is  included,  a  preliminary  analysis  indicates 
that  Cullis  and  Franklin’s  data  are  no  longer 
inconsistent  with  the  assumption  of  an  overall 


reaction  order  of  1.5.  Results  from  other  authors 
listed  in  Table  3  do  not  support  the  “well-known” 
reaction  order  of  2. 

Thus,  with  the  exception  of  the  very-low- 
temperature  data  of  Silcocks  and  the  value  from 
Ogura,  the  assumption  of  an  overall  reaction  order 
ot  1.5  as  suggested  by  acetone  initiation  is  not 
unreasonable.  In  fact,  the  experimental  data  seem 
to  support  a  value  of  1.5  as  well  as  or  perhaps 
better  than  the  order  of  2. 

If  acetone  is  the  initiator,  then  why  should 
experiments  performed  in  many  different  laborato¬ 
ries,  presumably  with  different  initial  acetone 
concentrations,  all  result  in  similar  overall  rate 
constants?  This  is  a  question  that  we  could  only 
partially  address.  A  variety  of  purification  tech¬ 
niques  have  been  used.  Often  the  level  of  impuri¬ 
ties  was  not  determined  or  not  determined  with 
sufficient  accuracy.  Based  on  our  own  experience, 
the  work  of  Hamins  et  al  [23],  and  reported  purity 
levels,  we  estimate  that  typical  purified  samples  of 
acetylene  contain  about  0.1  %  acetone,  with  a  total 
range  of  0.05%  to  0.2%.  This  variation  of  a  factor 
of  4  leads  to  a  factor  of  two  variation  in  overall 
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decomposition  rates  (because  of  the  half-order 
dependence  on  acetone  concentration)  and  is 
within  the  measured  experimental  uncertainty. 

DEVELOPMENT  OF  KINETIC  MODEL 

The  kinetic  model  used  in  this  analysis  is  based  on 
(1)  the  initiation  sequence  as  indicated  in  Table  2 
largely  using  published  rate  constants  and  (2)  an 
acetylene  pyrolysis  model  as  refined  by  Colket 
[16]  to  describe  results  obtained  in  a  single-pulse 
shock  tube.  Model  results  (using  a  program  based 
on  CHEMKIN  [30],  LSODE  [31],  and  isothermal, 
isobaric  assumptions)  were  compared  with  the 
flow  reactor  data  of  Munson  and  Anderson  [15]. 
The  mechanism  and  rate  constants  used  to  model 
the  flow  reactor  data  are  shown  in  Tables  2  and  4. 
Bimolecular  initiation  by  acetylene  (Rl)  was 
included  in  the  model  using  a  rate  constant  that  can 
be  expected  to  be  an  upper  limit.  In  the  modeling 
of  the  Munson  and  Anderson  data,  this  reaction 
always  had  a  negligible  role.  Alternatively,  Reac¬ 
tion  lb  was  not  included.  Assuming  kit,  = 
IQitg-M.ooo/RT,  jjjjj  reaction  could  be  a  few  times 
faster  than  Reaction  la  at  1000  K.  Nevertheless 
this  step  still  plays  a  negligible  role  in  radical 
initiation,  because  A:,b[C2H2]^  <  A:5[CH3- 
COCH3]  even  with  acetone  one-thousandth  the 
concentration  of  acetylene.  With  £ib  =  54  kcal, 
consistent  with  Bittner’s  thermodynamics  for  i- 
C4H3,  the  rate  of  Reaction  lb  is  still  an  order  of 
magnitude  less  than  that  of  Reaction  R5. 

A  few  rate  expressions  have  been  adjusted  in 
order  to  provide  better  fits  to  the  flow  reactor  data. 
Selection  of  rate  constants  for  modification  was 
accomplished  using  a  combination  of  a  reaction 
pathway  analysis  (with  computer  progiair.s  devel¬ 
oped  at  UTRC)  and  a  sensitivity  analysis  using 
CHEMSEN  [32].  For  each  rate  expression  modi¬ 
fied,  attempts  were  made  to  maintain  consistency 
with  the  SPST  modeling  results  [16]  by  adjusting 
both  the  A-factor  and  the  activation  energy.  In 
most  cases  this  was  possible  because  the  flow 
reactor  data  were  most  dependent  on  kinetics  over 
the  range  900-1000  K,  whereas  for  the  SPST 
results  the  range  was  1200-1500  K.  In  the  case  of 
the  competitive  sinks  for  the  n-C4H5  radical, 

(32) 


n-C4H5-l-C2H2-'C6H6-(-H,  (33) 

this  was  not  possible,  but  an  attempt  was  made 
instead  to  maintain  the  ratio  of  rate  constants  as 
determined  in  modeling  the  SPST  data.  Reaction 
33  has  been  simplified  from  the  previous  three- 
step  reversible  sequence  involving  intermediate 
formation  of  I-C6H7  and  c-CbHv,  since  a  QRRK 
analysis  by  Westmoreland  [33]  indicates  that  at  the 
conditions  of  the  present  study  the  reaction  rapidly 
proceeds  as 

C2H2  +  n-C4H 5  -  l-CftH  *  -  c-CftH  *  -  CsH^  +  H . 

Modeling  results  near  1000  K  using  the  irre¬ 
versible  one-step  or  the  three-step  reversible  se¬ 
quence  were  virtually  indistinguishable. 

Other  differences  between  the  present  and 
previous  [16]  model  include  rate  expressions  for 
H-atom  extraction  by  phenyl  or  vinyl  from  vinyla- 
cetylene.  These  were  equated  to  the  expression  for 
vinyl  plus  1,3-butadiene  suggested  by  Kiefer  et  al. 
[34].  Model  results  showed  little  sensitivity  to  rate 
constants  for  these  reactions.  Vinyl-vinyl  recom¬ 
bination  (R4)  was  taken  to  be  2  x  10'^  cmVmol/ 
sec  based  on  a  fit  to  the  Munson  and  Anderson 
data.  Rate  constants  fc  other  alkyl  radical  termi¬ 
nations  were  assumed  to  be  equal.  The  selected 
value  was  5  x  10‘- cmVmoI/sec.  After  Frenklach 
et  al.  [5],  the  two  isomers,  n-C4H3  and  i-C4H3, 
were  included;  however,  there  is  no  ready  sink  for 
i-C4H3  at  the  low  temperatures  of  the  Munson  and 
Anderson  study.  To  resolve  this  problem,  two 
reactions  were  included: 

n-C4H3  -  i-C4H3,  (61) 

i-C4H3-i-C4H4-C8Hb  +  H.  (66) 

The  first  of  these  was  found  to  be  negligible 
near  1000  K.  The  second  was  the  dominant  sink 
for  i-C4H3  while  affecting  phenylacetylene 
slightly.  The  first  reaction  was  assigned  a  rate 
constant  similar  to  a  value  [35]  for  the  i-C3H7 
n-C3H7  isomerization,  although  a  barrier  close  to 
the  strain  energy  in  cyclopropene  was  assumed. 
The  resultant  activation  energy,  £j,  =  53  kcal/ 
mol,  is  rather  high  (especially  if  tunneling  occurs, 
as  in  the  case  of  the  vinylidene-acetylene  isomer¬ 
ization)  so  we  tested  our  mechanism  with  a  much 


n-C4H5  H-t-C4H4, 


ACETONE-INITIATED  PYROLYSIS  OF  CjHj 

TABLE  4 

Proposed  Set  of  Reactions  and  Rate  Coefficictils  for  Acetylene  Pyrolysis  at  Low  and  Inlermedijte  Temperatures 

log  It  =  log  A  +  n  log  r  -  E/R/T/2.m‘ 


Reactions* 

Forward 

rate  constant 

log  A  n 

d 

£ 

Reverse 

rate  constant 

log  A  n 

If 

E 

Ref 

1 

2C2HJ  =  n-C4Hj  +  H 

13.00 

0.0 

67.0 

12.48 

0.0 

-2.0 

esi*^ 

2 

H  +  QHj  =  CjHj 

12.74 

0.0 

2.5 

13.01 

0.0 

43.7 

56 

3 

CjH,  +  CjHj  =  n-C4H5 

12.04 

0.0 

4.0 

13  88 

0.0 

38.4 

16 

4 

2C2H,  =  C4H* 

13.30 

0.0 

0.0 

16.71 

0.0 

101. 1 

pw^ 

32 

n-C4Hj  =  C4H4  +  H 

14.20 

0.0 

41.4 

13.60 

0.0 

3.0 

pw 

33 

n-C4Hj  +  C2H2  -  CJft  +  H 

12.65 

0.0 

10.0 

0.00 

0.0 

0.0 

pw 

34 

CjHj  +  C4H4  -  C»H«  +  H 

11.60 

0.0 

0.0 

0.00 

0.0 

0.0 

r 

35 

C4H4  +  n-C4H5  -  C,H,  +  H 

13.08 

0.0 

5.0 

0.00 

0.0 

0.0 

pw 

36 

C,H,  +  H  =  C,H,  +  Hj 

14.60 

0.0 

7.0 

13.65 

0.0 

12.0 

16 

37 

C.H,  =  C,H*  +  H 

12.85 

0.0 

45.0 

12.95 

0.0 

0.7 

esi 

38 

C2H4  +  C*Hi  =  C,H,  +  H 

11.27 

0.0 

2.1 

12.38 

0.0 

1.6 

/ 

39 

C»H*  +  H  =  H2 

14.40 

0.0 

16.0 

12.69 

0.0 

9.8 

29 

40 

CjHj  +  CJfj  =  C,H*  +  H 

12.00 

0.0 

4.0 

13.88 

0.0 

6  0 

t 

41 

CjH,  +  C.H4  =  C2H4  +  n-C4Hj 

13.30 

0.0 

14.5 

13.23 

0.0 

13.4 

34 

42 

C4H4  +  C^,  =  +  n-C4H, 

13.30 

0.0 

14.5 

13.60 

0.0 

18  9 

34 

43 

1-C»H,  =  n-C4H,  +  CjHj 

14.00 

0.0 

36.0 

11.69 

0.0 

1.0 

16 

44 

C»H,  =  l-C»Hj 

13.54 

0.0 

65.0 

9.91 

0.0 

1.1 

16 

45 

CjH,  +  =  CjH4  +  CJH, 

13.30 

0,0 

14.5 

12.94 

0.0 

9.1 

f 

46 

C:H,  +  H  =  C:H4 

13.60 

0.0 

0.0 

15.55 

0.0 

105,1 

est 

47 

C!H4  +  H  =  CjH,  +  H: 

14.84 

0.0 

14.5 

13.49 

0.0 

13.8 

29 

48 

CjH,  +  H  =  H,  +  C.Hj 

13.00 

0.0 

0,0 

13.34 

0.0 

63.1 

29 

49 

C4H4  +  H  =  n-C4H,  +  Hj 

13.90 

0.0 

14.5 

12.49 

0.0 

12.7 

29 

50 

COL  =  C^,  +  H 

15.70 

0.0 

107.9 

13.39 

0  0 

-2.6 

57 

51 

C*H,  +  CJL  =  C,jH,o  +  H 

11.80 

0.0 

11.0 

13.04 

0.0 

8.7 

e 

52 

2n-CiH5  =  C|H)o 

12.70 

0.0 

0.0 

16.92 

0.0 

106.9 

CSl 

53 

n-C4H,  +  H  =  C4H* 

13.60 

0.0 

0.0 

15.44 

0.0 

107.9 

cst 

54 

n-C4H,  ^  CjH,  =  C»H, 

12.70 

0.0 

0.0 

16.33 

0.0 

104.0 

est 

55 

C»H«  +  H  =  n-C4H,  +  Hj 

14.00 

0,0 

14.5 

12.76 

0.0 

10.9 

34 

56 

C4H«  +  CjH,  =  n-C4H,  +  C2H4 

13.30 

0.0 

14.5 

13.41 

0.0 

11.7 

f 

57 

CJI,  +  C:H,  =  C.H, 

12.70 

0,0 

0.0 

15.76 

0.0 

104.5 

cst 

58 

i-C4H,  +  H,  =  CjHj  +  CjH, 

10.70 

0.0 

20.0 

10.58 

0.0 

17.8 

16 

59 

i-C4H,  =  C4HJ  +  H 

12.00 

00 

49.0 

12.05 

0.0 

0.6 

16 

60 

n-C4H,  =  C4HJ  +  H 

12,60 

0.0 

40.0 

12.95 

0.0 

-0.4 

16 

61 

n-C.H,  =  i-C.H, 

13.00 

0.0 

53.0 

13.30 

0.0 

61.0 

cst 

62 

1-C4H,  +  H  =  C4HJ  +  Hj 

13.00 

0.0 

0  0 

13.65 

0.0 

56.0 

16 

63 

C4H4  +  H  =  i-C4H,  +  Hj 

14.49 

0.0 

14.5 

13.38 

0.0 

20  7 

16 

64 

C»H,  +  C4H4  =  C*H*  +  i-C4H, 

13.30 

0.0 

14.5 

13.90 

0.0 

26.9 

f 

65 

C2HJ  +  C4H4  =  C2H4  +  i-C4Hj 

13.30 

0.0 

14.5 

13.54 

0.0 

21.4 

f 

66 

C4H4  +  i-C4H,  -  C,H*  +  H 

11.60 

0.0 

0.0 

0.00 

0.0 

0.0 

cst 

67 

Hj  +  M  =  2H  +  M 

12.35 

-0  5 

92.5 

9.98 

0.0 

-  13.1 

29 

*  "  =  "  represents  forward  and  reverse  directions  included  in  model.  "  -  "  represents  forward  direction  only  included  in  model. 

*  Units  for  A:  cmVmol/sec,  Units  for  E:  kcaJ/mol. 

'  Estimate  based  on  thermodynamics  and/or  analogous  reactions. 

Rate  evaluated  from  present  work. 

'  See  citations  in  ref.  16.  (Reverse  rate  constants  may  differ  slightly  due  to  differences  in  the  thermodynamics  and  the  temperature 
range  of  the  fit  for  the  reverse  rate.) 

/  Rate  constant  half  of  that  cited  in  ref.  16 

*  See  text. 
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lower  activation  energy  of  5  kcal/mol.  Virtually 
no  change  was  observed  in  the  computations  for 
the  Munson  and  Anderson  conditions.  Despite  the 
insensitivity  of  our  calculations  to  this  rate  con¬ 
stant.  we  retain  the  high  £51  because  the  isomeriz- 
tion  process  is  more  likely  analogous  to  the 
rearrangement  of  the  vinylidene-acetylene  anions 
(because  of  the  unpaired  electron).  Calculations 
[22]  for  the  latter  process  indicate  a  barrier  of  45- 
50  kcal/mol  in  qualitative  agreement  with  the  first 
estimate  for  £41.  Reaction  66  has  been  proposed 
[36]  as  an  alternative  mechanism  for  forming 
phenylacetylene.  It  is  a  multistep  process  requiring 
a  1,4  H-atom  shift,  cyclization,  and  loss  of  an  H- 
atom  and  was  assigned  a  rate  expression  identical 
to  that  of  Reaction  4  which  requires  a  similar 
transformation.  Neglecting  the  back  reaction  is 
justified  since  the  reverse  process  is  slow  at  these 
temperatures.  A  QRRK  analysis  [37]  has  been 
performed  for  Reaction  66  and  indicates  that  near 
1000  K  approximately  50%  of  the  C4H4  +  i-C4H3 
adducts  will  go  directly  to  products  prior  to 
stabilization  to  intermediates.  Therefore  elimina¬ 
tion  of  the  intermediate  steps  is  reasonable. 

A  third  possible  sink  for  i-C4H3  not  included  in 
the  modeling  is  the  addition  of  i-C4H3  to  acety¬ 
lene,  cyclization  to  form  a  five-membered  ring  and 
then  isomerization  to  phenyl. 

Third-body  effects  have  been  assumed  to  be 
small  for  this  analysis.  The  greatest  deviation  from 
this  assumption  will  occur  for  the  H-atom  addition 
to  acetylene.  Fall-off  curves  given  by  Wamatz 
[29]  indicate  that  the  k/k^  for  this  reaction  is 
approximately  0.5  at  1  atm  of  argon  and  1000  K. 
In  the  present  study,  20%  of  the  gas  is  acetylene,  a 
more  efficient  third  body,  suggesting  the  k/k„  is 
larger  than  0.5;  therefore,  neglecting  pressure 
effects  is  believed  to  result  in  only  a  small  error.  A 
reduction  of  ki  by  a  factor  of  two  results  in  a 
reduction  of  the  overall  decomposition  rate  of 
about  5%-10%. 

Some  modeling  was  also  performed  for  compar¬ 
ison  to  shock  tube  results.  For  these  calculations, 
the  kinetic  mechanism  was  expanded  to  include 
ethynyl  and  related  radicals  (Table  5). 

Thermodynamic  values  are  reported  in  Table  6 
and  are  the  same  as  those  used  previously  [16], 
except  for  those  of  the  C4H3  isomers,  phenyl 


radical  and  the  Cg-species,  which  were  obtained 
from  Stein  [6].  The  latter  are  identical  to  those 
used  in  the  study  by  Frenklach  et  al.  [5]. 

KINETIC  MODEL  VS.  THEORY 

In  Figs.  2-5,  calculations  of  acetylene,  benzene, 
vinylacctylene,  and  ethene  yields  are  compared  to 
experimental  profiles  as  measured  by  Munson  and 
Anderson  [15]  in  a  flow  reactor.  For  these 
calculations  it  was  assumed  that  the  acetone- 
acetylene  ratio  was  0.001.  Overall  the  agreement 
is  quite  satisfactory.  Benzene  is  overpredicted  at 
long  times  for  the  highest  temperature  run,  pre¬ 
sumably  because  of  the  conversion  of  benzene  to 
higher-molecular-weight  species.  Vinylacetylene 
is  undcrpredicted  at  the  low  temperatures.  Hydro¬ 
gen,  not  shown,  is  underpredicted  by  about  a 
factor  of  5-10.  Although  substantial  methane 
formation  was  observed,  a  model  for  methane 
formation  was  not  included.  Methane  produced 
from  methyl  directly  or  via  odd-carbon  species 
derived  from  acetone  could  only  account  for  a 
very  small  ponion  of  the  methane  observed. 
Heterogeneous  catalysis  seem  the  most  likely 
source  of  methane  although  gas-phase  reactions 
are  possible.  For  example,  consider  addition  of  the 
i-C4H5  radical  (H2CCHCCH2)  to  acetylene  fol¬ 
lowed  by  cyclization  to  a  five-membered  ring  and 
isomerization  to  methylcyclopentadienyl.  Subse¬ 
quent  reactions  involving  methylcyclopentadienyl 
may  then  form  methyl  radicals  and/or  methane. 
The  i-C4H5  radical  is  generally  ignored  in  kinetic 
modeling  but  may  play  a  significant  role  because  it 
.should  be  formed  readily  via  H-addition  to  vinyla- 
cetylenc. 

The  formation  of  polyaromatic  hydrocarbons  is 
not  included  in  the  present  model.  When  reactions 
including  the  formation  of  species  up  to  and 
including  pyrene  are  added  to  the  reactions  in 
Tables  2  and  4,  computed  results  for  acetylene, 
benzene,  and  vinylacetylene  fit  the  experimental 
data  almost  perfectly  for  the  Munson  and  Ander¬ 
son  data  of  1073  K.  There  is  relatively  little 
change  to  the  model  predictions  for  the  lower 
temperature  runs.  Reactions  involving  the  high- 
molecular-weight  species  are  similar  to  those 
propo.scd  by  Frenklach  et  al.  [5].  Specific  reac- 
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TABLE  5 

Proposed  Set  of  Reactions  and  Rate  Coefneients  for  Higb-Tempcralure  Acetylene  Pyrolysis 
log  *  =  log  A  +  n  log  r  -  E/R/T/2.m 


Forward  Reverse 

Rale  Constant'’  Rate  Consunt 


Reactions* 

log  A 

n 

£ 

log  A 

n 

£ 

Ref 

68 

CjHi  +  M  =  CjH  +  H  +  M 

16.62 

0.0 

107.0 

15.21 

0.0 

-  18.2 

C 

69 

CjH  +  C4H4  =  C2H2  +  i-C.H, 

13.60 

0.0 

0.0 

13.29 

0.0 

27.0 

c 

70 

n-C4H3  =  CjHj  +  CjH 

14.30 

0.0 

57.0 

13.42 

0.0 

0.8 

16 

71 

n-C4Hi  +  H  =  C4H2  +  H, 

12.00 

0.0 

0.0 

12.95 

0.0 

64.0 

16 

72 

i-CHj  +  H  =  C4H4 

13.78 

0.0 

0.0 

15.49 

0.0 

98.2 

est 

73 

CjH  +  Hj  =  H  +  CjH: 

12.85 

0.0 

0.0 

13.65 

0.0 

20.8 

c 

74 

CjH  +  C2H2  =  C4H2  +  H 

13.60 

0.0 

0.0 

14.84 

0.0 

15.8 

c 

75 

C4H2  =  C4H  +  H 

14.89 

0.0 

120.0 

13.25 

0.0 

1.9 

c 

76 

C2H  +  C4H2  =  C»H2  +  H 

13.60 

0.0 

0.0 

15.04 

0.0 

15.5 

c 

77 

C4H  +  C2H2  =  C*H2  +  H 

13.30 

0.0 

0.0 

14.97 

0.0 

8.5 

16 

78 

CtH,  =  C4H  +  H 

14.89 

0.0 

120.0 

13.02 

0.0 

5.2 

c 

79 

C*H  +  H2  =  H  +  C4H2 

13.30 

0.0 

0.0 

14.34 

0.0 

13.8 

c 

80 

C*H  +  H2  =  H  +  c:»Hj 

13.30 

0.0 

0,0 

14.56 

0.0 

10.5 

c 

81 

C2H  +  CJL  =  CtH,  +  C2H2 

13.30 

0.0 

0.0 

12.40 

0.0 

14.7 

c 

82 

C4H  +  CtH*  =  C»H,  +  C4HJ 

13.30 

0.0 

0.0 

12.63 

0.0 

7.6 

< 

83 

CjHj  +  C4H2  =  C4H4  +  C2H 

13.48 

0.0 

23.0 

13.47 

0,0 

3.1 

16 

84 

C4H4  =  C2H  +  C2H, 

16.00 

0.0 

105.0 

13.37 

0,0 

-  16.2 

16 

•  "  =  ”  represents  forward  and  reverse  directions  included  in  model. 

‘  Units  for  A:  cmVmol/sec.  Units  for  E:  kcal/mol. 

'  See  citations  in  ref.  16.  (Reverse  rate  constants  may  differ  slightly  due  to  differences  in  the  thermodynamics  and  the  temperature 
range  of  the  fit  for  the  reverse  rate.) 

^  Estimate  based  on  thermodynamics  and/or  analogous  reactions. 
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Fig.  2.  Comparison  of  data  and  model  predictions  of  acetylene  decay.  Symbols  and  solid  lines 
are  data  and  best  fits  from  Munson  and  Anderson  (ref.  15).  Dotted  lines  are  model  predictions 
from  the  present  study  where  lOtX)  ppm  acetone  in  acetylene  is  assumed. 
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TABLE  6 


Selected  Thermodynamics  at  300  K* 


Symbol 

Name/structurc 

s* 

eu 

AH, 

kcal/mol 

C12H 10 

biphenyl 

93.6 

43.6 

C,H,o 

octatetraene 

90.4 

53.2 

C.H, 

styrene 

82.6 

35.3 

C,H, 

CJijCCHj 

83.8 

83.2 

C.H* 

phenylacetylene 

76.4 

75.2 

CtH. 

hexatriene 

79.4 

39.6 

benzene 

65.2 

19.8 

C*H5 

phenyl 

69.1 

78.2 

C*H2 

iriacetylene 

71.1 

169.7 

C*H 

C»H 

74.3 

233.2 

CA 

1 ,3-butadiene 

66.6 

26.1 

n-C.Hs 

HiCHCHCH 

69.1 

82.5 

CH4 

vinylacetylene 

66.1 

69.4 

i-C.H, 

HCCCCHj 

68.1 

116.1 

n-CJtj 

hccchCh 

69.4 

124.1 

C.Hj 

diacetylene 

59.9 

111.7 

C.H 

C.H 

62.8 

179.0 

CjHj 

vinyl 

54.5 

65.7 

C,H 

ethynyl 

49.6 

128.5 

C,H, 

1 ,4-pentadiene 

76.6 

18.7 

C5H, 

1,4-pentadien-l-yl 

79.3 

78.2 

cyclopentadiene 

64.6 

32.0 

CH)CO 

acetyl 

62.0 

-6.0 

CjH, 

ethyl 

54.6 

26.4 

CHjCHCH 

ch,chCh 

67.5 

61.6 

C,H, 

allyl 

62.1 

39.4 

C,H4 

methylacetylene 

59.3 

44.4 

C,H, 

propargyl 

60.0 

81.5 

CH: 

methylene 

46.4 

92.4 

CHjCO 

ketene 

57.8 

-  12.4 

'  Thermodynamics  of  species  not  listed  in  this  table  are 
essentially  identical  to  those  listed  in  ref.  35. 


lions  and  rate  constants  have  been  recently  de¬ 
scribed  (36]. 

Model  calculations  have  also  been  compared  to 
data  obtained  in  a  single-pulse  shock  tube  by 
Colket  [16]  at  3.5%  and  4.9%  acetylene  in  5  and  8 
atm  of  argon,  respectively.  Experimental  proce¬ 
dures  for  the  unpublished  4.9%  data  were  similar 
to  those  previously  described  [16].  The  acetone 
concentrations  were  equated  to  the  measured 
values,  that  is,  0.1%  and  0.2%  the  concentration 
of  acetylene  for  the  3.5%  and  4.9%  data,  respec¬ 
tively.  Calculations  were  performed  using  a  shock 
tube  code  modified  to  account  for  quenching 


processes.  Again  satisfactory  agreement  (compa¬ 
rable  to  that  achieved  in  ref.  16)  is  obtained  for 
most  species,  including  the  minor  products, 
styrene  and  ethylene.  The  main  descrepancy  is  that 
the  model  underpredicts  product  formation  at  low 
temperatures.  Of  particular  interest  is  Fig.  6,  in 
which  only  odd-carbon  products  derivable  from 
the  impurity  acetone  are  shown  for  the  4.9%  series 
of  runs.  Also  plotted  is  the  sum  of  the  concentra¬ 
tions  of  acetone  plus  half  of  all  species  containing 
an  odd  number  of  carbons.  Acetone  produces  two 
methyl  radicals,  each  of  which  could  produce  an 
odd-carbon  species.  The  fact  that  this  sum  is 
essentially  constant  over  the  range  of  temperature 
of  experiments  lends  strong  support  to  a  major 
thesis  of  this  work,  namely  that  acetone  produces 
methyl  radicals  that  convert  to  H-atoms  via  addi¬ 
tion  to  acetylene,  whereas  ethane  formation  via 
methyl  recombination  is  small  compared  to  methyl 
addition  to  acetylene.  Model  predictions  are  not 
shown  in  Fig.  6  because  prediction  of  the  odd- 
carbon  products  is  generally  poor.  This  limitation 
is  probably  a  direct  result  of  the  fact  that  there  are 
very  little  experimental  data  on  these  species,  and 
many  of  the  rate  constants  in  Table  2  were 
estimated. 

Ogura  also  obtained  single-pulse  shock  tube 
data  at  5  %  acetylene  in  about  two  atmospheres  of 
argon.  Model  predictions  are  compared  to  the 
Ogura  [3]  data  in  Fig.  7  for  vinylacetylene, 
diacetylene,  and  hydrogen.  Again  comparison 
between  experiments  and  theory  seem  quite  satis¬ 
factory.  For  the  model  calculations,  dwell  times  of 
1  msec  were  assumed  and  quenching  rates  compa¬ 
rable  to  those  measured  by  Colket  [16]  were  used. 
The  acetone-acetylene  ratio  was  assumed  to  be 
0.001.  Attempts  to  fit  Ogura’s  experimental  pro¬ 
files  of  methane,  methylacetylene,  allene,  and 
ethylene  for  experiments  with  10%  acetylene  in 
argon  were  also  made.  The  agreement  was  good 
for  methane  and  satisfactory  for  ethylene  except  at 
low  temperatures  when  the  model  significantly 
underpredicted  the  experimental  data.  The  model 
overpredicts  the  Cs-hydrocarbons  by  about  a 
factor  of  two  to  three.  This  difference  does  not 
seem  unreasonable  considering  the  general  lack  of 
knowledge  regarding  pyrolysis  of  odd-carbon, 
acetylenic  species. 
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TIME  (SEC)  time  (sec) 

Fig.  3.  Comparison  of  data  and  model  predictions  of  vinylacetylene  formation  and  decay.  (See 
caption  of  Fig.  2.) 

VINYLIDENE  VS.  CHAIN  radicals  during  acetylene  pyrolysis.  Calculations 

for  their  experimental  conditions  indicate  that 
Perhaps  the  most  significant  pieces  of  experimen-  during  acetylene  pyrolysis  the  total  radical  concen- 

tal  evidence  against  a  radical  chain-mechanism  are  (rations  were  always  substantially  less  than  one 

the  results  by  Hou  and  Anderson  [38]  and  more  part  per  million.  It  appears  unlikely  that  their 

recently  by  Duran  et  al.  [10].  The  earlier  study  facility  had  the  sensitivity  to  observe  species  at 

involved  a  molecular-beam,  mass-spectrometric  such  low  concentrations,  whereas  methyl  radicals 

analysis  of  acetylene  pyrolysis  products  from  a  readily  produced  from  the  test  compounds  would 

flow  reactor  near  1000  K.  The  usefulness  of  their  be  expected  to  achieve  substantial  concentrations, 

facility  was  "proved”  by  its  ability  to  detect  We  think  that  these  early  experiments  do  not 

methyl  radicals  during  pyrolysis  of  tetramethyl  disprove  a  radical  chain  mechanism, 

lead  and  di-t-butylperoxide.  They  detected  no  Recent  pyrolysis  experiments  by  Duran  et  al. 
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Fig.  5.  Comparison  of  data  and  model  predictions  of  ethylene  formation.  (See  caption  of  Fig. 
2.) 


[10]  were  also  performed  in  a  flow  reactor  with 
mass  spectrometric  analysis.  Additions  of  benzene 
and  toluene  showed,  respectively,  no  inhibition 
and  an  acceleration.  The  acceleration,  unexplained 
by  the  vinylidene  mechanism,  is  easily  explained 
via  chain  processes  since  initiation  via  toluene 
(C7H8  -*  C7H7  +  H)  using  any  of  the  recently 
proposed  rate  expressions  (see  refs.  39  and  40  and 
citations  contained  therein)  can  be  shown  to  be 
substantially  faster  than  initiation  by  acetone.  In 
addition,  Duran  et  al.  observed  a  variety  of 


vinylbenzenes  when  benzene  and  acetylene  were 
copyrolyzed.  They  interpreted  these  results  as 
evidence  for  insertion  of  the  carbene  (vinylidene) 
into  the  benzene  ring.  Alternatively,  these  vinyl- 
benzenes  may  also  be  the  result  of  the  displace¬ 
ment  (C2H3  +  CftHs  H  -F  C6H5C2H3),  which 
is  similar  to  the  known  [41]  methyl-substitution  in 
the  benzene-toluene  system.  Copyrolysis  of  to¬ 
luene  and  acetylene  produced  principally  a  species 
at  mass  118,  interpreted  to  be  methyl  styrene. 
Instead  this  species  may  be  phenylpropene  formed 
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Fig.  7.  Comparison  of  data  and  model  for  vinylacetylene,  diacetylene,  and  hydrogen  for 
pyrolysis  of  5%  acetylene.  Single-pulse  shock  tube  data  is  from  Ogura  (ref.  3).  Solid  lines  are 
model  predictions  assuming  KXX)  ppm  acetone  in  the  acetylene. 


by  acetylene  addition  to  the  stable  ben2yl  radical. 
Consequently,  we  question  their  conclusion  of 
proof  for  the  vinylidene  mechanism. 


OVERALL  REACTION  ORDER  AND 
ACTIVATION  ENERGY 


Reaction  order  was  determined  from  model  pre¬ 
dictions  by  varying  the  initial  concentration  (while 
holding  temperature  constant  at  973  K).  Plots  of 
log(t/[C2H2]/</f)o  vs.  log[C2H2]o  yielded  an  over¬ 
all  reaction  order  of  1.43  ±  0.10.  The  initial  rate 
d[C2H2]/d(]Q  was  taken  to  be  the  value  immedi¬ 
ately  after  a  short  induction  period  (  —  0.2  sec  at 
20%  acetylene,  873  K).  Reaction  orders  were  also 
determined  from  plots  of  predicted  log(d[C2H2]/ 
dt)  vs.  [C2H2]  as  a  function  of  time  over  the 
temperature  range  873-1073  K.  These  calcula¬ 
tions  resulted  in  a  reaction  order  of  1.43,  1.66, 
and  1.67  at  973,  1023,  and  1073  K,  respectively. 
These  values  are  consistent  with  the  overall  order 
of  1.5  obtained  from  the  steady-state  analysis. 
Assuming  this  latter  value,  an  overall  rate  expres¬ 
sion  was  determined  using  the  model  prediction 
for  Munson  and  Anderson’s  conditions.  The  rate 
expression  determined  was 


dlC2H2] 

dt 


=  A:,„(C2H2]'^ 


where 


i(rov=  10"e"'‘°*'‘'*""’®''^(cmVmol)''^  sec“‘ 


at  873-1073  K.  This  expression  can  be  compared 
to  that  obtained  from  the  simplified  steady-slate 
expression  (sss),  Eq.  B.  Equating  k-i^  to  ^3  and 
using  other  rales  from  Tables  2  and  4,  this 
simplified  expression  becomes 


t^IC2H2]\>» 
dt  /o* 


=  ^“[C2H2]'^ 


where 


/t“=  3.4x  10"e'^'‘'*""’°‘'’^^'(cmVmol)''^  sec " 

The  lower  A-facior  from  the  detailed  model  is 
due  to  a  combination  of  the  reverse  of  R3, 
additional  termination  reactions  (e.g.,  R54),  and 
the  fact  that  some  acetone  decomposes  via  a  chain 
process  rather  than  R5  (see  reactions  21-23  in 
Table  2). 


HIGHER  TEMPERATURE  MODELING 

Preliminary  modeling  has  been  performed  for 
shock  tube  conditions  at  higher  temperatures.  The 
scenario  regarding  initiation,  termination,  and  the 
fate  of  acetone  and  methyl  radicals  changes 


C-18 


M.  B.  COLKET  ET  AL. 


somewhat.  The  changes  include  a  greater  fraction 
of  acetone  decomposing  via  radical  attack.  This 
leads  to  the  production  of  ketene  which  also  is  an 
effective  radical  initiator.  The  resulting  sequence. 


H  +  CH3COCH3  «  CH2COCH3  +  Hz.  (22a) 
CH2COCH3  «  CH2CO  +  CH3.  (22b) 

CH2CO  «  CH2  +  CO,  (24) 

CHz  +  C2H2  C3H3  +  H ,  (28) 

C3H3  +  C3H3  «  C5H6.  (29) 


leads  to  the  net  production  of  only  one  methyl 
radical  (the  H-atom  is  regenerated)  compared  to 
the  two  radicals  from  initiation  directly  through 
acetone  decomposition  (Reaction  5).  This  se¬ 
quence  is  very  dependent  on  the  uncertain  high- 
temperature  chemistry  of  methylene  (CH2)  and 
propargyl  (C3H3)  radicals.  In  addition  to  recombi¬ 
nation,  the  latter  can  add  to  acetylene  to  form  a 
C5H5  radical  whose  stability  is  enhanced  by 
resonance.  Other  potential  complicating  features 
are  that  some  species  formed  as  by-products  of 
acetone  decomposition,  e.g.,  cyclopentadienes, 
have  weak  C-H  bonds  because  of  the  resonantly 
enhanced  stability  of  the  hydrocarbon  radical  (in 
this  case  cyclopentadienyl).  Consequently,  even 
though  acetone  may  be  destroyed  early  during 
reaction,  some  of  the  products  will  continue  to 
provide  radicals  and  maintain  the  decomposition 
of  acetylene.  At  high  temperatures,  vinyl  recombi¬ 
nation  no  longer  dominates  radical  termination. 
Methyl  radicals  and  radical  adducts  formed  by 
addition  to  acetylene  contribute  to  chain  termina¬ 
tion. 

Recent  time-of-flight/shock  tube  results  have 
also  been  obtained  by  Wu  et  al.  (8]  at  much  higher 
temperatures.  They  modeled  their  own  as  well  as 
other  high  temperature  data  using  the  bimolecular 
initiation  Reaction  lb  with  an  endothermicity  of  46 
kcal/mol  and  102  kcal/mol  for  the  heat  of  forma¬ 
tion  of  i-C4H3,  the  assumed  product.  Alterna¬ 
tively,  assuming  an  impurity  of  1000  ppm  acetone 
in  the  acetylene,  we  attempted  to  model  the  data 
using  the  reactions  in  Tables  2,  4,  and  5.  The 
experimental  decomposition  was  severely  under¬ 
estimated.  An  examination  of  the  model  results 


indicated  that  methyl  radicals  added  to  acetylene  to 
eventually  form  methylacetylene  and/or  allene  (via 
Reactions  7,  8,  11,  and  12).  Inclusion  of  Wu  and 
Kern’s  [42]  reactions  and  rate  coefficients  for  the 
thermal  decomposition  of  C4H3  resulted  in  a  closer 
description  of  the  T.O.F.  data.  A  very  close  match 
to  their  data  was  obtained  by  increasing  the  rate 
constant,  kT,  by  a  factor  of  3  and  assuming  the 
products  of  the  reaction  to  be  C3H4  -f-  H.  Using 
these  revisions,  acetylene  decay  can  be  predicted 
and  is  compared  to  experimental  data  in  Fig.  8a 
for  a  range  of  temperatures.  A  demonstration  of 
the  effect  of  different  assumed  levels  of  acetone  is 
shown  in  Fig.  8b  for  the  production  of  diacetylene. 
Although  a  pressure  dependence  is  expected  for 
R5,  this  was  not  included  in  the  calculations  shown 
in  Fig.  8.  This  omission  is  of  negligible  concern 
because  acetone  decomposes  much  faster  than  the 
time  scale  of  the  experiment  and  the  computed 
results  are  insensitive  to  the  value  of  k^.  Based  on 
uncertainties  in  thermochemistry,  the  acetone  con¬ 
centration,  and  the  revisions  to  Reaction  7,  we 
believe  that  the  good  comparison  between  the 
experiment  and  model  should  be  considered  quali¬ 
tative,  but  strongly  suggestive  of  the  fact  that 
impurities  also  effect  experiments  of  acetylene 
pyrolysis  above  2000  K. 

In  contrast  to  these  results,  Frank  and  Just  [43] 
found  that  "purified”  and  unpurified  acetylene 
produced  the  same  H-atom  profiles  at  50  ppm 
acetylene  and  2400  K.  Using  our  revised,  high- 
temperature  model,  H-atom  production  was  calcu¬ 
lated  assuming  both  0%  and  1.0%  (rather  than 
0. 1  %  because  purification  did  not  effect  the 
experimental  results)  acetone  in  the  acetylene.  The 
two  predicted  curves  were  indistinguishable  and 
matched  the  experimental  data.  However,  near 
2000  K  they  found  that  high-concentration  runs 
require  purification.  Again  this  result  is  consistent 
with  our  calculations.  In  fact,  our  model  predicts 
that  the  accelerated  rate  they  observe  for  the 
initiation  process  at  1800-2000  K  can  be  explained 
by  the  presence  of  small  amounts  of  acetone 
(-0.1%  in  acetylene).  The  relative  effects  of 
acetone  are  simply  explained  by  the  low  concen¬ 
trations  in  the  Frank  and  Just  experiments  as  well 
as  the  amount  of  H-atoms  produced  (relative  to  the 
concentration  of  acetylene).  The  low  concentra- 
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Fig.  8.  a.  A  comparison  of  T.O.F. /shock  tube  data  of  acetylene  decay  obtained  by  Wu,  Singh, 
and  Kem  (Fig.  2  of  ref.  8)  with  computations  using  a  revised  high-temperature  mechanism  (see 
text).  An  impurity  of  1000  ppm  acetone  in  the  3.2%  acetylene  is  assumed  for  the  calculations. 
Experimental  data:  □,  2032  K;  V  2147  K;  A  2534  K; - model. 


dons  favor  the  initiation  process,  reaction  68, 
which  is  first  order  in  acetylene.  Furthermore,  at 
2420  K,  the  experimental  results  show  a  produc¬ 
tion  of  H-atoms  (Fig.  3  in  Ref.  43)  which  is 
approximately  10%  of  the  initial  concentration  of 
acetylene.  From  1%  acetone  in  the  acetylene,  one 
could  expect  a  maximum  of  2  %  H-atoms  (one  H- 


atom  from  each  methyl  radical),  much  smaller 
than  the  measured  value.  This  situation  changes  at 
2050  K,  when  H-atoms  were  found  (experimen¬ 
tally)  to  reach  only  0.2%  of  the  initial  acetylene 
concentration.  Thus  under  these  conditions,  radi¬ 
cal  production  from  the  decomposition  of  impuri¬ 
ties  could  perturb  that  from  the  pure  pyrolysis. 


Fig.  8.  b.  A  comparison  of  experimeniaJ  (Fig.  3  of  ref,  8)  and  model  calculations  for  the 
production  of  diacetylene  during  the  pyrolysis  of  6.2%  acetylene,  at  2178  K,  and  0.48 

atmospheres.  Model  calculations  were  performed  assuming  0  (-■-),  500  ( — ),  1000  ( - ), 

and  2(X)0  (•  •  •)  ppm  acetone  in  the  acetylene. 
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ALTERNATIVE  INITIATION  STEPS 

The  focus  of  this  article  has  been  to  identify 
acetone  as  the  principal  initiator  of  radicals  during 
the  pyrolysis  of  acetylene  near  1000  K.  At  higher 
temperatures,  by-products  from  the  acetylene- 
acetone  system,  e.g.,  ketene  and  allene-methyl- 
acetylene,  contribute  to  initiation  of  radicals  (see 
previous  discussion  on  High-Temperature  Model¬ 
ing).  Efforts  to  perform  detailed  modeling  for  this 
aspect  of  the  acetylene  system  have  been  relatively 
minor  because  this  is  not  the  main  point  of  this 
article,  although  this  phenomenon  can  be  very 
important  depending  on  experimental  conditions. 
The  explanation  of  the  contribution  of  these 
species  despite  their  low  concentrations  is  related 
to  the  fact  that  the  energy  required  for  bond 
cleavage  in  these  and  similar  species  is  low  (70-90 
kcal/mol),  especially  relative  to  that  required  to 
break  the  HCC-H  bond,  i.e.,  125-130  kcal/mol. 
It  is  also  reasonable  to  expect  that  thermal  decom¬ 
position  of  other  product  species,  particularly 
olefins  such  as  vinylacetylene  and  styrene,  may 
contribute  to  initiation  once  they  attain  high 
concentrations.  However,  it  is  difficult  to  quantify 
their  role  in  acetylene  pyrolysis. 


SENSITIVITY  AND  REACTION 
PATHWAY  ANALYSES 

A  sensitivity  analysis  has  been  pierformed  using 
CHEMSEN  for  the  isobaric,  isothermal  conditions 
of  Munson  and  Anderson.  Results  of  these  calcula¬ 
tions  are  reproduced  in  Table  7  for  973  K  after 
10%  decomposition  for  the  six  reactions  most 
sensitive  to  the  computed  results  for  each  of  the 
five  species  listed  in  Table  7.  Rate  constants  for 
several  of  these  reactions  were  determined  in  this 
work.  Also,  a  reaction  pathway  analysis  was 
performed  using  programs  developed  at  UTRC. 
Reactant  decay  and  product  formation  as  deter¬ 
mined  using  this  analysis  are  described  below. 

Acetylene  was  calculatd  to  decay  principally  by 


three  reactions: 

H  +  C2H2  «  C2H3, 

(2) 

C2H3  +  C2H2  «  n-C4H5, 

(3) 

n-C4H5  +  C2H2  CftH^  +  H . 

(33) 

The  net  rates  of  the  first  two  are  essentially 
equal  whereas  the  rate  of  the  third  is  slower 
because  of  the  decay  of  n-C4H5  to  vinylacetylene. 


TABLE  7 

Normalized  Sensitivity  Coefficients  (dy/dk  x  k/yY 


CjHj  C2H4  C4rt.  C»H»  Hi 


2)  H  +  CjHi  ~  C:Hj 

3)  CjHj  +  CjHj  -  n-C.Hj 

4)  C,Hj  +  CjH,  «  C,H* 

5)  CHjCOCH,  CH,  +  CHjCO 

23)  CiH,  +  CHjCOCHj  -«  CjH4  +  CHjCO  +  CH, 

32)  n-C4Hj  «  C4H4  +  H 

33)  n-C4H,  +  C,Hj  ->  C»H«  +  H 

34)  CjH,  +  C4H4  -  C*H*  +  H 

35)  n-C4Hj  +  C4H4  C,H,  +  H 

36)  C,H,  +  H  -  C,H,  +  H, 

41)  C,H,  +  C4H4  CjH4  +  n-C4H, 

45)  CiH,  +  C*H*  «  C1H4  +  CJi, 


-0.014 

0.06 

0.13 

0.09 

-0.52 

-0.030 

0.10 

0.18 

0.21 

0.61 

0.025 

-0.42 

-Oil 

-0.23 

-0.55 

-0.050 

0.85 

0.22 

0.47 

1.11 

— 

0.20 

-0.01 

-O.OI 

-0.02 

-0.020 

0.35 

0.68 

-0.14 

0.64 

-0.035 

0.13 

-0.13 

0.73 

-0.04 

— 

-0.02 

-0.12 

0.10 

-0.01 

— 

-0.01 

-0,30 

-O.IO 

0.54 

— 

-0.08 

-0.02 

-0.04 

0.54 

— 

0.18 

— 

— 

— 

— 

0.22 

— 

— 

— 

*  As  calculated  by  CHEMSEN  at  10%  decomposition  of  acetylene,  973  K,  20%  acetylene  in  1  atm  of  argon. 
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Vinylacetylene  is  formed  almost  entirely  by  H- 
atom  elimination  from  n-C4H5  (Reaction  32).  Its 
decay  is  due  predominantly  to 

n-C4Hj-t-C4H4-CgHg  +  H.  (35) 

which  in  turn  forms  styrene,  although  H-atom 
abstraction  from  vinylacetylene  by  phenyl  or  vinyl 
(R41-42  ,  64-65)  also  contributes  to  vinylacety¬ 
lene  decay.  The  formation  of  benzene  is  from  n- 
C4H5  addition  to  acetylene  (R33).  As  the  reaction 
proceeds,  additional  benzene  is  produced  by 
phenyl  (from  R38)  abstracting  H-atoms  from 
molecular  hydrogen,  vinylacetylene  or  ethene. 
Ethene  arises  mainly  from  H-atom  attack  on 
styrene  (R38),  although  some  is  produced  by 
abstraction  of  H-atoms  from  hydrogen  and  vinyla¬ 
cetylene  by  vinyl  radicals  (R47  and  R65).  Accord¬ 
ing  to  the  proposed  model,  which  underpredicts 
hydrogen  by  a  factor  of  5  to  10,  hydrogen  is 
formed  by  H-atom  attack  on  styrene  (R36). 

DISCUSSION  OF  RATE  CONSTANTS 
AND  UNCERTAINTIES 

In  order  to  fit  the  Munson  and  Anderson  data, 
values  of  four  rate  constants  have  been  modified 
from  values  previously  proposed.  The  reactions 


are 

CjHj-t-CjHj  «  C4H6, 

(4) 

n-C4H5  C4H4  +  H, 

(32) 

n-C4H5-(-C2H2-CeH6  +  H, 

(33) 

n-C4H5  -t-  C4H4~*C8Hg  +  H, 

(35) 

kt,  the  rate  constant  for  vinyl-vinyl  recombina¬ 
tion,  was  selected  to  be  2  X  10'^  cm^  mol"' sec" ' 
in  this  work.  This  value  although  slightly  high  is 
consistent  with  rate  coefficients  for  recombination 
of  other  low-molecular-weight  hydrocarbons  (see 
ref.  44).  The  value  used  here  is  about  50%  larger 
then  the  expression  suggested  by  Ebert  et  al.  (45 J. 
The  rate  constant  for  R32  was  based  on  the 
thermochemcial  estimate  by  Weissman  and  Ben¬ 
son  [46]  but  the  A-factor  was  increased  by  60%  to 
enhance  the  fit  to  experimental  data.  Equilibrium 


calculations  give  a  reverse  rate  constant  of  about  8 
X  10'^  cm^  mol" '  sec" '  at  1000  K.  This  value  is 
about  a  factor  of  5  higher  than  that  used  for  the  H 
+  C2H2  reactions.  Errors  in  the  thermochemistry 
and  possible  contributions  from  reactions  involv¬ 
ing  vinylidene  may  be  the  cause  of  the  high  rates 
for  R32  and  R(-32)  as  required  by  this  analysis. 
At  1000  K,  the  rate  constant  for  the  irreversible 
reaction,  R33,  is  about  a  factor  of  4  [46],  2  [16, 
33],  and  1.5  [47]  lower  than  values  previously 
proposed.  The  major  difference  of  this  work  is  the 
high  activation  energy  (10  kcal/mol)  compared  to 
other  evaluations:  3.7  [47],  4.9  (from  a  fit  to 
reported  [33]  calculations),  6.9  [46],  and  9  [16]. 
Thus,  the  extrapolated  value  from  the  present 
work  is  essentially  identical  to  other  values  [33, 
47]  near  1200-1300  K.  Atjj  =  1.1  x  lO'^e-scoo^'^T 
as  used  in  this  study  is  substantially  lower  than  the 
high  value  suggested  by  Colket  [16]  but  is  about 
four  times  larger  than  the  value  reported  by  Cole 
et  al.  [47].  As  discussed  previously,  the  i-C4H5 
radical  may  play  a  role  or  C4H4-C4H4  reactions 
may  contribute. 

Rate  constants  as  adjusted  in  this  work  are 
generally  within  a  factor  of  4  of  other  values 
recently  proposed  and  experimental  data  for  these 
reactions  are  extremely  limited:  usually,  rate 
expressions  have  been  estimated  or  determined 
from  detailed  modeling.  Consequently,  we  believe 
that  the  evaluations  of  rate  constants  in  this  study 
are  very  reasonable  and  lend  strong  support  for  the 
general  thesis  of  this  article  regarding  the  impor¬ 
tance  of  acetone. 

Uncertainties  in  the  analysis  include  possible 
contributions  because  of  reactions  involving  viny¬ 
lidene  and  related  species.  As  discussed  previ¬ 
ously,  we  are  unable  to  disprove  the  existence  of 
such  mechanisms.  The  generally  good  agreement 
between  the  modeling  (without  invoking  such 
mechanisms)  and  experimental  data  supports  argu¬ 
ments  that  contributions  from  these  reactions  are 
negligible;  instead,  the  chain  mechanism  seems  to 
dominate  in  most  cases.  It  is  possible  that  the 
vinylidene  mechanism  contributes  to  vinylacety¬ 
lene  production  at  873  K  (Fig.  3)  when  the  chain 
underpredicts  the  experimental  data.  Other  uncer¬ 
tainties  are  related  to  the  lack  of  a  model  for 


C-22 


M.  B.  COLKET  ET  AL. 


methane  formation  at  low  temperatures.  At  973  K, 
Munson  and  Anderson  found  that  about  1.2%  of 
the  initial  carbon  was  converted  to  methane.  Often 
it  is  suggested  that  this  methane  formation  in¬ 
volves  heterogeneous  catalysis.  An  alternative  and 
speculative  mechanism,  described  in  this  paper  but 
not  included  in  the  model  involves  the  intermedi¬ 
ate  formation  of  the  i-C4H5  radical.  It  may  be 
formed  by  H-atom  addition  to  vinylacetylene  or 
isomerization  of  the  excited  n-C4H5  radical  after 
its  formation  by  vinyl  addition  to  acetylene.  The 
formation  of  excited  radicals  via  radical  additions 
to  acetylenes  is  another  area  of  uncertainty  in  the 
proposed  model.  A  limited  number  of  QRRK  [25] 
calculations  have  been  performed  and  indicate  that 
these  effects  should  be  included  at  elevated  tem¬ 
peratures.  An  example  is  the  system 

CzHj-hCjHz  «  n-C4H*, 

n-C4H*  «  n-C4H3, 

n-C4H*  C4H4  +  H. 

n-C4H*  «  i-C4H5, 

but  for  simplicity,  these  effects  have  not  been 
included  in  the  present  analysis. 

Significant  uncertainties  exist  in  the  kinetics  and 
mechanisms  related  to  the  odd-carbon,  acetylenic 
species.  This  belief  is  based  on  the  facts  that  very 
little  experimental  data  or  estimates  on  such 
species  exist;  large  uncertainties  exist  in  the 
thermodynaniics  of  the  odd-carbon  radicals;  and 
predictions  of  profiles  of  odd-carbon  species, 
although  qualitative,  are  not  quantitively  accurate. 
Despite  these  drawbacks,  the  major  conclusion  of 
this  article,  that  is,  that  acetylene  pyrolysis  is 
initiated  by  a  small  amount  of  acteone,  should  not 
be  altered. 

CONCLUSIONS 

We  are  unable  to  disprove  contributions  of  the 
vinylidene  mechanism  or  the  importance  of  the 
initiation  Reaction  lb  to  acetylene  decomposition. 
Recent  thermodynamic  estimates  rule  out  Reaction 
la  as  an  important  contributor  to  radical  initiation. 
Based  on  the  analysis  and  detailed  modeling 


presented  in  this  article,  we  conclude  that  the 
pyrolysis  of  pure  acetylene  has  not  been  investi¬ 
gated  at  temperatures  below  2000  K.  The  very 
high  C-H  bond  strengths  in  acetylene  prevent  a 
rapid  initiation  of  radicals;  instead  very  low 
concentrations  of  impurities  are  the  source  of 
radicals  which  initiate  and  sustain  a  chain  decom¬ 
position.  Acetone,  a  contaminant  present  in  exper¬ 
imental  studies  of  acetylene  pyrolysis,  has  been 
shown  via  chemcial  kinetic  modeling  to  dominate 
radical  initiation.  The  overall  reaction  rate  deter¬ 
mined  both  from  a  simple  steady-state  analysis  and 
dtailed  chemcial  kinetic  modeling  has  an  activation 
energy  close  to  40  kcal/mol  and  a  reaction  order 
close  to  1.5.  The  overall  activation  enery  is 
consistent  with  previous  experiments  whereas  the 
reaction  order  is  lower  than  the  usually  assumed 
value  of  2.  A  detailed  review  of  the  literature 
shows  in  fact  minimal  support  for  a  reaction  order 
of  2  while  experiments  are  reasonably  consistent 
with  an  order  of  1.5.  A  few  rate  constants  have 
been  selected  to  provide  a  fit  to  the  Munson  and 
Anderson  data.  Values  of  these  rate  constants  are 
in  general  consistent  with  values  previously  pub¬ 
lished.  The  modeling  makes  no  attempt  to  include 
reactions  involving  the  vinylidene  radical.  Al¬ 
though  there  is  no  proof  against  mechanisms 
involving  this  carbene,  there  appears  now  to  be 
ample  evidence  in  support  of  a  chain  mechanism— 
firstly,  experimental  evidence  from  a  wide  variety 
of  sources  and  conditions  support  the  presence  of  a 
chain  and,  secondly,  a  mechanism  initiated  by 
thermal  decomposition  of  acetone  adequately  de¬ 
scribes  experimental  results.  Vinylidene-related 
mechanisms  may  contribute  to  acetylene  pyrolysis 
but  we  believe  that  a  chain  mechanism  dominates 
all  acetylene  pyrolysis  above  900  K.  The  analysis 
performed  in  this  study  is  consistent  with  previous 
experiments  indicating  that  acetylene  pyrolysis  is 
much  less  sensitive  to  the  presence  of  impurities  at 
temperatures  above  2000  K,  but  only  at  very  low 
partial  pressures  of  acetylene.  At  higher  partial 
pressures  as  for  the  case  of  the  data  of  Wu  et  al., 
modeling  indicates  that  acetone  impurities  can  also 
play  a  significant  role. 
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Appendix  D 

The  Pyrolysis  of  Acetylene  and  Vinylacetylene 
In  a  Single-Pulse  Shock  Tube 
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THE  PYROLYSIS  OF  ACETYLENE  AND  VINYLACETYLENE  IN  A 
SINGLE-PULSE  SHOCK  TUBE 

MKRKnmi  1$  COLKKr.  Ill 

Ihiilfd  Tf(lnii>liif^es  l{f\i-inilt  C.niln 
Silvrt  Imuc.  Shill  Ship  3(1 
KrisI  Ihiiljiiiil.  (IT  (I6I()S 


Aci'l'li'iii'  iiiul  viiivUitflvlfiic  have  Ik'Cii  |>m<iK/c<1  in  a  sinnlc-puUc  shock  mix’  for  ihc  leiiipcra- 
luro  laiifjf  1100  to  2-100°K,  at  total  prcsstircs  ol  approximatcK  i-ighi  aliiiosphcrcs  and  (or  dwell 
limes  of  approximaielv  700  inicrosetonds.  Initial  coiueniialions  ol  the  hydrocarbon  in  argon 
ranged  from  alroiit  100  ppm  to  df?  .  (ias  samples  were  collected  and  analyzed  using  gas  chroma¬ 
tographs  foi  hsdrogen,  and  (!|  to  (^n.-hsdrtxai  Ixms.  Ihc  data  Ironi  the  pvrohsis  of  aceislene 
exhibit  substantial  production  of  \  ins  lacctylene.  Isenzene.  and  phenylacetylene.  but  agree  scell 
sciih  a  detailed  chemical  kinetic  motlel.  Data  IVom  vinvlacelslene  psrolssis  and  ihermochemical 
arguments  suggest  a  chain  mechanism  bs  svhich  11  adds  to  vinslacelslene  and  the  resultant 
adduct  decottiposcs  to  acetvletie  and  a  sinyl  radical.  Rate  constants  lor  the  reserse  steps  ot  those 
occurring  during  s  ins'lacetylene  and  benzene  psrolssis  have  been  calculated  usitig  ihertnody- 
tiamics  and  forsvard  rate  constants.  These  reverse  rate  constants  assist  iti  describing  the  proiluc- 
tion  ol  smslaceivlcne.  benzetie,  atid  phetivlaceislene  tiurmg  acetslene  |)srolssis. 

Font  separate  reaction  mechatiisms  for  the  itiilial  lormation  of  arotnatic  rings  have  been 
identihed.  The  relatise  impoiiaiuc  of  each  step  depends  on  the  ambient  tetnpet attire  and 
relatise  concentrations.  Ox'fitill  steps  can  be  svriilen  as 

(...II,  -i-  n-(;,lb.—  (■.,,11,.  +  II 
(  ,11,  -i-  n■(;,H.  —  (i,.ll-, 

(.,11,  -*■  ii-(:,ii,  — (:,.ii-,(;,ii,  ii 

(.,11,  -K  (;,ii,  —  (;„ii„  -t-  II 


Introduction 

Glassniati'  has  cons indiigis  argued  that  the 
cheniistrs  of  liiel  psrolssis  plass  an  iinportant 
lole  in  sotMtng  dtf  fusion  (lames.  Some  ps  rolysis 
steps  mas  he  t .tie-limiting  for  tlie  production  of 
iiKi|)ient  soot  particles,  and  therefoie  the  total 
.miount  ol  soot.'  .Xi  ined  svitlt  tins  infoi  ination. 
mans  lesearch  pitigrams  tiasc  hegtin  to  iden- 
tiis  tile  late-limitmg  pitzeesses.  for  example, 
liittnei  and  Hosvard.*  and  Bockhoiii.  et  al.' 
Iittve  prosided  experimental  confirmation  of  a 
sariets  ol  high  molecular  sveight  polstsclic 
.iromatics  present  in  licli  (lames,  (.’ole.  et  af.  ’ 
.md  Weissman  and  Benson''  h.ive  suggested 
mechanisms  sviih  rates  consistent  sviih  tlier  nio- 
chemical  analssis  .md  experimental  d.n.i  lor  the 
production  <if  single-ring  aromatics,  rrenkladi. 
et  al.'  has  modeled  soot  formation  liom  acety¬ 
lene  pvrolvsis  in  shock  tubes  using  a  detailed 
chemical  mechanism.  Ihc  mechanism  includes 
|>i<>duciion  of  aiiimatic  rings  and  contmuaf 
grosvth  of  polscslic  aromatics,  (jii.ilitat is e  suc¬ 


cess  has  been  obtained  in  this  substantial  cl  fort 
despite  the  lack  of  confimiing  experimental 
data  describing  profiles  of  intermediates  lot 
confirmation/supijort  of  tlie  proposed  model. 

It  is  the  objective  of  the  present  svork  to 
extend  tliese  earlier  ssorks  m  order  to  help 
elucidate  psrolssis  steps  and  chemical  mech.i- 
nisms  related  to  the  formation  and  bieak-iip  ol 
aiomatic  rings.  A  single-inilse  shock  tube''  '  sv.is 
selected  (or  this  stiidv,  not  onis  because  shock 
tubes  are  one  ol  the  less  devices  capable  ol 
geneialing  conditions  of  impoitance  to  soot 
foiination  in  diffusion  fl. tines,  hut  also  to 
extend  the  inloi  in.ition  on  soot  ptciduciion 
.ilreads  generated  through  optical  studies  in 
shoe  k  lubes. ' ' 

i’reliminat  s  sersions  ol  the  modeling  sstirk 
have  been  desciibed  sepat.itels  lot  acetslene" 
•md  vinvlacetslenc."  The  pieseni  svork  com¬ 
bines  the  earlier  models,  incitidmg  a  recent  one 
lor  benzene  pyrolssis,"'  into  one  model  svhich 
salisfac  lortlv  desc  ribes  both  deccimpositioii  (ol 
the  parent  hsclioc.irbons)  ,iiid  pioducl  (oi  ma- 
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lion.  Required  revisions  include  slight  adjusi- 
nients  to  rale  constanis  and  lo  assumed  ihernio- 
dynamic  parameiers. 

Description  of  Fa  dlities 

The  SI’ST  used  in  this  program  is  2E5  cm 
long  and  has  a  diameter  of  3.8  cm  (i.d.).  The 
driver  is  88  cm  in  length  and  can  be  turned  by 
shortening  its  length  in  3.8  cm  increments;  the 
driven  section  is  197  cm  long.  An  11.7  liter 
“dump  lank"  is  located  in  the  driver  (lower 
pressure)  section  30  cm  downstream  of  the 
diaphragm.  Pressure  profiles  were  determined 
using  Kistler  pressure  transducers  located 
15.25  and  2.50  cm  from  the  end  wall.  Arrival 
times  were  measured  lo  within  one  microsec¬ 
ond  using  digitized  pressure  traces.  Calculated 
quench  rales  are  typically  10^  K/sec  or  higher  in 
the  rarefaction  wave.  Starting  pressures  prior 
to  filling  are  0.2p.  and  leak  rates  are  less  than 
IpVmin.  Post-shock  temperatures  were  calcu¬ 
lated  based  on  the  measured  incident  shock 
velocity  and  normal  shock  wave  equations. 

The  procedures  for  performing  an  experi¬ 
ment  are  similar  to  those  described  by  Tsang.'^ 
except  for  an  automated  sampling  system.  The 
sample  is  collected  at  the  endwall  of  the  shock 
tube  using  0.0-15  inch  i  d.  tubing  heated  to  over 
SS^C.  Approximately  30  milliseconds  after  the 
gas  has  been  shock  heated  and  cooled,  a 
solenoid  valve  opens  to  the  evacuated  sample 
cell  and  then  closes  after  300  milliseconds.  The 
sample  storage  vessel  is  all  stainless  steel  with  an 
internal  volume  of  25  cc. 

The  sampling  volume  is  directly  coupled  to  a 
low  volume  (<3cc).  heated  inlet  system  of  a 
Hewlett  Packard  5880  A  gas  chromatograph. 
Valves,  detectors  and  software  integration  rou¬ 
tines  as  described  previously’^  enable  this  system 
to  provide  automatic  quantitative  detection  of 
hycirogen  and  hydrocarbon  species  up  to 
Cio-h\  drocartxjns.  Based  on  repealed  injections 
of  calibrated  samples,  overall  accuracies  arc 
estimated  to  be  three  percent.  Calibration  gases 
were  stored  in  stainless  steel  cylinders  with 
degreased  valves  and  were  heated  to  approxi¬ 
mately  60°C  prior  to  injection. 

Argon  (99.999%  pure)  was  obtained  from 
Maihcson  and  was  the  principal  diluent.  Com¬ 
pressed  acetylene,  also  from  Matheson,  con¬ 
tained  about  I  to  2%  acetone  (added  for 
stability)  depending  on  bottle  conditions.  Acety¬ 
lene  was  purified  by  repeated  freezing  and 
thawing  at  fiquid  temperatures  and  retain¬ 
ing  only  the  middle  50%.  Final  samples  still 
contained  0.1  to  0.2%  acetone.  Vinylacciylcnc 
was  obtained  from  Wiley  Organics  and  con¬ 


tained  an  imidcntificd  hydrocarbon  with  a 
concentration  of  afiproximaiely  8000  ppm. 

Model  Descri()lion 

Detailed  chemical  kinetic  calculations  have 
been  performed  usingCHF-MKIN,’”  lijODE,' ' 
and  a  version  of  a  shock  tube  code  originally 
developed  by  Mitchell  and  Kee'”’  but  modified 
to  include  quenching  effects  in  an  SPSl. 
Quenching  rates  are  determined  using  experi¬ 
mental  pressure  traces,  assuming  adiabatic  ex¬ 
pansion,  and  using  the  equation 

dr  T-i-\  dP 
dt  P  y  dt 

Calculated  quenching  rates  vary  as  a  function 
of  shock  strength  and  time.  Initial  quenching 
rates  are  as  high  as  2xl0®Kysec  for  shocks 
producing  initial  post-shock  temperatures  of 
2000°K,  but  only  25%  of  that  rate  for  shocks 
producing  reflected  shock  temperatures  near 
1200°K. 

The  chemical  kinetic  model  used  in  this  work 
is  reproduced  in  Table  I  and  is  based  on 
proposed  mechanisms  for  the  pyrolysis  of 
acetylene''  and  eihenc.'*  The  reaction  set  in¬ 
cludes  the  identification  of  the  C^Hi  isomers 
following  Frenklach,  et  al.'  The  radicals  arc 
denoted  n-C^Hj  for  the  (normal)  isomer  with 
the  unpaired  electron  on  the  terminal  vinylic 
carbon  and  i-CjHs  for  the  (iso)isomcr  with  the 
radical  site  on  the  interior  carbon  atom.  n-C^Hs 
(radical  on  end  carbon)  was  the  only  CjHj 
isomer  considered  in  this  work.  Also,  the 
mechanism  includes  the  forward  and  reverse 
processes  of  the  predominant  path  for  decom¬ 
position  of  benzene,,  i.e., 

QH.  -t-  H  -h  H,  (39) 

CfiHj  f-QH,  (5b) 

(-Ci.Hs  n-C^H,  -t-  (55) 

where  f-Cr.Hj  has  been  suggested"’  to  be  the 
l-hcxync-3,5-dien-()-yl  radical.  The  rate  con¬ 
stants  (R55,R5b)  have  been  resised  slightly 

from  the  previous  woik"’  in  order  to  be 

consistent  with  the  reverse  pitnesses. 

As  required,  thermodynamic  estimates  of 
some  species  were  made  using  group  addiovity 
techniques.-’’  Otherwise,  data  from  rciidih 
available  sources  were  used."'-'^  Both  estimates 
of  rate  parameters  and  those  calculated  using 
equilibrium  constants  are  dependent  on  the 
selected  thermodynamics.  Of  pariicuhir  con¬ 
cern  during  the  present  research  is  the  ;tj>par- 
eni  uncertaints  in  the  heat  of  formation  o( 
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l  AIU.E  I 

Proposed  Set  o(  Reactions  anti  Rate  (atelficients 
logk  =  logA  +  nlog  r  —  EyRy'r/2.!R)H  • 


Eoiward  Reverse 

Reactions  Rate  Constant  Rate  Constant 


logA 

t1 

E 

IngA 

11 

E 

Ref 

I 

C2H4  +  M=C2H3  +  H  +  M 

16.16 

0.0 

81.8 

14.24 

0.0 

-23.0 

22 

2 

C2H4  +  H  =C2H3  +  H2 

14.84 

0.0 

14.5 

13.53 

0.0 

14.0 

21 

3 

H  +  C2H2=C2H3 

12.74 

0.0 

2.5 

13.01 

0.0 

43.8 

21 

4 

C2H3  +  H  =  H2  +  C2H2 

13.00 

0.0 

0.0 

13.34 

0.0 

63.2 

21 

5 

H2  +  M=2H  +  M 

12.35 

“  .5 

92.5 

1  1.74 

—  .5 

-1  1.9 

21 

6 

C2H2  +  M=C2H  +  H  +  M 

16.62 

0.0 

107.0 

1 5.25 

0.0 

-  17.9 

21 

7 

2C2H2  =  n-C4H3  +  H 

12.30 

0.0 

45.9 

1 1.67 

0.0 

-25.0 

21 

8 

i-C4H3  +  H2=C2H2  +  C2H3 

10.70 

0.0 

20.0 

1  1.41 

0.0 

17.1 

P\V» 

9 

C4H4=i-C4H3  +  H 

15.20 

0.0 

85.0 

12.72 

0.0 

-12.1 

P\V 

10 

C4H4  =  n-C4H3  +  H 

15.00 

0.0 

100.0 

12.93 

0.0 

-7.7 

P\V 

1 1 

C2H  +  C4H4=C2H2  +  i-C4H3 

13.60 

0.0 

0.0 

12.48 

0.0 

27.9 

21 

12 

n-C4H3=C2H2  +  C2H 

14.30 

0.0 

57.0 

1 3.56 

0.0 

3.0 

P\V 

13 

i-C4H3  =  C4H2  +  H 

12.00 

0.0 

49.0 

12.86 

0.0 

-0.2 

P\V 

14 

n-C4H3=C4H2  +  H 

12.60 

0.0 

40.0 

13.04 

0.0 

1.4 

P\V 

13 

n-C4H3  +  H  =  i-C4H3  +  H 

13.48 

0.0 

0.0 

13.06 

0.0 

10.7 

PW 

16 

i-C4H3  +  H  =  C4H2  +  H2 

13.00 

0.0 

0.0 

14.47 

0.0 

55.2 

PW 

17 

n-C4H3  +  H=C4H2  +  H2 

13.00 

0.0 

0.0 

14.05 

0.0 

65.9 

PW 

18 

C4H4  +  H  =  i-C4H3  +  H2 

14.49 

0.0 

14.5 

12.62 

0.0 

21.9 

26. PW 

19 

C4H4  +  H  =  n-C4H3  +  H2 

13.90 

0.0 

14.5 

12.45 

0.0 

1  1.2 

20.  PW 

20 

C2H3  +  C2H2  =  n-C4H5 

12.04 

0.0 

4.0 

13.79 

0.0 

37.7 

PW 

21 

C2H  +  H2  =  H  +  C2H2 

12.85 

0.0 

0.0 

13.60 

0.0 

20.5 

22 

22 

C2H  +  C2H2=r4H2  +  H 

13.60 

0.0 

0.0 

14.78 

0.0 

15.4 

21 

23 

C4H2=C4H  +  H 

14.89 

0.0 

120.0 

13.28 

0.0 

2.1 

22 

24 

C2H  +  C4H2=C6H2  +  H 

13.60 

0.0 

0.0 

14.97 

0.0 

15.1 

21 

25 

C4H  +  C2H2=C6H2  +  H 

13.30 

0.0 

0.0 

14.91 

0.0 

8.1 

2I.PW 

26 

C6H2=C6H  +  H 

14.89 

0.0 

120.0 

13.05 

0.0 

5.3 

22 

27 

C4H  +  H2=H  +  C4H2 

13.30 

0.0 

0.0 

14.30 

0.0 

13.5 

22 

28 

C6H  +  H2=H  C6H2 

13.30 

0.0 

0.0 

14.53 

0.0 

10.3 

22 

29 

C2H  +  H=C2  +  H2 

12.00 

0.0 

23.0 

1 1 .55 

0.0 

3.8 

PW 

30 

C2H  +  M  =  C2  +  H  +  M 

16.67 

0.0 

124.0 

15.61 

0.0 

0.4 

36 

31 

C4H4  +  C0H5  =  C6H6  +  n-C4H3 

12.00 

0.0 

0.0 

12.56 

0.0 

3.1 

PW 

32 

C4H4  +  C6H5=C6H6  +  i-C4H3 

12.00 

0.0 

0.0 

12.14 

0.0 

13.8 

PW 

33 

C2H  +  C6H6=C6H5  +  C2H2 

13.30 

0.0 

0.0 

12.05 

0.0 

14.1 

34 

34 

C4H  +  C6H6=C6H5  +  C4H2 

13.30 

0.0 

0.0 

12.29 

0.0 

7.0 

34 

35 

C2H3  +  C4H2=C4H4  +  C2H 

13.48 

0.0 

23.0 

13.46 

0.0 

3.1 

PW* 

36 

C2H3  C4H4  =  C2M4  +  n-C4H3 

1  1.70 

0  0 

16.3 

1  1  .55 

0.0 

13.4 

35. PW 

37 

C2H3  +  C4H4=C2H4  +  i-C4H3 

1  1.70 

0.0 

16.3 

11.13 

0.0 

24.1 

35.PW 

38 

C6H6=C6H5  +  H 

16.18 

0.0 

107.9 

1 3.55 

0.0 

-3.0 

16 

39 

C6H6  +  H=C6H5  +  H2 

14.40 

0.0 

16.0 

12.39 

0.0 

9.5 

34 

40 

2Cf)H5  =  C12H10 

12.48 

0.0 

0.0 

16.57 

0.0 

108.3 

17 

41 

C6H5  +  C6H6=C12H10  +  H 

1  1.80 

0.0 

1 1.0 

13.27 

0.0 

8.4 

37 

42 

n-C4H3  +  C6H5-C10H8 

13.00 

0.0 

0.0 

0.00 

0.0 

0.0 

PW 

43 

C2H2  +  C6H5=C8M6  +  H 

12.00 

0.0 

4  (» 

13.50 

0.0 

3.1 

38.  PW 

44 

C2H4  +  C0H5  =  C8H8  +  M 

1  1 .57 

0.0 

2.1 

12.99 

0.0 

1.6 

38 

45 

C2H3  +  C4H4-C6H6  +  H 

1  1 .60 

0.0 

0.0 

0.00 

0.0 

0.0 

PW 

46 

n-C4H5  =  C4H4  +  H 

14.00 

0.0 

414 

13.42 

0.0 

3.2 

6 

47 

n-C4ll5  +  C2H2  =  /-C6H7 

12.81 

0.0 

9.0 

14.80 

0.0 

44.4 

6,PW 

48 

n-C4H5  +  H  =  C4H4  +  H2 

13.00 

0.0 

0.0 

13.03 

0.0 

66.2 

PW 

49 

C4HG  +  H  =  n-C4H5  +  H2 

14.00 

0  0 

14  5 

12.80 

0.0 

1  1.2 

26.  PW 

50 

C6H6  +  C2H=C8H6  +  H 

12.00 

0.0 

0.0 

12.25 

0.0 

13.2 

PW 

51 

C4H4  =  C2H  +  C2M3 

16.00 

0.0 

105.0 

13.46 

0  (1 

-  15. 5 

PW 

RKACTION  KINETICS 


I'AIU.E  I  (Cuu(iiiurd) 


E'orward 

Reverse 

Reactions 

Rate  flonslant 

Rale  Uoiistaiu 

logA 

n 

E 

logA 

n 

F, 

Ref 

52 

C4fH  -r  n-C4U5-C8U8  -r  M 

1 3.90 

0.0 

3.0 

0.00 

0.0 

0.0 

I’W 

5.3 

C8H8  +  H-C8M6  -t-  H  -r  H2 

14.60 

0,0 

7,0 

0.00 

0.0 

0.0 

1>\\ 

5-1 

C6U5  -I-  C2H3  =  C8M8 

1 3.00 

0.0 

0.0 

16.34 

0.0 

104.4 

1>\V 

55 

/-C6H5=n-C4H3  -I-  C;2H2 

14.00 

0.0 

36.0 

1  1.71 

0.0 

-0.1 

I’W 

56 

C6H5=/-C6H5 

13.54 

0.0 

65.0 

10.22 

0.0 

1.4 

I’W 

57 

C6H6  -t-  H=c-C6H7 

13.60 

0.0 

4.3 

13.12 

0.0 

24.6 

I’W 

58 

c-C6H7  =  /-C6H7 

14.48 

0.0 

50.0 

1 1.36 

0.0 

0.4 

I’W 

‘NOTES:  Units  for  A:  cc.nioles.sec. 

Units  for  E;  kcal/mole. 

"  =  "  represents  forward  and  reverse  dircttioiis  included  in  model. 

represents  forward  direction  only  included  in  motlel 
PW  indicates  rate  evaluted  from  the  present  work. 


vinyl-'’-'  and  related  (e.g.  n-QHs  and  n-QHs) 
radicals.  Uncertainties  in  the  heats  of  formation 
and  entropies  of  these  radicals  directly  translate 
into  errors  in  proposed  rate  constants.  Ther¬ 
modynamic  parameters  at  room  temperature 
for  species  included  in  this  work  are  presented 
in  Table  11. 

RcMilts  and  DLaciL'isioit 

Several  scries  of  runs  of  the  pyrolysis  of 
acetylene  and  vinylacetylenc  have  been  com¬ 
pleted.  Each  series  represents  approximately  10 
to  15  separate  experiments  and  each  experi¬ 
ment  in  a  series  has  the  same  initial  concentra¬ 
tion  in  argon.  For  each  run,  a  chemical  analysis 
is  performed  of  final  product  distribution. 
Initial  post-shock  temperatures  for  runs  from  a 
given  series  typically  range  from  1100  to 
‘24(K)'’K.  Detailed  chemical  kinetic  model  calcu¬ 
lations  have  been  performed  for  most  series 
and  the  numerical  results  for  selected  .series  are 
compared  to  the  experimental  data  in  Figs, 
1,3—5.  Mass  balance  data  were  |)resemed 
previously.'^  A  detailed  discussion  of  the  ex¬ 
perimental  data  and  kinetic  mechanism  follows. 

Acetylene 

The  dcc(jm[jo.siiion  of  3.79?  aceiyleiic  in 
argon  ami  pioduction  of  the  major  species,  i.e. 

H'j,  and  Ui.H^,  agree  well  with  existing 
kinetic  models.  As  seen  in  Fig.  1,  vinslacets- 
Icne  (VA).  ben/ene  and  ])henylacetylenc  are 
also  obseived  anti  may  play  a  critical  role  in 
the  growth  and  |)roduciion  of  polycyclic  aro¬ 
matics.  The  bimolecular  reaction  — ► 

(bH|  is  ihouglit  to  occur  at  low  temperatures. 


Extrapolation  of  low  temperature  rate  con¬ 
stants  provided  quantitative  agreement  with 
the  experimental  piofile.  In  addition,  the 
activation  energy  for  this  process  (~-I0  kcal/ 
mole)  is  in  reasonable  agreement  with  the 
assumption  of  intermediate  formation  of  the 
C4H4  diradical,  which  is  approxitnatclv  3.5  to 


TABLE  II 

Selected  Tliermodv  naiiiics  at  3t)l)K 


Sf/ecirs 

Unit  tij  Fiirmatnin 

(kcal'iiiole) 

t'.Utlllp\ 

1  111 

43.6 

93.7 

in 

36.0 

79.5 

ChM., 

78.3 

79.6 

97.4 

81.8 

c-(;„ll: 

49.9 

72.1 

19,9 

65.2 

c;.,H, 

/  H . ') 

69  4 

f-c;.,ii.. 

1 39.2 

so.n 

(;.ii. 

169.7 

71.1 

(;.,n 

233  2 

74  3 

C.ll., 

26. 1 

()(),(') 

n-C.H, 

82.5 

69.1 

C.ll. 

69.4 

66  1 

125.1 

68  7 

i-C.II. 

1  15.2 

71,6 

C.II, 

1  1  1.7 

59,9 

C,ll 

179.0 

62  8 

Cdl, 

12  (i 

52,4 

(;..ib 

.') 

1.,, 

54.2 

48  1 

(;,ii 

128.5 

49  (i 

(b 

200.2 

47  7 
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-1!^  kcal  more  encrj;cii(  (han  iwo  moleeiiles  ol 
acetylene.''^  Unlbriiinaiely.  iliis  medtaiiism 
does  not  provide  a  plausil>le.  parallel  |)aili  for 
the  formation  of  benzene;  \el  the  similarity  in 
profiles  of  l)en/.ene  and  \'A  stroni'ly  Mij'i^esl  a 
paiallel  meehanisiti.  Stein"'  sni^i^esied  an  alter¬ 
native  chain  process  itiitiated  l)\  H-aiotn  atldi- 
tion  to  acetylene  to  fprin  VA,  i.e., 

C.Hjf  +  H)—  C..H,(  -t-  CwH..) 
n-(;,Hr,(-H)  —  Ci.H, 

'I'liis  setitience  is  identical  to  the  fust  steps  in 
the  t  ing  fbnnation  process  proposed  h\  Frenk- 
lach,  ct  al.'  The  final  steps  of  the  proposed 
process  include 

— ►  phenyl 

followed  by  subsequent  formation  of  bettzene 
or  phenylacetylene. 

Detailed  ciiemical  kinetic  inodelitig  of  the 
present  experimental  data,  however,  stionglv 
suggests  that  early  (i.e.  low  lempertmtre)  ben¬ 
zene  fortnatiott  aiises  principally  from  atetv- 
lene  addition  to  the  tioi  inal-butadien\ I  radical, 
not  tlie  n-C4H:i  radical,  i.e.. 

n-CbH-,(  -t-  -»  — 

c-(;„Ht(-H)- 

where  f-CiiH;  is  delitied  to  be  the  l.:k.^-he\a- 
trien-l-yl  radical. 

I’relitninary  ctilculatiotis  using  this  mecha- 
tiism  have  also  beett  perfortned  lot  comparison 
to  the  flow  reactor  data  on  aceivlene  pvroivsis 
oljtained  by  Mittison  a.  d  .-Xnderson.^''  Seini- 
(|uaniitative  agreetnent  lor  the  production  of 
betizetie  atid  \'A  was  achic\e<l  and  thus  the 
bitnoleciilar  reaction  involv  ing  inlet  inedi. tie 
formation  ol  a  (!|H|  diradical  is  not  nevess.irv  to 
descrilve  lower  temperature  aceivlene  |)vrolvsis. 

.Above  ir»<)()‘’K.  cyclic  compounds  are  formetl 
principallv  by  aceivlene  addition  to  n-(;,U(  in 
agreement  with  the  proposal  bv  Frenkkich.  el 
al.'  Previously  it  was  assumed,  however,  that 
n-(itH:t  vvas  formed  through  path  .A.  wheieas 
(R-l'i)  was  the  dominant  louie  undei  the 
pre.sent  conditions 

n-C,Hd-H)  -»  (hH,(-H)  -*  n-t;,Hi  A 
C.H,.  -r  (;,.H  —  n-(;,H,  (-12) 

This  conclusion  arises  from  results  dis- 
(itssed  later  in  this  papei  which  show  that  V.A 
decomposition  piincipally  involves  H-alom  ad¬ 
dition  to  CiH^.  not  H-aiom  absiiaciion.  .As 
leinperttture  inc/eases,  the  .ibsirac lion  route 


will  become  more  signifK.mi;  howevei.  the 
toiuenlralion  ol  (!_.ll  will  also  mcie.ise  so 
(K-12)  still  remains  competitive,  i- in  i  bet  nioi  e 
(R-I’J)  is  ;i  SI  I  .light  lot  vv  at  cl  pi  oc  ess  lor  pio- 
duc  ing  the  not  inal  l  adic  al  and.  in  f  ig.  I. 
assists  in  .idecjiialely  desc  iibing  the  piodnclion 
of  benzene  and  ))henv  lac  elv  lene.  Demonsii  .i- 
lion  ol  the  relative  sources  of  n-(!|lli  tail  be 
seem  in  fig.  2.  ivhere  icMction  souiccs  .md 
sinks  ate  plotted  as  a  lunciion  of  time  lor  two 
dillei  ent  lent  pet  al  tires. 

figure  I  shows  that  the  model  oven  predicts 
observed  benzene  and  pbenviac  el  viene  c  one  en- 
tr.ilions  tibovo  l()()()°R.  ( ioiiu  idenitillv .  this  is 
the  s.ime  tenipeiatuie  at  which  a  significant 
deficit  in  recoveied  mass  first  appeals.''  .Some 
phenylacet viene  and  benzene  above  |(')Ol)°K  is 
probablv  converted  to  bigber  molec  ular  weight 
species  which  .ire  not  observed  and  not  .ic- 
connied  lot  bv  the  model. 

Viinltirii\lciir 

Nfixtiires  ol  O.t)  1 . 0. 1  I .-)  and  I  A)'/i  \’.A  in  argon 
weie  pviolvzed  ovei  the  lemperaiure  r.mge 
I  lot)  to  2.’)(>I)‘’K.  File  data  lor  the  lowest  and 
highest  concentrations  are  |)i  esenied  in  f  igs.  A. 
4  ;md  .n.  Fhe  oihei  species  not  shovvii  in  these 
figures  vvere  observed  .it  conceniraiioiis  less 
ih.m  2'f  of  the  p;ireiii.  I'hese  species  include 
met h. me.  eihv  lene.  allene/meihv  lac  eiv  lene.  sev  - 
eral  unidentified  (i-,  spec  ies.  ( ;,,Ff toluene,  and 
.1  (i.i-hvclrocai  bon  (possiblv  inclene).  Dai.i  on  the 
seiies  with  interniecli.iie  cimcenir.iiion  h;ive 
fjeen  piesenieci  prev  iouslv . '  '  Fhe  most  signili- 
c;ml  infot  nialioii  is  that  al  low  lemperal iires  and 


INITIAL  post-shcx;k  temperature  (K) 

f  ic.  I .  f.\|>ci  iiiic'iii.il  .111(1  .Model  Kesulis  loi  I’v  i  olv  - 
SIS  ol  .t.7V  .Ac  elv  lene'  in  .\igoii.  DvviTI  I  line  — 
7l)l)x  10  ''  sec.  I  ol.il  |nessnie  =  M  .iini 
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REACTION  KINETICS 


Fic.  2.  Nei  Formaiion  Rate  for  Produciion/Dc- 
struciion  of  the  n-C<Hj  Radical.  2a)  at  I600K.  2b)  at 
1800K. 

Contributing  reactions: 

R7  C2H2  +  C2H...sin-C<H,  +  H 
R12  n-c/Hj^iCjH;  +  C;H 
R14  n.C<Hj;:tC4H2  +  H 
RI9C^H,  +  Hs2n-C4H,  +  ff.. 

R55  l-CeHs^in-C^H,  +  C2H2 

Quenching  wave  arrives  at  700 x  10“*’  sec.  Total 
pressure  =•  8  aim.  Pyrolysis  of  3.7Sf  acetylene. 

low  initial  concentrations,  acetylene  is  the  pre¬ 
dominant  product,  while  diacetylene  (as  well  as 
H2)  is  produced  at  levels  approximately  ten 
times  less  than  that  of  acetylene.  Aromatics  arc 
produced  at  higher  temperatures  and  concen¬ 
trations,  but  not  in  quantities  suiricieni  to 
encourage  rapid  growth  of  polycyclics  and 
severe  mass  imbalance.  Relatively  little  data  have 
been  presented  previously  for  comparison  to 
this  SPST  data.  Lundgard’s^'  review  has  identi¬ 
fied  one  previous  work  which  overlaps  the 
present  tcm|)erature  range.  Yampol'skii,  ct  al.’^ 
examined  VA  pyroivsis  between  I02;f  and 


Fic.  3.  Experimental  and  model  results  for  100 
ppm  Vinylacetylene  in  Argon.  Dwell  time  = 
700xl0"‘  total  pressure  =■  8  atm. 


Fic.  4.  Experimental  and  .Model  Results  for  Ali¬ 
phatic  Species  during  Pyrolysis  of  17c  Vinylacetylene 
in  Argon.  Dwell  time  =  700xl0'®sec,  total  pressure 
=■  8  atm. 


INITIAL  POST-SHOCK  TEMPERATURE  (K) 


Fig.  Experimental  and  Model  Results  for  Aro¬ 
matic  Products  during  Pyrolysis  of  1%  Vinylacetylene 
in  Argon.  Dwell  lime  =■  700xl0"®sec,  total  pressure 
=»  8  atm. 
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1273°K  and  produced  qualitatively  similar  te- 
sults  to  the  SPST  work.  The  main  low  tempera¬ 
ture  aliphatic  found  b\-  Yampol’skii  was  acety¬ 
lene  with  traces  of  methane  and  ethene,  al¬ 
though  diacetylene  was  not  detected.  Benzene 
was  also  observed  previously,  as  well  as  a 
polymer  whose  precursor  may  be  styrene.  The 
ov  erall  decomposition  of  \'A  observed  by  \  ani- 
pol’skii.  et  al.,/r=  1.6x  lO''  exp  (  — 52,800cal/RT) 
sec” '  is  approximately  five  times  higher  than  the 
SPST  results. 

Kinetic  modeling  requires  the  identification 
of  the  initiation  step  and  a  chain  process  which 
describes  the  predominant  formation  of  acetv- 
lene  vvith  minor  production  of  diacetylene, 
styrene,  benzene  and  phenylacetylene.  The 
three  possible  initiation  processes  are 

R9  C^H4  ^  i-C^H3  H  -  98  kcal/mole 

RIO  C^H^  ->•  n-QHj  -t-  H  -  1 1 1  kcal/mole 

R51  C^H^  — CoH  CoHj  -122  kcal/mole 

where  the  estimated  endothermicities  mav  each 
be  in  error  by  as  much  as  10  kcal/mole  due  to 
uncertainties  in  heats  of  formation  of  the 
hydrocarbon  radicals.  In  this  work,  the  princi¬ 
pal  initiation  process  was  assumed  to  be  Reac¬ 
tion  9.  which  required  a  rate  expression  of 
10'-’ *  exp  {-42800/T)  sec”',  although  the  .-\- 
factor  seems  high  and  E.^,.  low,  for  the  C-H 
bond  scission.  It  is  important  to  note  that  the 
model  results  were  relatively  insensitive  to  the 
absolute  magnitude  of  the  initiation  rate  or  its 
temperature  dependence.  A  factor  of  three 
change  in  the  initiation  rate  lesulted  in  ap¬ 
proximately  a  20-309J  change  in  the  overall 
rate  of  decomposition  and  the  formation  of 
products.  The  explanation  for  this  phenome¬ 
non  is  that  the  decomposition  of  \'.A  in  the 
range  1200-1400°K  is  controlled  bv  a  chain 
mechanism; 


R9 

t:... 

R-*(i 

(mHiCI  ChU,  -  H  13.2 

8,3000 

ih.tin 

R-'_>() 

H  *  CjHiSn-C.PL,  1.3  (1 

1370 

chjtn 

R-3 

ii-C4H-,ciC..H,  +  C.H,.  H,(l 

33.3110. 

c  luin 

CjH.srsH  -  car.  13.0 

33710 

cc.iiidIc.scc 

c.il  riioli* 

where  the  last  three  reactions  account 

for  the 

majoriiv 

of  the  decomposition  of  the 

jjareni 

and  the 

predominance  of  acetylene 

in  the 

low-temperature  products.  This  chain  is  not 
onlv  important  specifically  to  understanding 
the  pv  roivsis  of  v  inylacetylene;  but.  as  discussetl 


earlier  in  this  paper,  also  to  the  mechanism  ol 
ring  formation  during  ticetylene  pyrolysis.  Ex¬ 
planations  for  eliminatitvn  of  other  mechanisms 
are  presented  below. 

Alternative  chain  mechanisms  for  decompo¬ 
sition  of  would  include  H-atom  abstraction 
from  C^H4  resulting  in  C:CCH:CH.., 
HCrCC.CHj,  or  HCiCCH.CH.  Formation  of 
the  first  of  these  can  be  neglected  from  thermo¬ 
dynamic  considerations.  In  addition,  direct 
formation  of  two  acetylene  molecules  from  this 
structure  seems  unlikely.  The  second  isomer. 
i-C4H;i.  should  decompose  principally  into  dia¬ 
cetylene,  since  breakage  of  the  C-C  bond  would 
have  a  high  activation  barrier  due  to  the 
formation  of  vinvlidene  as  an  intermediate. 
The  last  isomer.  n-CjHs,  may  decompose  into 
C4H2  +  H  or  C^H  C^Ho.  Estimated  rates 


1 

Log  k 

E 

RI2 

n-C4H:,i=C..H  -  C..H , 

3,72 

14.3 

37000. 

RI4 

n-C4H,=:(;4H..  -  H 

5.88 

12.6 

40000, 

suggest  that  Reaction  14,  i.e.  production  of 
C4H2,  is  faster  by  a  factor  of  14  near  1300°K 
where  \’.A  decomposition  is  observed,  .At 
higher  temperatures  the  twf)  lates  approacli 
one  another.  These  relative  rate  estimates  are 
supported  bv  experimental  decomposition  tlata 
of  benzene  pv  lolysis  from  Kern,  et  al.'*'  Those 
researchers  fouitd  that  the  initial  ratio  of 
products  C2H2/C4H2  at  1700°K.  is  approxi- 
niately  two  to  one.  .Assuming  the  overall  decom¬ 
position  path 

C.iH„(-H);=:C„H3^n-C4H  ,  -r  C.H2. 

then  the  C2H2/C4H2  product  ratio  suggests  iliai 
the  branching  ratio  A14/A12  is  also  approximatelv 
two  to  one.  The  calculated  ratio  using  the  above 
rates  is  three  (3)  at  1700°K.  Thus,  if  n-C4H4  is 
the  principal  intermediate  from  \'.A  decompo¬ 
sition,  then  the  initial  production  rate  of  C4H2 
would  be  similar  to  or  higher  than  that  ol 
during  \’.A  pvroivsis.  Such  a  prediction  is  not 
substantiated  bv  the  SFST  experiments. 

H-atoms  mav  add  to  \’.A  at  locations  other 
than  to  the  secondary  acetvienic  carbon.  .Addi¬ 
tion  to  the  other  carbons  would  form  H  >C: 
CCHCH.,  HCfCCHCH:,.  or  HC. CCH.CHj. 
However,  each  of  these  radicals  would  be  ex¬ 
pected  to  reform  V.A,  decompose  into  products 
other  than  two  acetylenes  or  invf)lve  energetic 
intermediates.  Thus  it  ajjpears  that  the  decom¬ 
position  chain  involving  ti-C4H-,  v\ill  domin.iie. 
Detailetl  chemical  kinetic  calc  illations  using  the 
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ri;a(;  i  k)\  kixl  i  k.s 


spoi  ilic  l  iitc  ( (iiisMilis  111  l  iihle  I  iii'o  <  <>ii.stsieiii 
u iili  iliis  .iiiaK SIN. 

I’riKliu'tioii  oi  (liaciMvlciif  is  dcsci  ihecl  l>\ 

R  +  C.H,  —  +  RH 

R  +  CiUj  —  n-(;,ll:,  +  Rll 

Inllowcd  lj\  dcioinposition  (i|  (dUi  into  <lia- 
ii-uleiH-  pills  H-.itoins.  R  ina\  lx*  an  H-atoin  or 
a  indrotai  1)011  r.idiial. 

AsslioiMi  in  Fij'.  a,  reasonalile  ai^i  eeineiu  will) 
expt'i  imi'iiial  jiroHlcs  ol  ai  oinaiic  species  were 
oliiained  iisini;  die  re.iciion  sei)uence  of  Tal)le  I. 
f  ile  appareiii  underprediction  of  \'A  at  ele- 
\aied  teiiiperainres  is  piesuiiialiK  an  experi¬ 
mental  problem  due  to  partial  gas  sampling  of 
houndarv  lasers  uliicli  contain  the  unlieated 
parent  li\ drocai bon.  Ahose  l(i()()°K.  |)lienvl  is 
lornied  Iroiii  the  n-(.\,f  {t  (ijHj  i ecombination 
and  is  lollimed  b\  com  ei  sion  to  ben/eneor  ])he- 
in laceiv lene.  At  louei  teni|)eratures.  the  |)re- 
doiiiinate  lonii.ition  routes  o(  ttroniiitics  ;ue 

R-la  C.dl,  +  —  (:„H: 

1 1  .:h.')-hexatrien-:<-\  I) 

(;„H:  —  c-(;„H: 
c-(:„H:  (;„H„  ^  H 

loi  hen/etie  lorniaiion.  tiiiil 

R.VJ  CdHi  -  —  Cslh. 

—  c-c;jH., 

c-CJh,  —  CJK  -  11 

to  produie  si\ reiie.  I’A  is  pioduied  principalh  b\ 

R.'i.'i  C-H,  *  M  ->  C-H:  -  H  . 

C:-H:  —  C-H,,  *  II 

near  1  alKl'K  .md  helms,  bui  at  higher  lenipera- 
lures.  aietslene  addition  to  pheiisl  alsocxcurs. 
Re.utioiis  -i.*).  .'I'J.  and  s\eie  assumed  to  be 
oserall  proiesses  and  iioiii  esei  sible.  For  die 
I  Olid II  ions  ( radii  al  <  OIK  em  I  at  ions  and  tempera¬ 
tures)  at  sihicli  these  le.iciions  comrilxiie.  pre- 
liniinars  i  ali  illations  i  onlirmed  the  salidiis  o| 
these  assuniptions.  Despite  the  good  agreemem 
M  ill)  I  he  ex))ei  inienial  data.  I  here  are  sesei  al  un- 
seiding  lealuiesol  these  pi  oposals.  Reaition  d.'i 
reijuires  a  1  or  I  .-1  H-aioni  shift  prior  to  cscli- 
/.iiinii  and  Reaction  -17.  whiili  has  been  pro- 
posi'd  presiousK’  should  he  .i  preferred  route 
(ill  at oinaii/.iiioii,  Hovi es ei .  =  li)"cc'mole- 

sei  i\  ould  be  i  et|uired  loexplainiheexiterimeii- 
lal  data.  1  his  I  ate  is  signilicaiiils  higher  than  the 
late  in  Fable  1.  L’ x  1 1)"(  t.  niole-sei .  rei|uired  bs 
the  ai  els  lene  modeling,  and  pres  ions  determi¬ 
nations  (T  .'ix  in'"  (Kel.  .'))  and  '1.1  x  in"(Re(.  li) 
at  l:ln(FKl, 


Another  concern  of  the  niodeliiig  results  is 
that  the  saliie  of  k-,-..  reijuired  to  match  the 
experittieiiial  dtita  is  a])proxiiiiaiels  100  times 
higher  than  a  presious  determimition  using 
data  from  a  low  pressure,  premixetl  huttidiene 
flattie.  Rate  deteriiiinations  from  the  SP.ST  dtita 
are  subject  to  errors  due  to  boundtirs-ltisers, 
i|uenching  effects,  tmd  modeling  coniple.xits'. 
Flosveser,  the  flame  determinttiion  is  expected 
to  be  a  loner  limit  since,  schen  esiiluating  the 
rate  constaiu.  Cole,  et  al.’  assumed  that  the 
n-CjH-,  tadical  seas  the  dominant  C^H^  species. 
Other  isomers,  pariicularls  FljC :  CHCCHj. 
sshich  is  less  rettciise  and  more  stable  bs  .S 
kcal/inole,  can  contribute  sigiiificatuly  to  mtiss 
5:f. 

The  iissuiiiption  of  the  predominance  of 
n-CdH;,  in  the  flame  work  seould  onls  ptn  ti.ills 
explain  the  difference  betseeen  the  iseo  estilua- 
tions.  Lncertainties  in  heats  of  formation  and 
entropies  for  n-C^H-,  and/or  sinsl,  nhich  mrs 
affect  the  n-(i^H-,  concentration  sia 

C.H,.  -  C.H:,==n-CiH:, 

tints  also  help  to  allesiate  the  difference;  but 
prolntbls  bs  no  more  than  ti  factor  of  fise. 
••Uterntitisely,  reactions  not  considered  in  the 
present  svork  mas  describe  production  of 
iiromaiics  or  significanils  enhance  the  concen¬ 
tration  of  n-C^H•,. 


Conclusions 

Sitigle-pulse  shock  tube  cUtta  base  been 
obt. lined  lor  the  psrolsses  of  aceislene  tmd 
s  ins  laceis  lene.  Flie  dtitti  base  been  used  to 
support  inesious  jtroposnls  for  acetslene  tmd 
ben/ene  jisrolssis  sciili  some  resisions.  de¬ 
composition  model  lor  s  ins  laieis  lene  insols  ing 
H-;iioni  addition  to  s  ins  Ittcets  lene  h.is  been 
|)ioposed  in  oilier  to  ex))liiiii  produit  fotiii.i- 
tioii.  Sesere  constraints  imposed  bs  sinuiliane- 
ousls  modeling  decomposition  of  eticli  of  the 
patent  hs drocai  bons  as  ssell  tis  fise  or  more 
pmiluct  species  httse  been  satisfied  ssith  the 
pi  oposed  model.  \’ins  Ittcets  lene.  Iien/ene.  tmd 
phens  liu  ets  lene  are  obsei  s  ed  during  acets  lene 
ps  rob  sis.  'File  se(|uence 

Cdf.  -I-  FI  ^  (AH. 

C;,ff,.  -r  CAH,  ;=t  ii-ChHs 
n-CoH,j:i(oH4  +  H 

+  n-Cdls  —  f-Ci.H: 

( 1  ..'b.'i-hexati  ien- 1  -s  1) 

(-(;.,H:  —  c-(;„H- 
c-(.„H:  (;„F1„  -  H 
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ACETYLENE  AND  VINYLACETVLENE  PYROLYSIS 


describes  the  low  icmpcraiurc  formation  of 
vinylacetylciic  and  benzene  (1  100— 1400°K). 
Reverse  rates  calculated  from  thermodynamics 
for  the  first  three  steps  describe  the  decomposi¬ 
tion  of  vinylacetylcne.  At  temperatures  above 
1500°K,  phenyl  is  formed  by 

C2H2  +  H  si  C.H  +  Ho 
CoHo  +  CoH  Si  n-QHj 
CaHo  +  n-C^H,  si  f-CfiHj 
( I  -hexyne-3,5-dien-6-yl) 
f-CfiHiSiphenyl 

and  is  followed  by  formation  of  benzene  or 
phenylacetylene.  Reverse  rates  calculated  from 
thermodynamics  are  consistent  with  a  previous 
proposal  for  decomposition  of  phenyl  during 
benzene  pyrolysis. 

Aromatic  formation  during  vinylacetylcne 
pyrolysis  at  low  temperatures  is  described  by 

C.Hj  +  C4H4  CeH: 

( 1 .3,5-hexatrien-3-yl) 

GiH;  —  C-CfiH; 

c-C,iH;  -  CoHo  +  H 
and  n-C^H,  +  ->• 

—  C-CsHr, 

C-C^Hg  —  ChHh  +  H 

The  proposed  rate  constant  for  the  first  step  is 
low,  but  the  reaction  requires  an  H-atom  shift, 
'fhe  retjuired  rate  for  the  addition  of  the  normal 
butandienyl  radical  to  vinylacetylcne  is  approxi¬ 
mately  100  times  higher  than  a  previous  mea¬ 
surement.  Consequently,  it  is  possible  that  calcu¬ 
lated  concentrations  of  n-CgHs  are  inaccurate  or 
that  other  mechanisms  are  significant. 
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COMMENTS 


F.  Temps,  MPI  F.  Stromungsforschung,  West  Germany. 
Would  the  author  please  comment  on  the  possible 
role  of  vinylidene  (singlet  or  triplet)  or  triplet  acety¬ 
lene  on  the  acetylene  pyrolysis? 

Author’s  Reply.  Vinylidene  appears  to  be  a  strong 
candidate  for  initiation  of  acetylene,  since  it  lies  only 
about  44  kcal/mole  above  the  ground  state  of  acety¬ 
lene  and  since  the  intermediate  transition  state  is  only 
two  (2)  to  four  (4)  kcal/mole  above  vinylidene'.  A 
subsequent  step: 

HjC  -  C:  -t-  HCCH  —  H,C(:CHCH  —  i-C,H,  -f  H 

is  endothermic  (overall)  by  about  15  kcal/mole  (or  69 
kcal/mole  above  two  acetylene  molecules).  The  overall 
endothcrmicity  of  this  sequence  is  lower  than  oiherpro- 
f)OScd  bimolecular  initiators  which  formCjH  +  C7Hjor 
n-C4Hj  +  H,  which  arc  about  86  and  79  kcal  endother¬ 
mic,  respectively.  Unfortunately,  the  endothermicityof 


the  sequence  involving  vinviidene  is  still  about  13  kcal 
higher  than  the  activation  energy  of  the  bimolecular  in¬ 
itiator  (Reaction  7)  used  in  the  present  work. 

Vinylidene  could  also  play  a  role  if  the  above  inter¬ 
mediate.  HjCCCHCH,  undergoes  a  1,4  H-atom  shift 
to  produce  vinylacetylene.  Although  this  route  was 
considered  in  the  present  work,  there  is  no  apparent 
parallel  pathway  to  benzene,  which  adequate!)  de¬ 
scribes  experimental  profiles.  It  is.  of  course,  possible 
that  this  mechanism  contributes  to  the  production  of 
vinylacetylene.  A  possible  role  of  triplet  acetylene  w  as 
not  considered  in  the  above  analysis. 
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ACETVLLN't:  AND  VINYLACETVLENE  I’VKOL^  SIS 


S.E.  Sleiii,  Cliewiciil  Kinetics  Dw,  \aliiinal  Hurrau  oj 
Sinntlards,  L'.S.A.  How  do  you  intei  prei  ihe  unusualh 
higl)  rate  constant  for  styrene  formation- 

Author’s  Kellis.  As  simnnarized  in  the  conclusions  of 
the  manuscript,  the  rate  expression  htj  -  7.9x  10’’ 
exp(-3000/RT)  cc/mole-scc,  is  about  a  factor  of  100 
too  liigh,  when  compared  to  rate  constants  of  compa¬ 
rable  addition  processes.  The  interpretation  is  that 
either  other  mechanisms  dominate  styrene  formation 
or  that  rate  coefficients  and/or  thermodynamics  for 
several  of  the  reactions/species  are  in  error.  The  most 
likely  alternative  mechanism  is  a  dimerization  pro¬ 
cess;  the  most  likely  errors  associated  with  rate 
coefficients  and  thermodynamics  probabli  involve 
vinyl  and  butadienyl  radicals. 

There  is  presently  insufficient  information  avail¬ 
able  to  ascertain  which  of  the  aboie  (or  if  all) 
contribute  to  the  discrepancy.  Consequent!) .  the  data 
are  presented  as  is. 


K.fl.  Homcinn,  Technische  Hochschule  Darmstadt,  UVjt 
Germany.  Did  you  consider  a  dimerization  step  as  the 
first  reaction  in  the  pyrolysis  of  CjHj  and  C(H<-  In  a 
recent  study  of  CjHj  pvrolysis  at  lOO'C  we  found  that 
the  homogeneous  initiation  step  was  2  C^H2  — •  C»H^. 
Addition  of  H  atoms  to  the  reaction  ssstem  had  no 
effect  on  the  products'. 

1.  Homan.x  K.H.  .a.vd  Pidoll  Ber.  Bun- 

senges.  Phys.  Chem.  1986.  in  press. 

Author's  Reply.  In  an  early  analysis  of  this  work,  the 
dimerization  of  acetslene  was  considered  as  the 
predominant  mechanism  for  vinylacetylene  forma¬ 
tion.  The  rate  constant  {k  =  3.6=  10'^  expI-J*!  lOO/RT) 
cc/mole-sec)  reported  b>  Bradlei  and  Kisiiakowskv' 
was  found  to  describe  nicely  the  low  temperature 
formation  of  C,H4.  Coincidentally,  the  overall  activa¬ 
tion  energv  of  this  process  is  similar  to  the  endother- 
micitv  for  the  formation  of  HC^CHCHCH.  the  diradi¬ 
cal  intermediate.  The  endothermicitv  of  this  process 
was  estimated  to  be  3.3  to  *13  kcal/mole  bv  interpreting 
results  from  Kollmar.  ei  af.  Despite  the  good 
agreement  with  experiment  and  reasonable  mecha¬ 
nistic  explanations  (see  also  the  response  to  Prof.  F. 
Temps  regarding  the  role  of  vinviidenc)  several 
problems  arose.  First,  the  Uradlev  and  Kistiakowski 
data  were  sh'.ivvn'  to  be  misinterpreted,  leav  ing  vtoubi 
regarding  their  reported  rate  constants.  Secondlv. 
theie  appears  to  be  no  parallel  paihvvav  to  the 
formation  of  benzene  at  low  temperature  vvhich 
describes  the  experimental  results.  In  support  of  the 
chain  mechanism  during  acetylene  pyroivsis  is  that 
production  of  (hH4  and  C^Hr,  can  be  explained  using 


the  proposetl  chain  which  includes  rate  constants  very 
much  consistent  with  previously  proposed  rate  con¬ 
stants  and  thermochemistrv .  .Additionally,  many 
other  products  were  not  reported  due  to  their 
relatively  low  concentrations;  but  their  existence  is 
strongly  indicative  of  radical,  chain  mechanism.  Fur¬ 
ther  support  for  the  chain  mechanism  is  drawn  from 
the  fact  that  the  proposed  mechanism  qualitatively 
predicts  the  product  distributions  observed  during 
lovv  temperature  pyrolysis  of  acetylene^  (allowing  for 
a  small  change  in  the  initiation  rate  and  for  pressure 
dependence  of  addition  reactions). 

Finally,  Callear  and  Smith'’  have  shown  convinc¬ 
ingly  that  processes  very  similar  to  those  described  in 
this  work  occur  at  room  temperature  when  H-aioms 
are  added  to  acetylene.  The  conclusion  of  the  present 
work  is  that  a  radical,  chain  mechanism  adequately 
describes  the  observed  dominant  products  using 
mechanisms  and  rates  which  are  reasonable.  Dimeri¬ 
zation  processes  cannot  as  yet  be  ruled  out  absolutely 
(and  they  mav  plav  a  contributing  role),  but  it  does 
not  seem  necessary  to  invoke  their  existence  to 
explain  the  observed  results. 

-As  strong  a  case  against  dimerization  processes 
cannot  be  made  in  the  case  of  vinylacetylene.  since  it 
is  necessary  to  add  a  reaction  (#45)  vvhich  had  not 
been  proposed  prev  iously  and  (most  disconcerting)  to 
use  an  unusually  high  rate  for  the  addition  of  n-C4Hs 
to  C,H4  (R52)(see  response  to  question  by  Dr.  S. 
Stein). ^  However,  with  the  lack  of  evidence  of 
dimerization  processes  occurring  during  high  tem¬ 
perature  acetviene  processes,  the  author  favors  a 
chain  mechanism  for  explaining  the  product  distribu¬ 
tion.  The  vinylacetylene  data  can  be  fitted  reasonably 
if  one  assumes  dimerization  processes  dominate 
aromatic  formation.  For  example,  by  reducing  Ag,  *45, 
kyi.  and  *ij  by  factors  of  2.  2,  100,  and  2  respectively 
and  including: 

2C4H4  —  CsH,(a)  and  2C4H4  CsHs  -r  2H  (b) 

(as  irreversible  steps)  with  respective  rate  constants 
of: 

k.  -  9x  10”  exp(-40000/RT) 

and 

-  3x  10'^  exp(-35000/RT)  cc/mole-sec 

then  experimental  profiles  can  be  predicted  reasona- 
blv  (although  benzene  is  overpredicted  at  low  tem¬ 
peratures).  k..  in  comparison,  is  approximately  two  to 
three  times  higher  than  an  extrapolation  of  Lund- 
gard  and  Heicklcn's(6)  expression: 

*•  -  3x10"  exp(- 29200  cahmolc/Rl )  cc/mole/sec 

for  the  net  formation  of  “dimerized"  products, 
formetl  via  a  combination  of  molecular  and  radical 
processes.  3  he  radical  process  at  low  temperatures 
was  assumed  to  involve  formation  of  a  diradical 
intermediate,  although  possible  contributions  from 
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the  chain  process  as  suggested  in  this  paper  sliotild 
not  be  necessarily  ignored. 

REFERENCES 

1.  Bradley.  |.N.  and  Ki.stiakowskv,  G.B.;  J. 
Client.  Fhys.  35,  2(id  (I9fil). 

2.  Koll.mar.  H..  Carrion,  F.,  Dewar,  M.J.S.. 
Binoma.m,  R.C.;  J.  Am.  Cheni.  Soc.  103,  5292 
(1981). 

3.  Gay,  I.D..  Kistiakowsky,  G.B.,  MtciiAEL,  J.V. 
and  Niki.  FI.:  J.  Chem.  Phys.  43,  1720  (1965). 

4.  Ref.  #30  of  manuscript. 

5.  Callear.  A.B.  and  Smith,  G.B.:  J.  Phys.  Chem. 
90.  3229  (1986). 

6.  Lundcard.  R.  and  Heicklen,  J.:  Int.  J.  Chem. 
Kin.  16,  125  (1984). 

• 

S.J.  Hams,  General  Molars  Research  Labs,  U.S.A.  Are 
C.iHj  and  CjFlj  equilibrated  in  your  system?  At  what 
temperatures?  If  they  are  equilibrated,  then  you 
cannot  obtain  mechanistic  information  from  the 
coniractions.  If  they  are  not  equilibrated,  are  they  far 
enough  away  from  equilibrium  so  that  you  have 
confidence  that  the  mecfianisttt  is  controlling  the 
products? 

Author's  Reply.  For  the  dwell  times  (500-800  micro¬ 
seconds)  of  these  pyrolysis  experiments,  the  acety¬ 
lenes  are  equilibrated  only  aboYC  1500-1700  K. 
depending  on  the  concentration  and  identity  of  the 
reactant.  Major  coiicltisioiis  of  the  present  work  are 
dependent  strongly  on  the  experimental  restilt  that 
products  of  low  temperature  ( 1 1 00- 1 300  K)  pyrolysis 
of  vinylacctvlene  are  principally  acetylene,  diacety- 
lenc.  and  hvdrogen  with  the  product  ratios  C;Hj/ 
C<M'j  •'•nd  C..H.../H..  approximately  ten  to  one.  Pertur¬ 
bations  of  ihese  priKluci  ratios  By  e(|uilibraiion 
reactions  would  sigiiificaiiily  alter  the  concittsions  of 
this  pa))ci  .  The  "non-e<|uilibratioii"  of  acetslene  at 
low  icniperatures  and  shoii  limes  c;iit  be  ilemon- 
slraied  bv  the  following  iiKKleliiig  studv  and  bv 
contpaiistm  to  ex|KM  imenls  on  ben/.cne  pyrolysis. 

In  an  artificial  modeling  exercise  ol  (ij  1 1  pvroivsis, 
rale  coefficients  for  Rcaciioiis  4(i  and  20  wcie 
decreased  while  coefficients  of  Reactions  13  and  18 
were  increased  be  factors  often  (10)  or  more.  Tlie  net 
lesitli  was  to  convert  the  overall  det (iin()osiiion 
palhwai  from. 

G,ll,  -*  2C..H  , 
to 

—  (.,11,  +  If, 

Fhe  iiKKlelmg  results  below  140(1  K  indiiated  final 
G,Hy(.,/f,  ralio.s  of  appioxiinaielv  one  hiindieil 


(100)  to  one  which  is  ‘2'A  orders  of  niagnitudc  above 
ecpiilibrium.  Gon.scc|ucntly,  this  modeling  exercise 
indicates  the  relative  uniinportance  of  equilibrating 
reactions  under  the  low  temperature  conditions  at 
which  the  dccom[M)sition  mechanism  was  cfeduccd. 

F?xpcrimciital  evidence  for  non-equilibration  at 
low  temperatures  can  be  inferred  from  recent  data  on 
benzene  pyrolysis.  At  )500K.  several  hundred  de¬ 
grees  higher  than  the  temperatures  at  which  mecha¬ 
nistic  information  was  obtained  in  the  present  work, 
Colket'  found  an  initial  C2Fl2/C,H2  ratio  of  approxi¬ 
mately  three  (3)  to  one  (1),  a  value  lower  than  the 
equilibrated  ratio  of  alxiut  five  (5)  to  one  ( 1 )  (depend¬ 
ing  on  assumed  thermochemistry).  The  experimental 
ratio  is  similar  to  the  value  (two  to  one)  obtained  by 
Kern,  et  al.^,  at  yet  higher  temperatures  (1704  K). 
Consequently,  both  modeling  and  experimental  data 
support  arguments  that  earlv  (i.e.  initial)  foriiiaiion  of 
acetylenic  products  are  not  equilibrated. 
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• 

J  H.  Kiefer,  L'liiv.  of  lllitiot.'.  at  Ghirago,  L  .S  A.  In 
answer  to  Dr.  Harris'  question,  current  models  ol 
high-ienipcralure  acetylene  (syrolysis  require  (i,H 
radicals  for  C,H,  equilibration,  and  the  formation  ol 
this  radical  seems  unlikely  at  these  temperatures.  1 
also  have  questions  for  the  speaker.  First,  for  \otir 
chain  production  of  C,!!.,  to  proceed  sufficiently  1 
suspect  you  may  well  need  the  high  rate  of  vine  I 
dissociation  you  assumed.  Otherwise  the  fast  process 
C.,H3-eH  — *■  C2H2-(-H-,  would  very  effectitely  termi¬ 
nate  the  chain.  Do  voti  think  voti  cotild  stand  the  low 
rate  of  vinyl  di.s.socialion  of  votir  butadiene  woik'- 
Finally,  if  the  heat  of  formation  of  C,H  ,  were  alHiui 
130  kcal/inol.  as  has  Iteen  stiggested  to  me  b\  Kail 
Melius,  could  your  model  still  provide  sulluieiu 
conversion  of  C,H,? 

RKFF.RKNCE 

1,  Kiekek,  ).K..  Wei.  H.(f.  Kern  R.D..  and  Wt  , 
(Ml.:  Int.  |.  (ihem.Kin.,  17.  22;')  (1985). 

Author's  Rrpl\.  Fhe  pliilosopin  lor  selection  c.l 
reactions  and  coi  lespondiiig  rate  coiistaiils  pi  eseiited 
in  Table  I  was  to  tise  well-established  mechanisms  and 
rate  constants  while  selecting  other  mechanisms  and 
rate  coitslants  as  re(|uiied  by  the  experimental  le- 
sults.  but  consistent  (in  most  c.ises)  with  thei  nioche- 
mistry.  Vinvl  detonqiosition  w.is  assutned  to  be 
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pressure  iiiilependeni  ;in(l  ils  rule  (iinsuitil  uiis  iiiii 
varied  in  the  aiiaivsis  desei  il)erl  by  the  iuaniiseri|ii.  A 
preliiuitiarv  exaiiiiiialiou  ol  the  elleei  ol  the  use  ol  a 
lower  rale  loi  \iiivl  deeoiiiposilion  was  exainmerl  h\ 
usiiin  Kieler’s  tale  expression  wliieli  is  more  than  tutl 
times  slowei  than  that  used  in  ih  present  anaivsis.  In 
the  modeling  ol  vinvlaeelylene.  the  oveiall  rler  onipo- 
silion  rate  ol  the  parent  deercased  hv  onlv  l^tl'.;. 
although  the  "tpiasi-sleady"  eoneentralion  ol  vinvl 
radicals  inei  eased  hv  neatly  a  lattor  ol  100.  1  he 
overall  routeol  vinviaeetylenedeeoinpositionehanged 
onlv  slighih  Irom  that  de.serila-d  in  the  manusrripi. 
.Main  pit 'tleins  that  at  ise  ate:  (I)  diacetviene  and 
ethene  are  overpredicted  hy  ahout  a  laetor  ol  two  aiul 
five,  respetiivelv;  and  (2)  early  formation  of  laiih 
vinviacelvlene  and  benzene  froin  acetviene  pvroKsis 
arc  underpredicied  each  hv  alanit  an  order  ol 
magnitude.  Some  features  suggest  that  an  intermedi¬ 
ate  value  lor  vinvl  deioinposilion  might  he  leason- 
ahle.  file  d.ita  on  vinviacelvlene  ;in<l  ateivk  ne  Iroin 
this  work,  however,  do  not  support  the  low  rale 
constant  as  tleicrmined  in  your  anaivsis  ol  inv 
butadiene  data. 

If  (i(Hi  radicals  have  a  heat  ol  formation  ol  l.'fo 
kcal/mole.  it  would  tend  to  support  present  arguments 
that  ll-aloms  add  to  vinviacelvlene  laiher  than  ab¬ 
stract  M-aloms  from  vinviacelvlene.  Hv  usinglhis  high 
heal  of  (drmalion  in  modeling,  however,  two  major 
problems  will  arise,  first,  the  elfeciive  activation 
energv  ol  the  initiation  process  ((i.l I ,  — •  i-( i,! 1 1  Hi 
was  S.'i  kcal/mole  in  this  work,  i  his  low  v.due  is 
(pieslionahle  as  it  stands,  hut  would  he  nearlv  impossi¬ 
ble  to  reconcile  il  the  overall  endolhermicli  v  were  I  I  :f 
kcal/mole.  .-\  second  prohlem  results  Irom  the  laci  that 
the  high  heal  of  formation  (vs.  12*)  kcal/mole  assumed 
in  this  work)  would  suhslanliallv  reduce  the  '(piasi- 
sleadv  slate’  conceiiiralion  of  n-Cdlli.  finis,  the 
lor  111.11  ion  ol  phenv  I  v  ia  n-(i|l  I ,  addition  to  .k  elv  lene 
would  Ik'  suhst.mti.div  reduced  and  heii/ene  .md 
]>henv  l.iceivlene  signiluanllv  tinder ))redu ted 

• 

Iff.'  (ffintinrt.  nf  I  f  liuhisioii 

of  il  revel  sihle  steps  m  |>olvmei  loi  in.iiion  met  h.i- 
nisms  has  to  Ih-  reg.iitlerl  with  spei  i.il  <.ne.  p.iriiin- 
l.n  Iv  III  the  i  .ise  of  soot  I  oi  illation,  whet  e  the  divei  siiv 
ol  p.iihw.ivs  ,md  uiisiispeded  leversihililv  ol  ipiiie 
exoergic  steps  make  the  points  ol  no  lelnrn  hist 
appear  at  moleciil.n  sizes  well  hevoml  those  loii- 
sidereil  in  this  p.ipei '. 

In  oidei  to  make  v.ilid  coinpai  isons  heiween  com¬ 
puted  .md  me.ismed  vields  oi  one-img  ,iiid  iwo-img 
arom.ilu  (ompomids.  the  le.uiion  mec  h.mism  miisi  he 
extended — under  provision  lor  reverse  le.iilioii — lo 
slill  larger  species,  .is  a  miniimim  lo  .ii  en.iplh.dene. 
for  this  re.ison  it  would  appe.n  ih.ii  the  i  .iie  i  onsi.mis 
used  in  this  p.ipei  to  model  lorm.iiioii  ol  .nom.iiii 
rings  h.ive  to  Ih'  leg.nded  asaiiiliii.il 
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I.  fKKNKi.vc  11.  M..  Ci.vKV.  D.W'..  (ivKMxnc.  \\  (f. 
.nid  SiKiN.  .S.f..:  ‘JDih  .Svmp  (Ini.)  on  ('omliiis- 
lion.  p.  KS7.  file  (  omhustion  Instiiute.  1‘ilts 
burgh.  lO.S.a. 

Aiithn}\  lirj>l\.  In  gener.il.  vom  lommeni  is  veiv 
reason. ihle  and  should  he  considered  lor  long  re.ii- 
lion  limes  or  conditions  whiih  .illow  conversion  ol 
signilicani  mailer  lo  l.irgei  species.  Ilowevei.  lor  the 
conditions  at  which  the  modeling  w.is  .ipplied  in  this 
wink,  this  elleci  is  rel.ilivelv  nnnn|)oi  lani.  In  the 
present  siudv.  the  kinetics  an.ilvsis  w.is  peiloiiiud 
for  short  reaction  times  (701)  mii  i  osei  onds)  and  low 
leinperaiuies  (<l(‘)()l)  K)  lor  wliiih  eai  lv  loi  ni.iiion 
kinetics  (ol  arom.ilirs)  could  he  observed.  I  ndei 
these  conditions,  rel.ilivelv  sin. ill  Iraiiionsol  the  ob¬ 
served  aiomalics  inidergo  growili  piiHCsses.  l  oi  ex¬ 
ample.  in  the  case  ol  :(..*>'/  .nelvlene  pviolvsic.  ID 
ppm  ol  benzene  w.is  ohserved  .il  I  I0()  K.  the  lowesi 
lemperaiure  .it  whnh  the  mixtuie  vvas  |)violvzcil 
(amcenttaiions  of  phenvi.ueivlene  .it  this  lempei.i- 
lure  can  l>e  estimated  to  he  I  W  to  2  orders  ol  m.igni- 
lude  lower  and  do  noi  .ii  hieve  a  comp.ir.ihle  lotueii- 
lialion  until  l.'iOO  K  lloi  the  s.iinc' dwell  lime).  H.ised 
on  data  using  higher  mili.il  lonienli.iiions  ol  .neiv- 
Icuc  (-f.fl'/c).  ct.cpluU.ileive  coin  emi  alions  .ne  mote 
ih.in  an  order  ol  niagniinde  below  ihoseol  phenv  l.i- 
celvlene  at  these  lemper.iiiires  An  .in.ilvsisol  ni.iss 
balance  indicates  good  leioveiv  ol  initi.il  c.irhon  he- 
lovv  |l‘>*>l)  K.  I  hese  resulis  m.iv  not  he  proof  hiil  iliev 
indicate  sironglv  ili.il  the  jn  oi  esses  observed  in  this 
siudv  are  desc  l  ipliv  e  ol  ihe  nnli.il  loi  m.ilion  su  ps  lol 
arom.ilics)  and  |>eriui  h.iiion  due  lo  kick  ol  inthisioii 
ol  siihsecpieni  growth  proi  esses  aie  negligible. 

• 

W-  I'trfifihitit.  Sliilr  I'liiw,  (  S  A  ()ni 

iiiei  haiiisin  ll.2l  lo  vsliuli  the  .nillioi  leleis.  v\.is 
c  om|M>sed  ol  .1  \  ei  V  l.n  ge  v  .n  let  v  ol  i  e.ii  i  ion  p.il  liw  .iv  s 
iiicliidiiig  .ill  those  distiissed  hv  tin  .iinlioi  No 
.issumpiioMs  .ihotil  the  nnpoii.nue  ol  niilividii.il 
elemeni.iiv  le.uiioiis  in  the  mei  h.imsin  v\eie  m.ide: 
oiii  I  one  liisions  .ihoni  donnn.ini  u.itiion  p.iihw.ivs 
le.idiiig  lo  polvc  V I  III  .11  OIII.  1 1 II  hv  di  Ol  .11  lions  .nid  sooi 
vveie  di.ivvn  nisie.id  lioin  die  lesnlis  ol  lonipniei 
expel  imeiils  desi  I  ihed  in  llieiiled  leleieiiies  I  Inis, 
loiiii.iiv  to  die  .niihoiss  .isseinon.  die  p.idiw.iv 

11  -  (  ,11.  —  (  ill,  •  II 
(:,ll,  *  II  ^  n  -  (  ,11,  ♦  If  ih 

W.IS  iiol  fivvio/iei/  hill  i/o/iioif  lo  he  the  donnti.nil 
p.ilhw.iv  loi  die  piodnition  ol  n  -  (  ,11,  loi  die 
1 1 II II  III  Mil  IS  ol  "soul  "  lol  1 1  l.n  ion  III  shot  k-t  nlie  pv  i  oh  sis 
III  .iietvlene.  .-Mdioiigli  we  lontid  die  le.iilion 

(.,11,  *  (  ,11  —  11  (  ,11,  r.'i 
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u>  l)C  l;iM  ill  vei  l  sliorl  rciiciioii  liiiicb,  il  t|iiii  kly  g.ivc- 
way  lo  llic  iCiiUL'iitc  (1)  tliiii  liciiimc  (loiiiitiani  loi 
soot  appearance.  I  lie  compelilion  between  (1)  iiiul 
(1?)  is  sensitive  lo  experinientitl  conditions;  e.g., 
reaction  (2)  plays  a  more  prominent  role  in  ilie 
presence  ol  <)xyj;cn  Hence,  dillerences  in  experi- 
meniiil  conditions  and  dilleieni  assumptions  alxiiii 
lliermixliemical  values  ;ire  the  likeK  sources  ol  the 
diflerenl  relative  lliixes  of  (1)  and  (2)  noted  in  this 
paper. 
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Aiilliiir's  Reply.  As  siij'gested  b\  Dr,  Fienklacli.  the 
manuscript  should  be  chant;ed  to  l  ead;  "Previouslv,  i; 
\\, Mill'd  need  .  .  .  thill  n-(!-l  H5  wiis  formed  ihrou>;h  path 
A,  whereas  (R-12)  w;is  the  domiii.mi  route  under  the 
present  conditioiis."  In  boili  the  present  and  previous 
work,  it  should  be  remembered  ih.it  the  authors  have 
l>cen  forced  to  esliniaie — or  assiinie — values  of  abso¬ 
lute  rate  consianis  for  which  either  no  data  or  little 
data  are  available.  Oftentimes,  siibsliiiuial  tinceriaintv 
exists  in  estimation  lechminies  C.onseijuently.  ibis 


piocediiie,  bv  iis  very  nature,  leiids  lo  n.'isumpliim.s 
aboui  the  imporiance  of  iiitlividiiiil  elemeniarv  reac¬ 
tions.  (.'oncliisions  reached  or  dediued  in  aiiv  work 
usually  ilie  onlv  iis  f;ood  as  the  iniliiil  iissiimpiions 

• 

A. A)  Demi.  Exxuti  Res  cs”  Co..  L'.S..\.  Civeii  the 
last  rale  lonsi.mi  for  C-.dl,  +  C<H,  (with  no  aciiviiiion 
energy)  ihiil  w.is  needed  lo  explain  C,.li..  lormaiion,  is 
il  not  |)ossible  lo  conclude  that  the  rapid  latc  ol 
aromatics  lormaiion  is  siill  a  puzzler 

.Aul/iiii'i  Reph.  Flic  rale  constant  lor  k,s  was 
determined  to  be  -4X10"  cc/mole-scc,  which  is  iilxiut 
two  orders  ol  magnitude  below  collision  Irctiiicncv 
and  is  similar  m  magnitude  to  rate  constants  deter¬ 
mined  lor  other  addition  reactions  The  onlv  unusual 
feature  regarding  Reaction  45  is  that  the  linear 
intermediate.  1 .8..5-hcxatricn-3-vl,  must  undergo  an 
H-iiiom  shill  and  convert  to  the  1 .3.5-hcxaii ien- 1 -\  1 
radical  prior  lo  cyclizaiion.  It  is  believed,  however, 
that  this  H  aiom  shift  is  not  rate  controlling,  Tbesi/c 
of  the  rate  constant,  in  the  opinion  of  the  authors 
ihercfore.  does  not  appear  to  be  unreasonablv  large. 
File  activation  energv  could  not  be  determined  in  this 
experimeni;  and  so  none  vvas  assigned,  although  it 
would  not  lie  unreasonable  to  assume  an  ellective 
activation  energv  of  at  least  several  kcal/mole. 

Certain  teaiures  ol  aromatics  formation  certainlv 
can  be  consitlered  lo  be  unresolved,  or  at  least  under 
di.scussion,  .See.  lor  example,  the  other  questions  and 
responses  to  this  paper.  Certainlv  more  expei  imental 
and  modeling  work  needs  to  Ire  performed.  One 
critical  area  of  reseaich  is  the  accurate  determinaiion 
of  enthalpies  and  entropies  (and  their  temperature 
ilependencies)  lor  important  radical  intermediates. 
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